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supply has been in a state of flux and under Congressional scrutiny 
Since 1985. This thesis analyzes and compares the U.S. Navy's 
current economic retention process to four mathematical Economic 
Retention Decision Models designed to assist in making retention 
determinations with respect to excess inventories. The motivation 
for this research was based on several factors, the two primary 
factors were; the Navy does not currently use a classical economic 
retention decision model when making retention/disposal decisions 
for "essential" material, and U.S. Navy inventories in long supply 
were estimated to be as high as 3.4 billion dollars in March 1993. 
A Pascal based simulation was developed to compare the Navy's 
retention process and the mathematical models. The comparison was 
based on performance with respect to the Measures Of Effectiveness 
(MOE) of Total Cost and Average Customer Wait Time. The simulation 
was designed to emulate the portions of the Navy's consumable item 
inventory management system (UICP) applicable to the demand process 
for a Navy managed consumable item. The goal of this research was 
to determine how effective the Navy's retention process was as 
compared with economic retention decision models for both a steady 
State and a declining demand environment. In general, results 
Showed that at least one mathematical model performed better than 
the Navy's process for all demand scenarios that were simulated and 
that the ideal model varies between demand scenarios and changes in 
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EXECUTIVE SUMMARY 


OVERVIEW: Retention and disposal policy for U. S. Navy 
wholesale inventories in long supply has been ina state of 
flux and under congressional scrutiny since 1985. Comments 
from the Chief of the Supply Corps on 19 July 1993 indicated 
that one of the preeminent issues regarding the future of the 
Supply Corps was inventory reduction. He stated that 
inventory reduction 1S "a congressionally mandated process and 
a fiscal necessity .... we must continue to aggressively 
pursue inventory reductions in an intelligent manner", and 
that it "demands our immediate and continuous attention." 

An important aspect of inventory reduction is the 
retention/disposal process for excess material. This thesis 
evaluated the effectiveness of the Navy's UICP economic 
retention model. The evaluation was performed by comparing 
several mathematical economic retention models with the Navy's 
existing retention model. 

There were three primary factors that motivated this 
thesis. First, the Navy Inventory Control Points (ICP) are 
not confident that eight years worth of forecasted annual 
demand 1S an appropriate inventory retention level. Second, 


with continued budget reductions and reductions in the size of 


ew 


‘Naval Supply Systems Command, Subject: Naval Supply 
Corps FLASH from the Chief, No. 7-93, 19 July 1993. 


yale 


the Fleet, excess inventories will continue to be a financial 
and administrative burden. For example, as of March 1993 the 
Meeeeneld $1.9 billion in Economic Retention Stock: and $1.5 
bellion in potential excess inventory for 1H, 3H and 7 COG’ 
material. Finally, DOD Regulation 4140.1-R recommends that 
better analysis supporting retention decisions be done through 
the use of economic retention decision models. The Navy does 
not currently use a classical economic retention decision 
model when making retention and disposal decisions for 
"essential" material. 

ANALYSIS: An analysis of the models was performed for a 
variety of demand scenarios in both steady state and declining 
demand situations. The analysis was designed with two 
objectives in mind. The first objective was to determine 
which model(s) were most effective in a demand environment 
Similar to the Navy's stochastic demand environment. The 
second objective was to evaluate how the Navy's retention 
process performed with respect to the mathematical models. 

A discrete event Monte Carlo simulation of the Navy's UICP 
demand process and the mathematical retention models was 
developed to evaluate the performance of the models. The 

‘Economic Retention Stock (ERS) is that material which is 
more economical to hold for future requirements as opposed to 
disposing and reprocuring in the future. 

“Cognizant symbols (COG) are two character alpha-numeric 
codes which identify and designate cognizant inventory 


managers who exercise supply management over a specific 
category of material. 


Te 


Simulation was developed by the author and LT Glenn 
Robillard’, and was designed to emulate the portions of the 
Navy's Uniform Inventory Control Program (UICP) applicable to 
this research. The Simulation represents the demand process 
of a hypothetical Navy managed consumable item. The 
evaluation of the models' performance was based on the 
measures of effectiveness (MOE) of total cost (TC) over a 
specified period of simulation time and average customer wait 
time (ACWT) per requisition for all requisitions which occur 
over a specified period of simulation time. 

The mathematical models chosen for this research were 
based on their applicability to the Navy's excess inventory 
problem and the simulation. The mathematical models chosen 
were Simpson's "Economic Retention Period Formula", Tersine 
and Toelle's simple "Net Benefit" model and present value "Net 
Benefit" model, and the simple "Net Benefit" model modified to 
account for the potential for stockouts associated with Navy 
managed items. 

The analysis and performance comparisons of the models 
were based on MOEs calculated from output data from the 
Simulation for six basic demand scenarios. The demand 
scenarios were based on varying combinations of unit price, 


mean quarterly demand and variance of mean quarterly demand. 


‘LT Robillard is a U.S. Navy Supply Officer =a 
graduate student at the Naval Postgraduate School studying 
Operations Research. 


For each demand scenario four retention scenarios were 
analyzed using the simulation. The four retention scenario 
analyses follow. A Total Cost Analysis was performed to 
determine what the true optimal amount of inventory to hold 
was for a given quantity of initial excess inventory. A 
Constant Demand Analysis was performed to compare the various 
models to the theoretically optimal retention quantity that 
was determined during the Total Cost Analysis. A Declining 
Demand Analysis was performed to compare the models under 
three scenarios of declining mean demand patterns. Finally, 
Sensitivity Analysis was performed for four combinations of 
demand scenarios and declining mean demand patterns. The 
parameters evaluated in the Sensitivity Analysis were 
inventory holding cost rate, obsolescence rate, administrative 
order cost rate and salvage rate. 

CONCLUSION: The findings of this research showed that none 
of the models analyzed consistently yielded the lowest total 
cost and ACWT for all of the demand and retention scenarios 
examined. AS a group, the "net benefit" models performed the 
best and generally performed better than the UICP retention 
model. Additionally, for most demand scenarios in both the 
Constant and Declining Demand Analysis, the decision on which 
model to chose could typically be determined by the MOE of 
Metal cost alone. This was due to the fact that the 
difference between the various models' ACWTs for each demand 


scenario, was generally insignificant. In summary, the above 


findings indicate that for Navy managed items the "optimal" 


retention quantity differs significantly from item to item 


based on variations in mean quarterly demand and unit 


price. 


eal 


er BACKGROUND 


A. INTRODUCTION 

Retention and disposal policy for U. S. Navy wholesale 
inventories in long supply has been ina state of flux and 
under congressional scrutiny since 1985. Comments from the 
Chief of the Supply Corps on 19 July 1993 indicated that one 
of the preeminent issues regarding the future of the Supply 
Corps was Inventory Management/Reduction. He stated that 
inventory reduction is "a congressionally mandated process and 
a fiscal necessity .... we must continue to aggressively 
pursue inventory reductions in an intelligent manner," and 
that it "demands our immediate and continuous attention" [Ref. 
Lalli 

A key aspect of inventory reduction is the process used to 
identify two types of inventories: Economic Retention Stock 
(ERS) and potential excess inventory. ERS (sometimes referred 
to aS Economic Retention Requirement (ERR)) is the portion of 
the inventory above current requirements which is determined 
to be more economical to retain for future use as opposed to 
disposing and reprocuring in the future. The sum of current 
requirements and ERS is called the Retention Level (RL) when 
it 1s defined in terms of years worth of annual demand and is 


called Retention Quantity (RQ) when it is defined in terms of 


the number of units. For this thesis the retention limit will 
generally be expressed in terms of years worth of annual 
demand and referred to as the RL. Potential excess inventory 
1s that portion of material on-hand and on order beyond the 
Ray 

In 1985 the DOD adopted a policy to retain all units of 
any ltem having application to a weapons system in active use 
by any of the U. S. military services [Ref. 2]. This disposal 
moratorium was established as a result of inconsistencies the 
GAO identified in U. S. Air Force economic retention policy. 
In effect, the moratorium eliminated the need for any economic 
retention models. Motivated by new GAO findings in 1988 and 
1990 regarding the growth of DOD secondary inventories [Refs. 
3 & 4], in 1990 the DOD lifted the disposal moratorium [Ref. 
Zila. NAVSUP Instruction 4500.13 [Ref. 5] was subsequently 
1ssued to provide policy on retention of wholesale Navy 
material. The retention limit was set at 20 years worth of 
forecasted annual demand for items that have been stocked in 
the supply system for more than seven years and coded as 
"essential" material. Here "essential" material is defined as 
an item whose failure would result in the loss or severe 
degradation of primary mission capability. As a result of the 
shrinking DOD budgets and continued congressional concern over 
large DOD secondary inventories the retention level for 
wholesale Navy material was further reduced in August 1992 to 


elght years worth of forecasted annual demand [Ref. 6]. 


This thesis contains an analysis and comparison of the 
U. S. Navy's current economic retention process to four 
mathematical/optimization models (Economic Retention Decision 
Models) designed to assist in making retention/disposal 
determinations with respect to excess inventories. The 
motivation for this research was based on three factors. 
famest, the Navy Inventory Control Points (ICP) are not 
confident that eight years worth of forecasted annual demand 
1S an appropriate RL. Second, with the ongoing budget 
reductions and reductions in the size of the Fleet, excess 
inventories will continue to be a financial and administrative 
burden. For example, as of March 1993 the Navy held $1.9 
feeeeton in ERS and $1.5 billion in potential excess inventory 
bere 1H, 3H and 7 COG’ material. Finally, DOD Regulation 
4140.1-R [Ref. 7:p. 4.5] recommends that better analysis 
Smeporting retention decisions be done through the use of 
economic retention decision models. The Navy does not 
currently use a classical economic retention decision model 
when making retention/disposal decisions for "essential" 
material. 

A simulation was developed in the Pascal programming 
language to compare the Navy's retention process and the 
mathematical models. The comparison is based on performance 

‘Cognizant symbols (COG) are two character alpha- 
numeric codes which identify and designate cognizant 


inventory managers who exercise supply management over a 
specific category of material. 


with respect to the measures of effectiveness (MOE) of total 
cost (TC) and average customer wait time (ACWT). The 
Simulation was co-developed by the author and LT Glenn 
Robillard, and was designed to emulate the portions of the 
Navy's Uniform Inventory Control Program (UICP) applicable to 
this research. The simulation represents the demand process 
of a hypothetical Navy managed consumable item. The period of 
time over which demand is simulated and the characteristics of 
the item are specified by the user during the initialization 
of the simulation. Measures of effectiveness to be used in 
the performance comparison will be calculated from the actual 
cost and customer wait time data generated by the simulation. 
The UICP retention process and the various retention decision 
models will be tested ina variety of simulation scenarios. 
The scenarios are based on combinations of: 

- unit price 

- mean quarterly demand 

- variance of quarterly demand 

- patterns of declining mean quarterly demand 

- levels of excess inventory 

- inventory holding cost rate 

- obsolescence rate 

- administrative order cost rate 

- salvage rate 

The goal of this thesis is to determine how effective the 


Navy's retention logic is as compared with the four economic 


retention decision models. 


B. U. S. NAVY ECONOMIC RETENTION POLICY 

meorsevescea In Ehe introduction to this chapter, the 
Navy's Economic Retention policy has been ina state of flux 
for approximately nine years. The current RL for "essential" 
eens (1.e., item Mission @ssentiality Codes (IMEC) 3, 4, 
and 5) is set at eight years worth of annual forecasted 
demand, with ERS constrained to a minimum retention quantity 
of five units. All material that has been stocked in the 
Supply system for less than seven years 1S not subject to a 
retention limit. This material is retained until the seven 
year walting period has passed before being subject to 
retention review. 

Retention and disposal requirements are reviewed by the 
ICP semi-annually in conjunction with the execution of the 
March and September inventory Stratification, UICP application 
B20 . Stratification is the process of matching current 
inventory to requirements and categorizing inventory based on 
the type of requirement. DOD Regulation 4140.1-R [Ref. 7:p. 
4.3) defines the Stratification categories as Authorized 
Acquisition Objective (AAO), Economic Retention Stock (ERS), 
Contingency Retention Stock (CRS), and Potential Reutilization 
mmock (PRS). The Authorized Acquisition Objective is a 
combination of the peace-time requirements for U.S. Forces 
through the end of the second fiscal year following the 
Current date and the approved stockage requirements for grant- 


aid and military assistance programs. Economic Retention 


Stock is inventory held beyond the Authorized Acquisition 
Objective which is determined to be more economical to hold 
for future requirements as opposed to disposing and 
reprocuring in the future. Contingency Retention Stock is 
inventory held for known or potential requirements not covered 
by Authorized Acquisition Objective, such as initial 
outfitting, mobilization and Foreign Military Sales (FMS). 
Potential Reutilization Stock (also known as Potential Excess 
(PE)) is all inventory beyond the sum of the Authorized 
Acquisition Objective, Economic Retention Stock and 
Contingency Retention Stock. 

The ICPs will make the final retention/disposal decisions 
on material categorized as Potential Reutilization Stock. 
When a disposal release order 1S issued by the ICP, the depot 
holding the Potential Reutilization Stock will transfer the 
material to Defense Reutilization Marketing Office (DRMO) for 
Salvage or reuse. For this research all Potential 
Reutilization Stock is assumed to be sent immediately to DRMO 
move olikeloro sell - 

The calculation of Economic Retention Stock (ERS) 
performed during the UICP Stratification application is 


summarized aspfollows ikem co 


ERS = Max {(RL-D1-D2-D3-M), 5} 


Where: 


RL = eight years worth of forecasted annual demand. 
DI = forecasted demand, remainder of current year. 
D2 = annual forecasted demand, appropriation year. 
D3 = annual forecasted demand, budget year. 

M = reorder Objective, which equals the sum of 


safety stock, leadtime demand, and an economic 
order guantity (EOQ). 


Mmmiewecalculation for Economic Retention Stock (Equation 
imers based On recurring demand and does not take into 
account the portions of the Authorized Acquisition Objective 
which are considered non-recurring demand, such as Preplanned 
Program Requirements (PPR), Prepositioned War Reserves (PWR), 
Other War Reserves (OWR) and outstanding backorders (Due-out). 
In addition, Equation 1.1 constrains the Economic Retention 
Stock to a minimum of five units, to ensure a minimal buffer 
or safety stock iS maintained for "essential" material. The 
actual amount of inventory held is equal to the sum of 
Authorized Acguisition Objective, Economic Retention Stock and 
Contingency Retention Stock (where Authorized Acquisition 
Objective plus Economic Retention Stock equals the System 
Retention Level). By placing the five unit minimum constraint 
Seeeeconomic Retention Stock, the System Retention Level is 
also constrained to a minimum of five units. For this thesis 
Planned Program Reguirements, Prepositioned War Reserves, 
Other War Reserves and Contingency Retention Stock were 
assumed to be zero. 

Because the key to the amount of inventory categorized as 


memomac Retention Stock and Potential Reutilization Stock is 


the RL, this research will] focus on alternative methods: 
calculating a RL through the use of Economic Retention 
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C. ORGANIZATION OF RESEARCH 

The remainder of this thesis will be devoted to the 
discussion of mathematical economic retention models, the 
development of the analytical approach and simulation, and the 
presentation of the simulation results and concluSions. 
Chapter II reviews various mathematical models and discusses 
selection of the models chosen for the research. Chapter III 
develops the analytical approach to be used in comparing the 
UICP retention process to the mathematical models chosen in 
Chapwen. Ii. Chapter IV provides a description of the 
Simulation, to include a discussion of the major procedures 
and algorithms used. Chapters V and VI present the simulation 
results. Finally, conclusions and recommendations are 


presented in Chapter VII. 


A i oe ECONOMIC RETENTION DECISION MODELS 


A. LITERATURE REVIEW 

Excess inventories are an administrative and economic 
burden which consume valuable warehouse space, deplete working 
Capital and help to reduce inventory accuracy. In general, 
there are two causes for excess inventory. First, the demand 
rate may be overestimated due to a forecasting error, a change 
in technology or a change in operating tempo. Second, the 
Navy may obtain more units than they intend in a given 
replenishment action. This can happen as a result of errors 
in procurement document quantities or because the supplier 
delivers more units then the Navy requested. 

Mathematical models designed to represent the excess 
inventory problem are known as Economic Retention Decision 
Models. The objective of an Economic Retention Decision Model 
1s to reduce the administrative and economic burden of 
Carrying excess inventory through disposal of surplus stock. 
The approach to determining how much excess inventory to carry 
and how much should be disposed of varies from model to model. 
The basic idea behind most Economic Retention Decision Models 
is to determine the trade-off between the cost to dispose of 
material and the cost to hold material. What differs between 


models is how to define the cost to dispose of material and 


the cost to hold mawerioie While considerable literature 
exists on determining inventory retention levels, few 
researchers have directly addressed the Navy's’ excess 
inventory problem. 
1. Heyvaert and Hurt 

Heyvaert and Hurt developed one of the first models 
that treated the situation in which mean demand is declining, 
which is one of the causes of excess inventory [Ref. 9]. The 
model was designed to provide a simple, fast and accurate 
method for determining optimal stocking levels for slow-moving 
items. A unique objective function based on material storage 
costs and the cost of non-satisfaction of a demand was 
derived, with the optimal inventory levels (available level) 


being determined by minimizing the total cost function (W): 


Zao 
W = alt+ BPP 
; P ie 
% = >, (s-d/2)p, + (s*/2d) Dy 
=0 =8+1 
B = (d-s) pg 2 


=98+1 


Where: 
a = long run mean stock level, assuming variations in 
demand are linear. 
I = total cCeOSst Cov erere one Unie inges 


replenishment period (t). 
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replenishment period (t). 
-weOedl Gost resulting from non-Satisfaction of a 
demand requirement. 
S = current inventory on hand and on order (available 
level). 
= demand during a replenishment period (t). 
Py = Probability that an issue of size d will have to 
be made, assumes d has a poisson distribution 
Peeieimiedn = iN, 0.1 < it < 10.0. 


B = expected number of shortages during a 
F 


Although this model does not treat the problem of 
excess stock generated from reduced demand rate, the concept 
of determining optimality based on cost and customer 
Satisfaction helped motivate the use of total cost and ACWT as 
the MOEFs to be used in the performance comparison phase of 
this research. 

2. Rothkopf and Fromovitz 

The Rothkop£f and Fromovitz model for a save-discard 
@eersiOn involves a bulk commodity that comes in a rented 
Memeainer [Ref. 10]. Although this model is too specific to 
adapt to the Navy problem, it is one of the few models which 
Geals with the stochastic nature of demand. It also applies 
Mees concept of discounting future costs. 

S. Hart 

Hart designed a procedure to calculate a procurement 
schedule and retention quantity for a selected inventory item 
oeert. 11). The procedure minimizes the sum of discounted 
relevant costs which vary in amount or in timing with changes 


in the retention quantity. Relevant costs include the cost of 


ipa 


holding the retaimed GQUWantity, cost’ of "Monn scrappitg Heme 
retained quantity, cost of delaying the write-oft Gieeeee 
retained quantity (write-off occurs when the material is 
elther sold or scrapped), cost of procured quantities, and 
cost of holding the procured quantities. The minimum cost 
retention quantity is determined using a sequential search 
procedure based on the "Golden Section" method. For each 
retention quantity considered, a procurement schedule is 
determined heuristically according to a set of rules based on 
Economic Order Quantities and Economic "Bridging" Quantities. 
While Hart's model provides an interesting approach to the 
excess inventory problem, the level of effort required to 
incorporate his model into the Navy's UICP levels software 
application was beyond the scope of this research. 
4. Simpson 

Simpson's "formula" is one of the most frequently 
cited works in recent literature dealing with the excess 
inventory problem [Ref. 12]. The formula provides a clear and 
easy-to-use procedure which was originally developed for 
possible implementation by the Navy. 

The formula compares the cost of storing material, 
considering the chance that it may become obsolete and the 
cost of repurchasing the material in the future when needed, 
1f present surpluses are sold by disposal action today. An 


economic retention period formula was derived which equals the 


eZ 


cost (per dollar value of material) of retaining X years of 
stock (C,) less the cost (per dollar value of material) of 
disposing of X years of stock (C;). In the derivation of the 
formula it was assumed that future demand was known and 
constant, all general price levels and rates were also 


—~weant. The derivation is a follows: 


epee (1—p) *+z( (1-p) (1+i)*+(1-p)?(1+i)*71+....+(1-p)*(1+2)) 


Zz 


SeeGOst OF LelLaining X years of stock. 
er -secost Of disposing Of X years of stock. 


D fraction of present unit price of material which 
will be realized in disposal sales (i.e. 15 cents 
em the dollar, D = .15). 

p = fraction of material which will become obsolete in 
any one year. 

r = annual storage cost rate per dollar of material. 

1 = annual interest rate. 

X = Retention Level (RL). 


Equation 2.4 (C,.) represents the obsolescence cost and 
storage cost incurred from holding material for X years. The 
obsolescence cost term (1-(1-p)*) calculates the dollar value 
of loss due to obsolescence (per dollar of material) 
compounded over X years. The storage cost represents the 


cumulative cost of holding inventory X years, where the dollar 


Is 


value of inventory is reduced by p each year due to 
obsolescence, and includes the cost (compounded annually) of 
lost interest revenue from money used for storage costs. 
Equation 2.5 (C,) represents the cost (per dollar of 
material) of furnishing a given quantity of an item at time t, 
given material was disposed of at time t,. The C@grae. 
disposal is reduced by the return from disposal sales, which 
1S increased in value at the compound interest rate until t,. 
The value for X, the peel number of years stock to 
be retained (RL) is obtained by equating C, to Cy, and solving 


for X. Simpson gives the following such solution: 


Log ieee = | 2.6 
xX = i+p+r(1-p) (1+1) 


Log =F 
deaeel 


5. Mohon and Garg 





The Mohon and Garg model expanded on Simpson's 
economic retention period formula by considering the case in 
which shelf Jife’ 1s prebaomiterie [Rere f.1- They also 
derived the specific case in which shelf life is exponentially 


distributed. While the Mohon and Garg model may offer some 


‘Mohan and Garg assume shelf life is a function of 
obsolescence and deterioration. The Navy uses a combination 
of shelf life codes to account for deterioration of material 
and an obsolescence factor included in the system (UICP) 
holding cost rate. 
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improvements over Simpson's basic formula, it would be 
Gifficult to apply their model in the Navy's. UICP. 
Determining the appropriate probability distributions for 
obsolescence and deterioration rates to use with the expanded 
model would be a complex task. Because of this, a retention 
model which has robust performance with respect to 
obsolescence rate might be more appropriate for the Navy. 
6. Tersine and Toelle 

Tersine and Toelle developed two "net benefit" models 
Molt tering complexity for determining inventory retention 
levels [Ref. 14]. The models indicate how much inventory 
Should be held (economic time supply or RL) and how much 
Should be disposed of at a specific salvage price for a given 
item. In the derivation of both "net benefit" models it was 
assumed that future demand was known and constant, all general 
price levels and rates were alsSo constant, and no stockouts 
were permitted. 

The first or simple net benefit (NB) model calculates 
the economic time supply of material to hold that maximizes 
net benefit (cost savings) resulting from the sale of excess 
Stock. The formulation of the NB equation and the economic 


time supply (t,) 1s as follows: 


iss) 


Net Benefit = Salvage Revenue + Holding Cost Savings Jae 
-Repurchase Cost - Reorder Cost 


Salvage Revenue =G@Px= P2(M=tR) meer Ma eere 2.8 
2 -~q) ? 2.9 
Holding Cost Savings = MUPF _(M-q)" PF 
2R 2R 
_ M*PF_ RPFtC2meMOPF NOPFE 
2R 2 2R 2 
Repurchase Cost = Pg = PM-PRt 2.10 
22 
Reorder Cost sea 8 CRE 
Q Q OQ 
Where: 
q =M -tR = amount of excess inventory that is 
disposed of, in units. 
t = time supply, in years worth of inventory 
retained. 
tp = economic time supply in years worth of inventory 
retained (RL). 
C = ordering cost per order. 
F = annie ene ding cestartroachlen. 
M = available stock in units. 
P = unit cost of the item. 
P. = unit salvage value of the item. 
Q = economic order size in units. 
R = annual demand in units. 
The resulting net benefit formulation is as follows: 
_ _ RPFt?* ORFs CR\,., M*PF_MOPF CM 
f(t) = -———+PR-P,R+ = —_——_— —- —=—— +P M- PM- —— 
eo) 2 { on aes “2R 2R°* QO 
Zee, 


HKG 





Note that f(t) describes a parabola and therefore has a single 
maximum. By taking the first derivative of f(t) with respect 
to t and setting it equal to zero, the economic time supply 


(t,.) equals: 


Since the second derivative of f(t) 1s negative, t, is located 
at the maximum point. 

The second model, a present value net benefit (NB-NPV) 
model, compensates for the fact that investments occur at 
different points in time by discounting them to their present 
value. Under continuous compounding, the present value of a 
future purchase of an item with a current price (P) at time t 
is Pe'*)*, where i is the annual inflation rate and k is the 
discount rate. For this thesis inflation was assumed to be 
zero and the discount rate was set to seven percent. 

The formulation of the objective function of the net 


present value version of the net benefit model is as follows: 


_ PFtR(e-*t-1) PFQ PO+C (i-k) €_ 
Be) = ee lars gear |e PRO 


eee PFO MPOrC | ti-nwe 
2k 2(i-k) e(li-k)0/R_4 


ie 


Although Equation 2.14 cannot be solved directly “femuam 
Newton's method can be used iteratively to obtain a solution. 


Where: 


a Op es fia) fons 
f(t.) 

For this thesis the t, obtained from the NB model was 
divided by two and then used as an initial estimate for the 
NB-NPV model t,. The NB model t, was divided by two to ensure 
that the initial approximation to the NB-NPV model t, was 
sufficiently close to the optimal solution so that Newton's 
method would converge upon a solution. This choice of initial 
Starting solution was particularly important for the demand 
scenarios with low unit price, because the RLs for the NB-NPV 
model were expected to be significantly less than the 
respective RLs for the NB model. Successive values for t were 
calculated umtil [t,4-t,| < 90.01. When this stopping 
condition was satisfied, the final t, for the NB-NPV model was 
set equal to t,,,. 

Although the Navy UICP assumes that demand is 
stochastic and allows for stockouts, Tersine and Toelle's "net 
benefit" models are well suited for application in the Navy's 
UICP. In an effort to account for the potential for stoGgkouE 
due to the stochastic nature of demand typically associated 
with a Navy managed item, a modified "net benefit" (NB-MOD) 


model was developed. 
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Disposal of Some quantity of excess inventory will 
Cause the inventory position (IP) to reach the reorder point 
(RO) prior to the time it would have reached the RO without 
the disposal of the excess inventory. Therefore, with 
disposal the inventory system will experience one or more 
additional reorder cycles, depending on the quantity disposed. 
Because of the stochastic nature of demand, every additional 
reorder cycle exposes the inventory system to an increase in 
the number of possible stockouts. In the modification of the 
NB model, for every additional reorder cycle that occurs due 
to disposal, the net benefit from disposal is reduced by the 
expected additional shortage costs. Themodified formulation 


(NB-MOD) is: 


Net Benefit(MOD) = Salvage Revenue + Holding Cost Savings 
- Repurchase Cost - Reorder Cost 
- Shortage Cost 


Urls 

The new term, shortage cost, is a linear function of 

the number of additional reorders (N) that are made due to the 
disposal of gq units worth of stock. We must first calculate 


iN: 





M_ (M-q) 
R R M-tkR 
N= = e ihn 
ie 0 
R 


ig 


N = number of additional reorders required due 
to the original disposal Of g units: 


M/R = mean time supply of material without 
disposal. 

(M-q)/R = mean time supply of material with disposal. 

C7 = mean time between reorders. 

E[X>RO] = expected number of shortages in a reorder 
cycle. 

RO = reorder polne 

A - shortage cost per unit. 

x = actual demand during a procurement 


leadtime. 


Now we may obtain the shortage cost: 


Shortage Cost = NA(E[X>RO] ) Dae 


The expected number of shortages (E[x>RO]) in a 
reorder cycle, assuming that X is normally distributed with 


mean, W and variance, o* is given by [Ref. 15]: 


E[x>RO] = (p-RO) xP (2> <) +oxe Zame] 


Where: 


P (2> ie = Probability of a stockout. 


20 


f | = Standard normal distribution function 


evaluated at wl 


RO] "Rb + OZ. 

fee Staencard normal distribution value which 
Sateoiies Ene UNCP “probability of a stockout"’* 
expression for a given values of R, L, u, ov, F, 
P, A, and E. 

Seiicen leadtime demand’. 

Voiaiitee me emlondtaime demand . 

= procurement leadtime demand in years. 


car 
II 


Because the term E[x>RO] in Equation 2.20 is not a function of 
t, the expected number of shortages in a reorder cycle is 
Created as a constant. 

Collecting these terms together, the objective 


meim@erionl Of the modified net benefit model is: 














2 2 
f(t) = -e +(PR-P,Rs ae Se 
2.20 
CM _|M-tR 
pm-pm-<i A(E[x> RO 
| a PRO RO) 


The UICP levels application calculates the probability 
of stockout using the following expression: FP/(FP+AE), where 
Mememene annual holding cost fraction, P is the unit cost of 
an item, A is the shortage cost per unit and E is the military 
essentiality. 


“In UICP this parameter is PPV. 


°In UICP this parameter is BO19A. 


ou 


Next we must determine if Equation 2.20 is a parabola. 
Note that Equation 2.20 can be expressed in the form at-+bt+c 
and thus 18 a parabola [Ref. 16,p.39]. By grouping terms 


appropriately we obtain the constants a, b, and c: 








2 = -~RPF) Jia 
2 
b = PR-p,R+ PF, CR Ra (pix Ro} ) 2 ee 
2 QO @Q 
M*PF MOQPF pp ye) 
SR aR *PeM- PM--S- aE Lx? RO) ) 


By taking the first derivative of f(t) (Equation 2.20) 
with respect to t, setting 1t G@qual to zero and solvingmtouges 


the modified economic time supply (t,) 1s obtained: 


P-P,, Q , C+A(E[X>RO] ) yi 
PF 2R QPF 


Since the second derivative of f(t) 1S negative, t, 1s located 
at the maximum: Deane 
7. Silver and Peterson 
Silver and Peterson developed a rule for the disposal 
of excess inventory which, while derived using a different 
approach from that of Tersine and Toelle, yields the same 


numerical results [Ref. 17:Chap. 9]. In a manner similar to 


GOD 


Simpson's approach, Silver and Peterson focused on the cost of 
no disposal (C,,) versus the cost of disposal (C,). Then, 


assuming an EOO strategy with deterministic demand, Silver and 





Peterson formulated an objective function of Cy, - C, , where: 
- peasy 9 2.25 

ep 2D 
DIS 





C, = -gw+( a (=—*) vr+e ( /ZADVE+Dv) 


= 
= 
OD 
ry 
OD 


@oct Of no disposed . 

COsteot ciuspoosa! . 

amount of excess inventory to dispose in units. 
on hand inventory in units. 

expected annual demand in units. 

unit price. 

Salvage value per unit. 

MoOlCnWiGmeOct erate S/S/YV . 

administrative order cost per order. 
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The last term in C, represents the inventory holding cost, the 
Peumistrative ordering cost and the repurchase cost of the 
Smeem Gisposed (W) incurred after the stock retained is 
exhausted (which occurs at time (I-W)/D and continues until 
time I/D). The inventory holding cost and the administrative 
ordering cost are calculated assuming an EOO strategy. The 
repurchase cost of the stock disposed (W) is. calculated 
assuming the repurchase unit cost equals the unit cost at the 
meme, OL disposal. 


By taking the first derivative of the objective 


ZS 


function (Cy, - C)) with respect to W and Setting 12t equemeeee 
zero we obtain Silver and Peterson's "decision rule for 


disposal," an expression for W, which maximizes C,, - C,. 
w = r-Eog-L\V-9) 2G 
VI 


Although Silver and Peterson used a different approach 
in the formulation of their model than Tersine and Toelle, it 
can be show that Silver and Peterson's "decision rule for 
disposal" and Tersine and Toelle's simple "net benefit" model 
yield the same results. Using Silver and Peterson's notation 
1t can be shown that Tersine and Toelle's economic time supply 
(to) multiplied by annual demand (D) equals Silver and 
Peterson's equation for the amount of inventory to retain (I- 
W), as follows: 


t)xD = D(v-g) , DA_, EOQ 
VI VIEOQ 2 


substituting {= for EOQ yeilds 


t.xp = 2 


) 4 
E =z * lave 


= Divied) +> 
VI 2VI 

= Piv'9) 4 Fog = I-w 
VI 


QED 
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Because the two derivations result in the same 
economic retention decision, only the notation from one 
derivation was used in the thesis. Tersine and Toelle's 
notation and approach was chosen, primarily because of the 
extensive background provided on the excess inventory problem 
and the thorough development of the derivation of their model. 

8. Rosenfield 

Rosenfield developed a model for the optimal number of 
jmeems to retain for slow moving or obsolete inventories under 
conditions of stochastic demand and perishability (shelf-life) 
[Ref. 18]. This model is one of the few that addresses the 
Drobabilistic nature of demand for the general excess 
inventory problem. Rosenfield's basic model assumes that 
episodes of demand can be represented by a renewal process. 
This allows for a variable number of units demanded per 
episode. The model determines the correct number of units to 
retain. In the model a unit 1S worth disposing of if its 
immediate salvage value (it's present resale value) exceeds 
it's expected discounted sales value (from a future sale if 
the unit is held in inventory) minus the expected holding 
costs to be incurred (until the time of sale). 

Because Rosenfield's final expression for the number 
of units to retain contains the moment generating function for 
the distribution of time between demand episodes, the model 


becomes complex when the distribution of demand episodes is 


Ze 


not a Poisson distribueven Although this model may have 
application to the Navy's excess inventory problem, the level 
of effort required to incorporate Rosenfield's model into the 
Navy's UICP levels software application was beyond the scope 


of Ehis reseemche 


B. SUMMARY 
The mathematical models chosen for this research were 
based on their applicability to the Navy's excess inventory 
problem, the UICP model, and the simulation. The models 
chosen were: 
- Simpson's “economic retention period formula" (TRAD). 
- Tersine and Toelle's simple "net benefit" model (NB) 


- Tersine and Toelle's present value "net benefit" model 
(NB-NPV). 


- The modified "net benefit" (NB-MOD), a version of the 
Simple "net benefit" model. 


These models, together with the Navy's UICP current retention 
logic, will be referred to as the "models" throughout the 
remainder of the thesis. 

Although the UICP model was developed under the assumption 
that demand is stochastic, all the mathematical models listed 
above were developed under the assumption that demand was 
deterministic (with the exception of NB-MOD). The decision to 
use primarily deterministic models was based on two factors. 
First, as Simpson [Ref. 12] discussed, the effect the 


deterministic assumption has on a Retention Level (RL) 1s not 
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Sumer cant. Secondly, the difficulty of incorporating into 
the UICP model and into the simulation the stochastic models 
reviewed does not justify the small improvement in accuracy 
which, according to Simpson, we would experience. Because a 
true stochastic economic retention model was not used in this 
research, a Total Cost Analysis (see Chapter III.C.1) was 
Gonducted to develop a baseline, with respect to cost, to 
evaluate how the deterministic models actually perform in a 


Beocmestic environment. 
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| ae Be RESEARCH APPROACH AND ANALYTICAL METHOD 


A. OVERVIEW 

The analysis that was done for this thesis made use of a 
Simulation that was written in Pascal. The simulation was 
developed to represent the Navy's UICP model as well as the 
mathematical models that were analyzed in this research. A 
complete discussion of the simulation program is contained in 
Chapter Iv: 

The analysis and performance comparisons of the models 
were based on MOEsS calculated from simulated data for six 
basic demand scenarios. For each demand scenario four 
retention scenarios were analyzed using the simulation. A 
Total Cost Analysis was performed to determine the optimal 
amount of inventory (from just the cost standpoint) to hold 
for a given quantity of initial excess inventory. A Constant 
Demand Analysis was performed to compare the various models to 
the theoretically optimal retention level that was determined 
during the Total Cost Analysis. The same input parameter 
values were used in the Constant Demand Analysis as in the 
Total Cost Analysis. A Declining Demand Analysis was 
performed to compare the models in three scenarios (patterns) 
of declining mean demand. Finally, Sensitivity Analysis was 


performed on various combinations of demand scenario, pattern 
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of declining mean demand, and the parameters of administrative 
reorder cost rate, salvage rate, inventory holding cost rate, 
and obsolescence rate. (A complete discussion of the 
Sensitivity Analysis 1S contained in Chapter VI.) 

Table 1 provides a summary of retention scenarios, cross 
rererenced by demand scenario and mean quarterly demand 
pattern. Pileheciimay) Imeetme table represents a set of 
Simulations and will be referred to as a simulation setting. 
The meanings of the demand scenario acronyms can be found in 
Table 2. A summary of the 16 specific settings to be 
Gensidered in the Sensitivity Analysis 1S provided in Chapter 
i Table 9. 

In the performance comparison phase of the research the 
models were ranked based on the MOEs of total cost and ACWT. 
The comparisons were done by demand scenario for the results 
from the analysis scenarios of Constant Demand Analysis, 
Declining Demand Analysis, and Sensitivity Analysis. Multi- 
Attribute Decision Making techniques and hypothesis tests 
Mesea On a paired difference t-test were used to compare the 


performance of the models. 


B. DEMAND SCENARIOS 

Items managed by the Navy are assigned a Navy Mark Code 
based on unit price and mean quarterly demand. The Mark Code 
indicates the probability distribution for leadtime demand and 


the inventory level setting method to be used in the UICP 
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model [Ref. 19:p. 3-9]. Six hypothetical) items basea cia 
Mark Code deSignation criteria were selected for use 
throughout the research. The hypothetical items, called 
demand scenarios, were chosen so that the effect of varying 
level setting computation methods, unit price and mean 
quarterly demand on economic retention decisions could be 
analyzed. The demand scenarios described in Table 2 are a 


function of the probability distribution of demand episodes, 


TABLE 1. SUMMARY OF SIMULATION SETTINGS 


DEMAND 
SCENARIO 

HDAVEP | Sor ie Oi Sei riot ie 
PATTERN 

| TCA TCA TCA 


CONSTANT 
MEAN 
DEMAND CDA CDA CDA 


| DECLINING 
|} MEAN 


DECLINING DDA DDA 
MEAN 
DEMAND SA (16) SA(16) 
LGONY Ea 


| DECLINING || DDA | | Ve DDA 
MEAN 
DEMAND SA(16) SA(16) 
"CONCAVE" 


TCA = Total Cost Analysis, CDA = Constant Demand 
Analysis, DDA = Declining Demand Analysis, SA = 
Sensitivity Analysis (16 simulation settings for 
each demand scenario and demand pattern combi- 
Hatliom:. 
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mean quarterly demand (high and low), variance of quarterly 
demand (high and low), and unit price (high and low). Demand 
Variance for the demand scenarios with a normal distribution 
are classified as high (with a standard deviation to mean 
ratio of 1.25) and low (with a standard deviation to mean 


memo of 0.30) [Ref. 20). 
TABLE 2. DEMAND SCENARIOS 


NAVY | PROBABILITY DEMAND UNIT ACRONYM 
MARK |DISTRIBUTION]|QUARTERLY| VARIANCE | PRICE(S) 


High: 20| High: 625|High:1500| HDHVHP_ 
High: 20| tow: 36|H#igh:1500| HDLVHP_ 


Normal High: 20 | Low: 36 | Low: HDLVLP 


2 
p3 | roisson | row: 2| wa |wign:1500 


High: 20 | High: 625 





C. ANALYSIS SCENARIOS 
1. Total Cost Analysis 

This analysis was performed to compute a total cost 
for 100 quarters of demand activity for a given demand 
scenario based on the following set of assumptions. Assume at 
time zero the inventory 1S i1n an excess position and an 
immediate retention/disposal decision is made. Next, assume 
that this is followed by 100 quarters of demand activity with 


a stationary quarterly mean demand. The initial on-hand 


al 


inventory selected for demand scenarios with high unit price 
was equal to 20 years of average annual demand. For demand 
scenarios with low unit price, the initial inventory was equal 
to 25 years of average annual demand. A total cost was 
calculated for various retention levels beginning with a level 
equal to 0.5 years of annual demand and continuing, in 
increasing increments of 0.5 years annual demand. Retention 
levels were not increased beyond the inventory on hand at time 
Zero. Based on an initial inventory of 20 years worth of 
annual demand for the demand scenarios with high unit price, 
40 total cost! data points (retention levels) were calculated. 
These data points were used to construct total cost curves for 
the demand scenarios with high unit price. Based on an 
initial inventory of 25 years worth of annual demand for the 
demand scenarios with low unit price, 50 total cost data 
points (retention levels) were calculated. These data points 
were used to construct total costs curves for the demand 
scenarios with low unit price. 

Each total cost data point 1s discounted to Guzers 
year dollars and is equal to the sum of material cost, 
administrative ordering cost, inventory holding cost, shortage 
cost and salvage revenue which accrue over a simulation period 


(See Equations 3.1 and 3.2). The total cost data points for 


'The total cost figure used for each data point is the 
average total cost over all replications of the respective 
Simulation. 
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each demand scenario were then plotted to form a total cost 
curve (See Appendix E, Graphs 13 through 24). The goal of the 
Total Cost Analysis was to determine if a minimum total cost 
associated with a single retention level existed in a 
stochastic demand environment in the same way as shown by 
Tersine for the deterministic case [Ref. 14]. The minimum of 
each total cost curve was used to obtain the optimal retention 
level for each demand scenario. These optimal retention 
levels were used as a benchmark for comparing the performance 
of the models in the Constant Demand Analysis phase. 
2. Constant Demand Analysis 

This analysis was designed to compare the performance 
of the models to the performance of the optimal retention 
level determined in the Total Cost Analysis. The comparison 
was done for all combinations of the demand scenarios and the 
models under the same simulation settings that were used in 
the Total Cost Analysis. The goal of this analysis was to 
determine, for each demand scenario, how the models performed 
in the Navy's stochastic demand environment with respect to 
the optimal retention level. 

3. Declining Demand Analysis 

This analysis was designed to compare the models under 
a scenario involving declining mean quarterly demand. Three 
patterns of declining demand where developed for this 


analysis. The declining demand patterns represent possible 
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effects the reduction in Naval Forces and budget might have on 
demand for Navy managed items. In Appendix E, Graphs 1 
through 6 depict the six patterns of declining demand that 
were used. Demand activity for these scenarios begins with a 
pattern of 30 quarters of stationary mean quarterly demand. 
This allows the simulation model to reach steady state as 
discussed in Chapter IV. This was followed by 20 quarters 
with declining mean quarterly demand and finished with 16 
quarters of constamt mean quarterly demand. The J/6S"@ieneme 
period was included to allow the determination of the long 
term effect that a specific retention policy might have on 
performance. Over the period of the decline of the mean 
quarterly demand, for demand scenarios with a high mean 
demand, the demand decreased from a mean of 20 units per 
quarter to a mean of 2 units per quarter. The mean quarterly 
demand for demand scenarios with low demand decreased from a 
mean of 2.0 unitsS per quarter to a mean of 0.2 units per 
Quarter. The comparison of model performance was done for 
all combinations of the demand scenarios, models, and decline 


patterns. 


D. PERFORMANCE COMPARISONS 

The concept behind the performance comparisons is to 
provide Navy inventory modelers with some quantitative data 
that will help them select the most suitable model to use in 


a given situation. The use of total cost and ACWT as the MOEs 
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was motivated by two factors. The first was Heyvaert and 
Hart's use of cost and customer satisfaction in the 
development of their model [Ref. 9], which in essence asserts 
that when evaluating a model total cost is not the only 
evaluation criteria to consider. Modelers should also 
consider how a model satisfies customer requirements. The 
second was the fact that total cost and ACWT are generally of 
primary concern to the managers at the Navy's inventory 
control points when they make inventory policy decisions. 
The total cost MOE (Equations 3.1 & 3.2) is based on the 
Navy's UICP model total cost objective function [Ref. 19:p. 3- 
A-4]. Total cost is discounted to current year dollars and is 
equal to the sum of material cost, administrative ordering 
COSt, inventory holding cost, shortage cost and salvage 
revenue which accrue over a simulation period. Costs were 
discounted because of the length of time (simulation period) 
over which the analysis was performed. Additionally, costs 
were discounted to evaluate the effect, over time, the models' 


Varying disposal decisions had on total cost. 





g — ( EHP sc 
TC(D) = (Q,P+C,At ( J |en$2-p PF 


S15) 


Where: 


TC(D) = total discounted cost for one®repiaicatienees 
a Simulation given D units disposed during 
the simulation period. 

GES COUnTA LACES . 

number of units ordered during quarter k. 
bliaulie worengee< 

administrative order cost. 

number of orders placed during quarter k. 
inventory on hand at the end of week j. 
holding Gost fraction (S/untr—yiaie 

time Weighted Units Short (TWUS) for quarter 
k, see Equation 3.4. 

shortage cost ($/unit-yr). 

number of units disposed of during quarter k. 
Salvage rate (a fraction of P). 

discount rate. 

number of quarters simulated. 

summation index for 13 weeks of a quarter. 
Summation index for the number of quarters 
Simulated. 


~~ 
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The ACWT measures the mean time required, in days, for the 
wholesale supply system to meet customer demands. ACWT for 
one replication of a Simulation equals the time weighted units 
short (TWUS) divided by the total demand (D) over the 
Simulation period (Equations 3.3 & 3.4). The simulation ACWT 


was equal to the average of all replication ACWTs. 





: Ba 
ACWT = a 
- 4 
TWUS = » [ (RD,-BOD,;) xAR,] 2 
1=0 
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n Se iiumoem@ et bag@korders (in units) for 
measurement period. 

RD, = receipt date of the i’ backorder. 

BOD, = date the it’ backorder occurred. 

Peeeeemamount Of i°° backorder (in units) filled on 


RD.. 


The actual performance comparisons were done using two 
methods. One method is the paired difference t-test and the 
other method is Multi-Attribute Decision Making (MADM). 

1. Paired Difference t-Test 

Hypothesis tests based on a paired difference t-test 
Statistic (Ref. 21:p. 572] were conducted on the results of 
the Constant Demand Analysis, Declining Demand Analysis, and 
Sensitivity Analysis simulations to determine which model(s) 
performed better than all others in each MOE category. Given 
that model "X" had the best result for a specific MOE, the 
null hypothesis was that the corresponding result, for every 
other model was equal. The alternative hypothesis was that 
the corresponding result, for every other model was not equal 
to the result for model "X." 

The paired difference t-test was used because there 
was dependence between the MOE results of the models for each 
setting simulated. The dependence was attributed to the fact 
that for each replication of a simulation, the randomly 
generated demand streams were identical for all the models 
within a setting. Further discussion of the relationship 


between random number generation and the dependency of results 
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1s contained in Chapter IV. 
2. Multi-Attribute Decision Making (MADM) 

In order to compare the models performance, the 
decision analysis technique known as Multi-Attribute Decision 
Making (MADM), a subset of the decision making processes known 
as Multi Criteria Decision Making (MCDM), was used. There are 
four characteristics which make this performance comparison a 
Multi-Criteria Decision Making problem [Ref. 22,p. 2]. First, 
there are multiple attributes (MOEs of total cost and average 
customer wait time). Second, there is conflict among the 
MOEs, 1.e. the higher the TC (which is bad) the lower the ACWT 
(which™@ie Gooai= Third, the MOEs have different units of 
measure (TC 1S per simulation period and ACWT is in terms of 
days per requisition). Fourth, the selection of the best 
model is to be made based on each model's level of achievement 
in the MOEFs of TC and ACWT [Ref. 22,p. 3]. The primary 
feature which makes the model selection decision a MADM 
process 1s that there are a limited number of predetermined 
alternatives [Ref. 22,p. 3]. In this case the alternatives 
are the retention models being analyzed. By using the MADM 
technique a final decision (model selection) can be made. 

The Simple Additive Weighting Method, one of the best 
known and widely used methods of MADM, was the method used for 
this thesis [Ref. 22,p. 99-103]. To determine a preferred 


model, a decision matrix must be constructed that incliugG@esaam 
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MOE values for each model. Because the Simple Additive 
Weighting Method requires a comparable scale for all elements 
in the decision matrix, a comparable scale matrix is obtained 
using Equation 3.7 to convert the MOE values to comparable 
Units. In addition to the comparable scale decision matrix, 
a set of importance weights are assigned to the MOES, w = 
Mieewayps- Lt Should be noted that w is normalized to sum to 
one. The weights should reflect the decision makers marginal 
worth assessment for each MOE. A total score (weighted 
average) for each model (A,;) and the most preferred model (A*) 


can be determined as follows: 


35 
A* = max {A,|Vi = 1,...,m} 
4 Size 
2, WF a3 
A; = — 
>», 3 229 
j=l 
r;,= min {x,,|Vi = 1, rm} / X;, 
Where 
m = the number of models being analyzed. 
i = the it? model of the m models. 
jmee=stme MORES of TC (J=1) and ACWT (J=2). 
w,; = the importance weight for the j‘? MOE. 
r,, = the comparable scale value for the j™ MOE of 


the i‘? model. 
foe= ene 7 MOE value for the i™ model. 


Although MOE results (x;;) are transformed onto a 


comparable scale (r,,;) by Equation 3.7, the decision makers 
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perspective regarding a difference of 0.2 between two model's 
r,, for the attribute of ACWT may not have the _ same 
Significance as a difference of 0.2 between the same model's 
r,, for the attribute of TC. For example, if the ACWT x,, is 
1.0 day in Model 1 and 0.8 days in Model 2 and the TC x,, 1s 
$80,000.00 in Model 1 and $100,000.00 in Model 2, a decision 
maker would probably consider the change in the TC x,,S to be 
more Significant. But if TC and ACWT are weighted equally 
Model 1 and Model 2 would have the same A;. The key to making 
effective use of MADM technigues is selecting proper MOE 
weights. Weights should be chosen to reflect the relative 
Significance of trade-offs between TC and ACWT. 

Because the selection of MOE weights 1S somewhat 
Subjective and could vary between decision makers, three sets 
of weights were used when comparing the performance of the 
models (see Table 3). The use of three sets of weights will 
Show the sensitivity of model selection to MOE weights. The 
sensitivity of model selection to changes in MOE weighting 
Should also identify models which perform better with respect 


EO total cose or 2e a 


TABLE 3. MADM MOE WEIGHT SETS 
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Due to the subjective nature of MOE weight selection 
and the difficulty of determining the relative significance of 
trade-offs between ACWT and TC between various models, the 
MADM results should not be considered a solution to the 
problem. For this thesis the results were used to help 
develop criteria for selecting a model based on demand 
scenario and the decision maker's emphasis on the MOEs of TC 


and ACWT. 
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Ve. SIMULATION 


A. SIMULATION STRUCTURE 

A discrete event Monte-Carlo simulation was used to obtain 
Statistical estimates of the values of the measures of 
effectiveness used in the thesis. The events of the 
Simulation occurred on a quarterly basis and were defined by 
the activities associated with the UICP demand process. 

The main routine of the simulation was representative of 
the actions which occur in the Navy's UICP model given the 
quarterly generated demand observations. Execution of these 
actions is controlled by two "for" loops. The OUG=ixiiaecds: 
loop controlled the number of replications of the simulation 
to be run. The inner "for" loop performed the functions of a 
Simulation clock and timing routine, where each increment of 
the inner "for" loop represented one quarter. The major 
procedures which are called in the timing routine are: Demand 
Observation Generation, Demand Forecasting, Inventory Level 
Setting (Levels), and Supply/Demand Review (SDR). A complete 
copy of the simulation is included in Appendix D (The Pascal 
code can be obtained from Navy Ships Parts Control Center, 


Code 046, Mechanicsburg, PA 17055-0788). 
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1. Demand Observation Generation 

Demand observations for the number of quarters 
Simulated, for each replication of a simulation, are generated 
using an appropriately transformed pseudo-random number 
generator. The resulting demand stream is a function of the 
Mee@eaoi lity distribution that is selected (Normal or Poisson), 
the mean quarterly demand, and the variance of demand. The 
probability distribution, mean quarterly demand, and variance 
of demand are specified during initialization of the 
Simulation. The method for generating a unique demand stream 
for each replication of a simulation is discussed later in 
this section. 

The algorithm for generating demand observations with 


EeeeOorsson(A) distribution was based on the relationship 
between the Poisson(A) and Exponential(1/A) distributions 


Seem 2s:p. 503): 


Meee 6”, b = 1, and i= 0. 
Generate U,,, ~ U(Q,1) and replace b by Dbu,,,. 
mio <a, reeurn X = 1. 
Otherwise, go to step 3. 
Be Replace 1 by 1+ 1 and go back to step 2. 
@ie algorithm returns X, when the Zj_,(-log(U,)) is less than 


A (equivalently, when ITj_, (U;) eee eee eee log (U,) Ss 


are exponential, they can be interpreted as the interarrival 


MmmeceOfL a Poisson process having rate 1. Therefore, X = X(A) 


1S a Poisson random variate equal to the number of events that 
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have occurred by time id. 


The algorithm for generating demand observations with 
a Normal distribution was based on the "polar method" 
[ReE.. 235 0eee oe 


1. Generate U, and U., as IID U(0,1), 
let Vos 2U, => 1 or Veranciae, 
and tet. =e an ae 
2. te WS], CC ode wemoren 
Otherwise, let Y = [(-2ln(W))/W]'”, 
X, = V\Y and X, = V.Y. 
Then X, and X, are IID N(0,1) random variates. 
The Uniform (U(0,1)) random number generator used in 


the Poisson and Normal random variate algorithms 1S a prime 
modulus multiplicative linear congruential generator Z[1i] = 
(630360016 * Z[i-1]) (mod 2147483647), based on Marse & 
Robert's portable FORTRAN random number generator UNIRAN [Ref. 
2 eee The simulation has the capability to produce 
20,000 unique seeds for the random number generator based on 
the NXSEED function, also from Marse & Roberts [Ref. 23:p. 
456]. Using the NXSEED function, a unique demand streams for 
each replication of a simulation 1s generated by reseeding the 
random number generator with a new seed prior to generating 
the next replication demand stream. A further discussion of 
seed selection and unique demand stream generation 1s 
contained in Section ieee 

Because the internal execution of the Supply/Demand 


Review procedure is on a weekly basis, each quarterly random 
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demand observation 1s subdivided into a 13 week demand stream 


oe 2@l lows: 


meeesOr 1 = 1 to 13, the demand observation for 
week(1) = 0. 

2. For i= 1 to current quarter's demand observation 
a. Generate a random uniform integer(X) from 1 to 


ie 
b. increment the demand observation for week (xX) 


by one. 


This routine randomly disperses one quarters worth of demand 
throughout the 13 weeks of a quarter. 

Am Option at simulation initialigatwon 1s to include 
one to five trend periods and/or one to five step changes in 
mean quarterly demand (D[t], where t equals a specific 
Quarter). The trend function follows an exponential growth 


pattern of the form [Ref. 24]: 


D[{t] = M * (1+A*t(0) 7) al 


Where: 

Mo = Initial Trend Mean, the mean quarterly demand 
at the beginning quarter of a trend period. 

A = trend coefficient. 

t(0) = at the beginning of each trend period this 
variable 1s reset to one and incremented by one 
at each quarter during a trend period. 

B =—eeema Dower function. 


The number of trend periods, the quarters in which a trend 


Starts and stops, and the parameters A and B for each trend 
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period are specified during initialization of the simulation. 
The step function applies a step multiplier (any non-negative 
number) to D[t-1] to determine D[t] [Ref. 24]. The number of 
steps, the quarter in which the step occurs (D[t]) and the 
step multiplier are specified during initialization of the 
Simulation. 
2. Forecasting and Inventory Levels Setting 

This part of the simulation was written to emulate, as 
closely as possible, the forecasting and cyclic levels 
application (D0O1) of the UICP model. 

a. Forecasting 

NAVSUP Publication 553 [Ref. 19:Chap. 3] G@mt@amaae 
general background information on the forecasting application 
in the D0O1 application. Single exponential smoothing or a 
moving average 1s used to forecast mean quarterly demand, 
depending on the results of step and trend tests. Single 
exponential smoothing or a power rule is used to forecast Mean 
absolute deviation of demand (MAD), depending on the results 
of step and trend tests. A smoothing constant of 0.01 was 
used for exponential smoothing in the simulation. 

Prior to actual computation of the next quarterly 
demand forecast, the most recent quarterly demand observation 
1S examined by two processes: "Step" filtering [Ref. 19:Chap. 
3]; and the Kendall trend detection test [Ref. 25]. These 


tests are used to determine if there has been a change in mean 
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quarterly demand that is significant enough to warrant 
discarding most of the historical demand data and to recompute 
the forecast using only recent data. When the process is "out 
Of filter" or a trend is detected a four quarter moving 
average 1S used to compute the next forecasted mean quarterly 
Meee The MAD is then forecasted using a power rule [Ref. 
270) - 
b. Levels Computation 

NAVSUP Publication 553 [Ref. 19:Chap. 3] contains 
a description of the Levels computation application in the 
DO1. The purpose of this part of the software 1S to compute, 
for a given Navy managed item, the economic order quantity and 
reorder point for the next quarter. The UICP calculations for 
inventory levels were developed within the guidelines of DOD 
Instruction 4140.39. Note that these guidelines follow an 
approach used by Hadley and Whitin [Ref. 27]. The optimal 
inventory levels are determined by minimizing an average 
annual variable cost equation composed of ordering, holding, 
and shortage costs. The level setting calculations in the 
Simulation are based on FMSO Level Setting Model Functional 
Description PD82 [Ref. 28] which was written by the Navy Fleet 
Material Support Office. Executable code obtained from the 
Navy Ships Parts Control Center (Code 046) was used in the 


Simulation to perform the actual level setting calculations. 
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Be Supply/Demand Review (SDR) 

The SDR routine of the simulation was coded to 
replicate the UICP model when processing material receipts, 
issues, and orders. In addition, a material disposal function 
was 1ncorporated in the routine. The disposal function occurs 
bi-annually in conjunction with inventory stratification and 
executes economic retention decisions. The events in the 
SDR routine are driven by the output from the Demand 
Observation Generation, Forecasting, and Levels routines for 
the respective quarter. The SDR routine 1s called once a week 
during each quarter and the events occur in the following 
sequence: material disposal (this disposal routine is used 
only during the first week of the first and third quarters of 
each year), receiving, iSsuing, and ordering. In addition, 
the SDR routine calculates and records data for TWUS, ACWT, 
aAnGdmMEOraleeost ; 

a. Material Disposals 

A semi-annual inventory stratification was 
performed to determine the "retention level" and to calculate 
the amount of "potential excess." The economic retention 
model specified during initialization of the simulation is 
used to perform these calculations. The models available in 
the simulation are: 

= sue © 
- Optimal 


- Traditional (TRAB) 
~ Net Benefit (NB) 
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- Net Benefit-Mod (NB-MOD) 
- Net Benefit-NPV (NB-NPV) 


For simulation purposes all "potential excess" is 
disposed of immediately and revenue from disposal is 
determined by multiplying the unit price of the item by the 
Quantity disposed and the salvage rate (salvage rate is 
Seeectmitead by the user during initialization of the 
Smelat1On). Total cost for the simulation period is reduced 
by the discounted revenue recognized from disposal. 

b. Material Receipt 

Outstanding reorders are maintained ina "priority 
heap" [Ref. 29:p. 149] in order of scheduled receipt date. If 
an outstanding reorder is due in the current week, the reorder 
1s removed from the outstanding reorder heap. The receipt 
M@emeity 1S applied to the outstanding backorders heap. 
Backorders are removed from the heap and filled until all the 
backorders were filled or the receipt quantity is exhausted. 
If all backorders are filled, the remaining receipt quantity 
is added to the current on-hand inventory. 

Cc. Material Issue 

If a demand is generated in the Demand Observation 
Generation routine for the current week and the current on- 
hand inventory is sufficient to meet the requirement, then 
material is issued and the on-hand inventory 1s decreased by 


the amount of the demand. When the requirement 1S greater 


than current on hand inventory, a backorder is created for the 
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amount of the requirement in excess of Current On-fand 
inventory. The backorder is inserted into the outstanding 
backorder heap, a FIFO priority heap [Ref. 29:p. 149], based 
on the date at which the backorder occurred. 
dad. Material Order 

At the end of each week the inventory position 
(IP) 1S examined to determine if a reorder 1s necessary [Ref. 
19:p. 3.24/25].' If IP is less than or equal to the reorder 
point (RO) then a reorder is placed. An RO is calculated for 
each quarter in the Levels routine prior to making the weekly 


calls to the SDR routine. The reorder quantity (ROQ) equals: 


ROQ = EOQ+RO+BO-OH-OS 4.2 
Where: 

IP = OH + OS - BO 

EFOQ = economic order quantity for current quarter, 
based on output from the Levels routine. 

RO = cree@rdenmocimes: 

BO = total backorders outstanding at the end of the 
current week. 

OH = total on hand inventory at the end of the 


current week. 
OS = total quantity of material on order at the end 
of the current week. 


A random procurement leadtime 1s generated at the 
time of reorder and a receipt date equal to the current date 


plus this generated procurement leadtime 1s assigned to the 


‘SDR is currently run somewhat less frequently and less 
regularly than once a week at the Navy Inventory Control 
Points. 
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reorder. The reorder is then inserted into the outstanding 
reorder heap. The random procurement leadtime is based on a 
MOmmal distribution with mean of eight quarters and variance 
of 64 quarters. The actual procurement leadtime used 1s 
constrained to a maximum of 14 quarters and a minimum of two 


quarters. 


B. SIMULATION SET-UP 
ar System Parameters 

The UICP model system parameters and their default 
settings are displayed in Table 4. The default values are the 
same as those used in the UICP, Computation and Research 
Evaluation System (CARES-D56) [Ref. 30].! Although any of 
these parameters may be changed during initialization of the 
Simulation, the default CARES values were used for Total Cost 
Analysis, Constant Demand Analysis, and Declining Demand 
Analysis simulations. The capability to change these default 
values was used in the Sensitivity Analysis simulations. 


TABLE 4. SYSTEM PARAMETERS 


Probability Break Point: 0 

fame RiSk(Prob of a stockout): 0.10 

Max Risk(Prob of a stockout): 0.35 

Shelf Life Code: 0 

Order Cost Rate: 400.00:S/order 
Obsolescence Rate: 0.12:S/unit-yr 
mie Price: P50 MOO. > tii 


MepeS 1S an application designed to provide ICP 
management with a tool to analyze and evaluate alternative 
inventory management policies prior to their implementation in 
erCe. 
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Time Preference Rate: OT07 a7 ye 

Salvage Rate: 0.02:-6/unie wares 
Storage Rate: 0.01:S/unit=-yr 
Procurement LeadTime: 8,00: @ers 

Shortage Cost: LHOOOR00:S/ Unser 
Military Essential: 0.50 

Requisition Size: 1:unit/requismemwem 


2. Random Number Seeds 

As discussed in Chapter IV.A.1 there 1s an array of 
20,000 seeds available to seed the random number generator for 
each replication of a simulation. During the initialization 
of the simulation any series of seeds in the array equal to 
the number of replications can be chosen. For example, ina 
100 replication simulation, the series of seeds from 1 to 100, 
900 to 999 or 10001 to 10100 can be specified, as long as the 
Starting seed position in the array is less than or equal to 
20,000 minus the number of replications for the simulation. 
The purpose of this feature is to allow for generation of 
dependent or independent output samples from two or more 
Simulations. The importance of this feature 1S theese 
affects the type of statistical test which may be performed 
when comparing the output from two or more simulations. 

For this thesis, dependent output samples were created 
for all simulations run within Gach ~setting. This was 
accomplished by specifying the same series of seeds for demand 
Stream generation for each simulation in a setting. Using 
dependent demand streams for performance comparisons allows 


for the comparison of the models in a similar demand 


BZ 


environment. However, the analysis must be done using a 
statistical test for dependent samples such as the paired 
difference t-test. If independent samples are desired, each 
Simulation would have to be run uSing a unique series of 


seeds. 


3. Number of Replications 
[ierorger tO obtain reasonable precision in the 
confidence intervals for the estimates of ACWT and total cost, 
the absolute error method [Ref. 23:p. 536] was used to 
determine the total number of replications to run. By using 
the absolute error method with a simulation run consisting of 
400 replications, absolute errors were obtained of no more 
than 20% of the true mean ACWT and no more than 7.5% of the 
true mean total cost with a probability of 0.95. Based on 
these results, 500 replications were used in all simulations. 
This yielded an absolute error of no more than 15% for the 
true mean ACWT and no more than 5% for the true mean Total 
Cost with a probability of 0.95. Although the error for ACWT 
may appear rather high, the error, when measured in days, was 
typically less than two days. 
4. Initial Conditions Warm-up Period for Declining Demand 
Analysis 
Mmnerent in the Simulation of a stochastic process is 


the initial transient or the start-up problem. The difficulty 


oe 


1s in determining the warm-up period for a model. The warm-up 
period covers the time it takes for the means of the random 
variables being measured in a simulation to converge to their 
Steady state values. 

We employed the "graphical procedure" that is due to 
Welch [Ref. 23:p.544] to identify when the simulation 
approached steady state. The Welch procedure is applied to 
each demand scenario. The Welch graphs (Appendix E, Graphs 7 
- 12) were generated from data that was obtained from a 100 
replication, 80 quarter simulation. The steady state random 
variable shown in the graphs is the investment (measured in 
units) in a given quarter, averaged over all replications. 
Investment in this case is the number of units on-hand plus 
the number of units in outstanding orders at the end of a 
Quarter. Investment was chosen because it most accurately 
reflects the balance between material issuing and ordering and 
when the inventory system has reached equilibrium or steady 
state. Based on Graphs 7 - 12 in Appendix E, it was 
determined that the simulated model reaches steady state with 
respect to investment by quarter 30 at the latest for all 
demand scenarios. 

The amount of time the random variable's mean remains 
in a transient state is affected by the initial conditions of 
the simulation. In an effort to reduce the warm-up period, 
the following logic was used to determine the initial on hand 


Quantity, and to schedule receipt dates and quantities for 
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Meeroders Outstanding at the start of the simulation. The 
Initial quantity of on hand inventory is set equal to EOOQO 
Givided by 2 plus safety stock [Ref. 17:p. 275]. Safety stock 
is set equal to the reorder point minus the forecasted 
leadtime demand [Ref. 19:Chap. 3]. The number of reorders 
outstanding at the start of the simulation is set equal to the 
expected number of reorders outstanding at any instant of time 
for the deterministic setting. This number equals the 
procurement leadtime divided by a reorder interval (uSing a 
0.5 rounding rule), where a reorder interval equals the EOOQ 
Givided by the forecasted quarterly demand [Ref 31:p. 93]. 
For all simulations the EOQ, reorder point, and forecast for 
Quarter one is used to calculate these initial conditions. 
The receipt dates of the reorders outstanding are uniformly 
distributed from simulation time zero to simulation time zero 
plus one procurement leadtime, and the quantity of each 


reorder outstanding was set equal to the EOO for quarter one. 
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ae SIMULATION RESULTS 


A. OVERVIEW 

This chapter will discuss the simulation results from the 
Total Cost, Constant Demand and Declining Demand Analysis. 
Total cost curves generated from the Total Cost Analysis are 
presented in Appendix E, Graphs 13-24. The simulation results 
and MADM analysis from the Constant Demand Analysis and the 
Declining Demand Analysis are presented in Appendices A and B, 
respectively. The remainder of this chapter will discuss the 
general results of each Analysis based on the goals of the 
Analysis. In addition, specific observations which deserve 


further analysis will be examined. 


B. TOTAL COST ANALYSIS 

The goal of this particular analysis was to determine if 
a minimum Total Cost (TC) associated with a single retention 
level ( symbolized by t, or RL ) existed in a stochastic 
demand environment as Tersine showed for the deterministic 
case [{Ref. 14]. Assuming a minimum TC exists, an optimal 
retention level (t,.) for each demand scenario in the Total 
Cost Analysis setting was determined that minimizes the 
respective TC. 

The results of the Total Cost Analysis simulations show 


that the TC curve for each demand scenario simulateaquwewee 
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Darabola (Appendix E, Graphs 13 to 24). While the high unit 
price demand scenario TC curves had an easily identifiable 
Siemum point, the low unit price demand scenario TC curves 
tended to be flat in the vicinity of the minimum. This 
indicates that for the low unit price settings there may be a 
range of retention levels that yield statistically equivalent 
mmm COtal costs. In addition, finding the best t, for the 
Tow unit price settings may involve other MOEs such as ACWT. 

Although all the total cost curves for the demand 
scenarios Simulated are parabolas, an interesting 
characteristic in the TC curve for the LDLP demand scenario 
can be observed (Appendix E, Graphs 18 and 24). There is a 
Memeo in the TC curve and specifically in the Total Order 
Cost curve at a retention level of approximately 3.5 years 
annual demand. The initial inventory position (IP) at time 
zero after disposal of excess inventory, for a retention level 
less than 3.5 years, was below the time zero reorde- point 
(RO) (the RO is depicted by the vertical line in G ne etc 
een 24). This caused an additional reorder to be placed 
during the simulation period for all retention levels less 
than 3.5 years. The "step" down in the total order cost curve 
occurred after the retention level exceeded 3.5 years because 
an additional reorder was not placed at time zero. The 
magnitude of the "step" down was due to the high 
administrative order cost ($850/order) in relation to the low 


unit price ($20/unit) and low mean quarterly demand (2 units/gtr). 
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Since the TC curves were parabolas, the next step in the 
Total Cost Analysis was to determine the respective optimal 
retention level (t,) that minimized TC for each demand 
scenario in the Total Cost Analysis settings. FOr elms 
analysis the optimal retention level was defined as the 
arithmetic mean of the retention levels which resulted in the 
minimum total cost for each of the 500 replications omen 
respective demand scenario simulation. The optimal inventory 


level t,, was calculated as follows: 


t¢] 
db: = il 





E = 1=1 
; n 
Where: 
1 = index for a replication of a simulation. 
n = total number of replications of a simulapaem 
t, = retention level which resulted in the minimum TC 


for a specific replication of a simulation. 
The t, values are presented in Table 5 under Alternative A. 
The t, values represent years worth of demand at the 
forecasted annual demand rate. 

In order to test the sensitivity of t, to different 
initial inventory amounts, the simulations for the Total Cost 
Analysis settings were rerun with an initial inventory of 75 
years worth of annual demand. The results of these 
Simulations are shown in Table 5 under Alternative B. The 
results presented in Table 5 indicate that t,1s very robust 


with respect to initial inventory. 


Bye. 


TABLE 5. TOTAL COST ANALYSIS OPTIMAL RETENTION LEVELS 






ALTERNATIVE 


Demand 


HDHVLE CP 


HDLVLP +90.25 +0.48 


LDHP i 
op [ass | 20.04 | as.e | 20.88 


1S a 95% aonndanee Tatecal iG ales, 

To summarize, the initial results indicate that a t, 
exists for each demand scenario simulated, and the value of t, 
varies considerably with respect to unit price, mean quarterly 
demand and variance of demand. The following correlation 
between t, and unit price, mean quarterly demand and variance 
of demand in a stochastic environment can be developed. As 
unit price increases t, decreases, aS mean quarterly demand 
increases t, decreases, and as variance of demand increases t, 


macreases. 


C. CONSTANT DEMAND ANALYSIS 

The goal of this analysis was to observe the performance 
of the various proposed models under the same conditions used 
in the Total Cost Analysis. We hoped to draw some conclusions 
about the performance of these models in ae stochastic 


environment by comparing the performance of the models to the 


SI, 


appropriate optimal retention levels (t,) obtained from the 
Total Cost Analysis. 

Simulation and performance comparison results are 
presented in Appendix A. ACWT and TC values that appear in 
bold print in Appendix A indicate these values are 
Statistically equal to or less than the respective optimal 
value, based on the paired difference t-tests conducted in the 
performance comparison. 

Table 6 summarizes the results of the performance 
comparison. The table is designed to be a decision tool to 
assist in determining which models might be appropriate for a 
specific demand scenario with respect to the relative weight 
that management places on the MOEs of TC and ACWT. Entries in 
Table 6 indicate which models were the best performers for a 
specific combination of demand scenario and MOE weighting. 


TABLE 6. CONSTANT DEMAND ANALYSIS SUMMARY RESULTS 
T DECISION TABLE 


DEMAND SCENARIO MOE WEIGHTING 


Mean Demand | Unit Total |Mostly}| Equal |Mostly| ACWT 
Demand| Variance] Price |Cost (TC) nie TC /AGCWT. lee CWE 
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While no single model's RL consistently matched the 
optimal retention level, the NB-MOD model performed the best 
across all demand scenarios. Additionally, there was 
typically at least one model's RL which matched the optimal 
for each demand scenario. 

The RL for the TRAD model remained constant for all demand 
scenarios because mean quarterly demand, unit price, and 
demand variance are not parameters in the calculation of the 
TRAD model's RL. The RLs for the "net benefit" models as a 
group behaved the same as the optimal with respect to changes 
in mean quarterly demand and unit price as discussed in the 
Total Cost Analysis results. Changes in demand variance had 
little effect on the RLs of the "net benefit" models, most 
likely because demand was assumed to be deterministic in the 
derivation of the basic net benefit equation. 

The following general observations can be made from the 
performance comparison results. Based solely on TC, there was 
usually one model which obtained the true optimal solution. 
The only exception was for the HDLVHP demand scenario in which 
no model had a TC which was statistically equal to the true 
optimal solution. This can most likely be explained by the 
fact that the total cost curve for the HDLVHP demand scenario 
(Appendix E, Graph 14) has the most distinct minimum point on 
its curve as compared to the other demand scenario total cost 
Gurves. This argument is also supported by the fact that the 


Confidence interval about the optimal retention level for the 
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HDLVHP demand scenario is the smaller than the confidence 
intervals of the other demand scenario optimal retention 
levels (Chapter V, Table 5). 

When taking into account ACWT and TC there were generally 
several models which performed as well as or better than the 
optimal, with the NB-MOD model being the most consistent top 
performer. The TRAD model consistently had a higher RL and 
was the best performer with respect to ACWT for all demand 
scenarios except HDLVLP and LDLP. For the latter two demand 
scenarios the difference between all the models' respective 
ACWTs' was insignificant. 

It 1S interesting to note that under the HDHVLP and LDLP 
demand scenarios the TRAD and NB-NPV models had lower average 
total costs than the respective optimal solution. The lower 
TC for the two models could be expected due to the fact that 
both the HDHVLP and the LDLP TC curves (Appendix E, Graphs 15 
and 18) from the Total Cost Analysis were flat in the vicinity 
of the minimum TC point on the curve. After further analysis 
it was determined that the calculated optimal retention level 
for the HDHVLP and the LDLP demand scenarios may vary 
depending on how optimality was defined in the Total Cost 
Analysis. In light of the HDHVLP and LDLP results an 
alternative definition of the optimal retention quantity was 
developed. 

In the Total Cost Analysis the optimal retention level, t, 


for each demand scenario in Chapter V Table 5 (Alternate A) 
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was defined as the arithmetic mean of the retention levels 
which resulted in the minimum total cost for each of the 500 
replications of the respective demand scenario simulation. 
Tae revised optimal retention level (t*) was defined as the 
retention level associated with the arithmetic mean of the 

fiomnem total costs of all the replications of the respective 
demand scenario simulation. The revised optimal retention 


Heyelt* was calculated as follows: 


5 





aye, 
Ces 
i = 1=1 
Ec n 
t* = argmin C, She 
Ee I 
Where: 
a = the average TC for a specific retention level 


across all replications of a simulation. 


c.; = the TC for a specific retention level anda 
specific replication of a simulation. 

t =a specific retention level simulated. 

T = the set of all retention levels simulated (0.0, 
G) 3S) SE 0) eo , mM) 

fie—etiitial On Nand inventory prior to disposal. 

Mmmme=eiincaex for a replication of a simulation. 

mM = total number of replications of a simulation. 


Table 7 presents the t, and t* values for all demand 
Streams. The values for t* tended to be greater for the HDHVLP 
and LDLP demand scenario, and were also closer to the 


respective retention levels obtained from the TRAD and NB-NPV 
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models than to the respective values for t,. For the HDHVLP 
demand scenario this quantity was 13 years and for the LDLP 
demand scenario this quantity was 17 years. It should be noted 
that the differences between the respective t* for the 
remaining demand scenarios and the optimal t, were not 


statistically, signi leone 


TABLE 7. OPTIMAL RETENTION LEVELS CALCULATION ANALYSIS 





(C.I. 18 a 95% confidence interval) 


The difference between t, and t* for the HDHVLP and LDLP 
demand scenarios can be attributed to backorders which 
occurred when the Total Cost Analysis optimal quantity, to, 
was retained and which did not occur when the t* quantity was 
retained. The backorders occurred in approximately 10% to 15% 
of the replications of the Constant Demand Analysis 


Simulations due to large spikes in observed demand between 
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Quarters 30 and 55. However, the extra stock held when t* was 
retained was sufficient to satisfy this increased demand. 
Because the two demand scenarios were low unit price 
(S20/unit) scenarios, the high shortage cost ($1500/unit year 
of shortage) tended to dominate TC. Therefore when these 
backorders occurred, the TC for the t, retention level 
increased by 120% to 150% and was significantly higher than 
the TC for the t* retention level. This tended to force the 
Simulation average minimum TC out to t*. 

It should be noted that for 85% to 90% of the Constant 
Demand Analysis simulation replications the t, retention level 
resulted in the minimum TC. Additionally, over an entire 
Simulation the average total costs for the HDHVLP and LDLP 
demand scenarios and the TRAD and NB-NPV models, respectively, 
were statistically equal to the respective average optimal 
total cost based on the t, retention level. 

fomcummary, it 1S difficult to conclude whether t, or t* 
better defines the optimal retention quantity for the HDHVLP 
and LDLP demand scenarios. Although there is a significant 
difference between t' and t* for the HDHVLP and LDLP demand 
Scenarios, the average total costs which result from the two 


retention levels are statistically equivalent. 


D. DECLINING DEMAND ANALYSIS 
The goal of this analysis was to compare the models ina 


scenario that involved declining mean quarterly demand. For 


SS 


this analysis, simulation and performance comparison results 
are presented in Appendix B. ACWT and TC values that appear 
in bold print in Appendix B indicate the values which were the 
best performers from among the five models. When more than 
one value 1S in bold print this indicates that the values were 
Statistically equivalent based on the paired difference t- 
GbEeSiSs2 

The values for TC and ACWT shown in Appendix B were 
accumulated over quarters 30 through 66 in the respective 
Declining Demand Analysis simulations. Data for TC and ACWT 
was originally collected for the full 66 quarters of each 
Declining Demand Analysis simulation. The results using the 
full 66 quarters of data were significantly affected by the TC 
and ACWT data collected during quarters 1 through 29 when mean 
quarterly demand was constant. In general, the results 
showed that the performance of all of the models was 
Statistically equal when the full 66 quarters of data were 
used. Therefore, in order to get a more accurate picture of 
the effect each model's RL had on the its TC and ACWT during 
the declining demand period, data for the performance 
comparison was collected for quarters 30 through 66 only. 

Table 8 summarizes the results of the performance 
comparison. The table is designed to be a decision tool to 
assist in determining which models might be appropriate for a 
specific demand scenario with respect to the relative weight 


management places on the MOEs of TC and ACWT. Entries in 
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Table 8 indicates which models were the best performers for a 
specific combination of demand scenario, pattern of declining 
demand and MOE weighting. 


TABLE 8. DECLINING DEMAND ANALYSIS SUMMARY RESULTS 
AND DECISION TABLE 


DEMAND SCENARIO | MOE WEIGHTING 


Mean jewels Mostly | Equal | Mostly | ACWT 
Demand Price TC/ACWT| ACWT 


High 
























Demand Decline 









Variance Pattern 





soy ea al ie le Nal 
concave pre eee ee eee 









eee | = TRAD, 2 = NB, 3 = NB-MOD, 4 


The following general observations can be made from the 


results of the performance comparison. No one model dominated 


67 


across all demand scenarios based on TC alone. For the 
"mostly TC" and "mostly ACWT" categories of management 
emphasis, the NB-MOD and the NB-NPV models were consistently 
top performers regardless of demand scenario and pattern of 
declining demand. For the "only TC" category of management 
emphasis, the NB-NPV model was consistently a top performer 
regardless of demand scenario and decline pattern. Similar to 
the correlation seen in the Total Cost Analysis between the 
changes in the RL and changes in demand, the RLs for the "net 
benefit" models increased as demand decreased during the 
Simulation's period of declining mean quarterly demand. The 
increases were most apparent for the low unit price scenarios. 
Because the RLs for the "Net Benefit" models were changing 
throughout the Declining Demand Analysis simulations, the 
retention levels shown in the Declining Demand Analysis 
results (Appendix B) represent the average RL over quarters 30 
through 66. Graphical illustrations of the change in the RLs 
for all of the demand scenarios and patterns of declining 
demand are shown in Appendix E, Graphs 25 to 42. 

There are several noticeable effects on the RL 
calculations made during periods of declining demand, using 
the "net benefit" models. The effects can be attributed to 
the demand forecasting method used in UICP and the use of the 
forecasted demand in the RL calculations. First, there is a 
lag between the time the declining demand period starts and 


the time the RL reacts to the changing demand. This lag is 
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directly correlated to the lag between the time the actual 
demand changes and the time the forecasted demand reflects 
this change. 

Second, the step-ups in RLs for the demand scenarios with 
high quarterly mean demand (Graphs 28 to 33 and 37 to 42) 
occurred when a "trend" (declining demand) was detected by the 
UICP demand forecasting application. When a "trend" is 
detected, demand forecasting switches from simple exponential 
smoothing to a four quarter moving average. This change in 
forecasting method caused the forecasted demand, reorder 
Seemeity (FOO) and reorder point to drop rapidly, which in 
turn resulted in the step increases in the RLs. The step is 
more prominent in the demand scenarios with a convex pattern 
of declining demand. This as due to the fact that the 
decrease in demand was more rapid for the convex pattern of 
declining demand and the final forecasted quarterly demand was 
approximately one unit per quarter less than the concave and 
Step patterns of declining demand. 

Third, the steps down in the RLs for the demand scenarios 
with low mean quarterly demand and high unit price (Graphs 25 
to 27) occurred when actual demand approached zero at the end 
of the declining demand period and the forecasted demand had 
not yet stabilized. For some simulation replications, several 
Quarters of zero demand, in sequence, were observed when 
actual mean quarterly demand was close to zero after the 


period of declining demand. For these replications and 
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quarters this caused the forecasted demand and the RLs to go 
bQssZerGr Therefore, the simulation average RLs for those 
quarters were lower than the average RLs for the remaining 
quarters. When the demand forecast stabilized about the final 
mean quarterly demand, the RLs also stabilized. 

Finally, the RLs for the NB-MOD model in the demand 
scenarios with high mean quarterly demand and high unit price 
did not increase as expected when demand decreased (Graphs 28 
to 33). This can be attributed to the decrease in expected 
number of shortages as demand decreased. The NB-MOD model RL 
(Equation 2.24) is a function of the NB model RL (Equation 
2.13) plus a term added to account for potential shortages. 
As seen in Graphs 28 through 33 the NB model RLS were 
increasing as demand decreased. Because the NB-MOD model RLs 
are decreasing in these same scenarios, this indicates that 
the increase in the RLS due to the decrease in demand was more 
than offset by the reduction in the RLs due to the decrease in 
expected number of shortages. 

A specific observation which warrants further discussion 
is the effect that the five unit minimum Retention Quantity 
(RO) constraint (used in the UICP retention logic) has on the 
results of simulations involving low mean quarterly demand. 
The Declining Demand Analysis simulations were originally run 
with only the UICP model constrained to a minimum RQ of five 
units. As a result, when forecasted annual demand approached 


zero at the end of the declining mean quarterly demand period, 
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the UICP RQ remained fixed at five units while the 
unconstrained RQs for all of the mathematical models 
approached zero. In essence, without the constraint the 
mathematical models' RO stayed at zero regardless of how large 
the respective RLS were. Additionally, while the UICP RQ 
remained a five units, the RL grew substantially. Based on 
preliminary results it became apparent that the five unit 
minimum retention quantity gave the UICP a significant 
advantage over the other models with regard to total cost and 
average customer wait time. The five unit minimum retention 
Quantity was then applied to all the models and the Declining 
Demand Analysis Simulations were rerun to determine what 
effect this constraint would have. We found that this minimum 
retention quantity improved the performance in both the TC and 
ACWT MOEs for all of the models and these results were used to 
make the final performance comparison presented in Appendix B 


and Table 8. 
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vi. SENSITIVITY ANALYSIS 


A. OVERVIEW 

The sensitivity analysis was designed to determine how 
changes in selected parameter values affect the retention 
levels of the respective models. The parameters used in this 
analysis were chosen because it is extremely difficult to 
accurately estimate the parameter values from available 
historical Gostes The estimates for these rates could be 
somewhat inaccurate because the historical costs associated 
with a given parameter are either not available or not easily 
allocated to the individual items. Therefore, it 1S important 
to determine how each model reacts to changes in these rates. 
The goal of the sensitivity analysis is to identify which 
model's RL calculations are robust with respect to changes in 
the various parameter values. This information she@ulda yaw 
decision makers in the selection of an appropriate model based 
on the level of uncertainty in the value of a specific 
parameter. In addition to the robustness of the RL's of the 
models based on changes in a given parameter, we will also 
look at the robustness of the model's performance, with 
respect to TC and ACWT for four specific scenarios from the 
Declining Demand Analysis. 


The sensitivity analysis was conducted for two demand 
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scenarios (HDHVHP and LDHP) and two declining demand patterns 
Peom the Declining Demand Analysis (convex and concave). For 
each combination of demand scenario and declining demand 
pattern, four parameters were analyzed. For each parameter 
meee values (including the UICP (CARES) default rates used in 
the Declining Demand Analysis) were used. Table 9 summarizes 
memo simulation settings which resulted from combinations of 
demand scenario, declining demand pattern and parameter 
values. For a specific setting all other parameters and 
Simulation characteristics were identical to those used in the 
Declining Demand Analysis for the respective demand scenario 


and declining demand pattern. 


TABLE 9. 16 SENSITIVITY ANALYSIS SIMULATION SETTINGS 











HOLDING COST ORDERS COs T 
RATE RATE 


SALVAGE 
RATE 






OBSOLESCENCE 
RATE 


0.06 0.01 0.01* 200 
S/UNIT-YR %/UNIT COST S/UNIT-YR S/ORDER 
0.09 0.02* 0.03 400 
S/UNIT-YR %/UNIT COST S/UNIT-YR S/ORDER 
0.12* 0.05 0.05 800* 
S/UNIT-YR %/UNIT COST S/UNIT-YR S/ORDER 
0.15 0.15 0.07 1200 
S/UNIT-YR %/UNIT COST S/UNIT-YR S/ORDER 


(* Denotes UILCP(CARES) default value) 







B. RESULTS 
Simulation and performance comparison results are 


presented in Appendix C. The ACWT and total cost in bold 


1S. 


print indicate the value which is the best performer in its 
respective MOE category. When more than one value is in bold 
print this indicates that the values were statistically 
equivalent based on the paired difference t-test. Table 10 
and Table 11 summarize the effects the varying rates had on 
each model's RL for the HDHVHP demand scenario and the LDHP 
demand scenario, respectively. 

In general, based on the results displayed in Tables 10 
and 11 the following observations can be made with regards to 
the sensitivity of the RL's of the models to changes in a 
given parameter. All models were robust with respect to 
changes in order cost rate and the three "net benefit" models 
were robust with respect to changes in the holding cost rate. 
The TRAD model was sensitive to changes in holding cost rate 
and all models showed sensitivity to changes in obsolescence 
rate. The type of demand scenario had little effect on the 
RL's for all of the models. 

Observations regarding the sensitivity of the models due 
to changes in a given parameter value are summarized in Table 
oe The observations in Table 12 indicate the effect of 
changes in a given parameter value for a specific demand 
scenario and pattern of declining demand on the performance of 
the various models. For each parameter, the respective UICP 
(CARES) default parameter value was used as the comparison 
baseline. The following types of observations were made. 


Observation type 0 means no significant change occurred ina 
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TABLE 10. RANGE OF AVERAGE RL - HDHVHP SCENARIO 


Convex 







Scenic 





Cost 


| Concave 


Convex 





} Concave 


Obsolete Convex Low 13.9 7.6 Seer roar 


Concave 









Concave 


Low 1 ey 5.4 |} 6.9 Sy # 8.0 
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TABLE 11. RANGE OF AVERAGE RL - LDHP SCENARIO 


Holding Convex 


Cost 


Convex 


Obsolete Convex 


Convex 


| Concave 
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model's performance. Observation type 1 occurred when a 
model's performance improved for parameter values greater than 
the respective UICP (CARES) default parameter value. 
Observation type 2 occurred when a model's’ performance 
Improved for parameter values less than the respective UICP 
(CARES) default parameter value. Observation type 3 occurred 
when a model's performance declined for parameter values 
greater than the respective UICP (CARES) default parameter 
value. Observation type 4 occurred when a model's performance 
declined for parameter values less than the respective UICP 
(CARES) default parameter value. 

Based on the results displayed in Table 12 the following 
general observations with regards to the sensitivity can be 
made. The performance of the NB and NB-MOD models was robust 
with respect to changes in all parameter values for all 
scenarios. The performance of the UICP model was sensitive to 
changes in all parameters values, except salvage rate, for all 
LDHP scenarios. The performance of the TRAD model tended to 
improve with both increases and decreases in the obsolescence 
rate and salvage rate parameter values for all HDHVHP 
scenarios. The NB-NPV model's performance tended to decline 
for salvage rate parameter values greater than the UICP 


(CARES) default value in both the LDHP and HDHVHP scenarios. 
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TABLE 12. SENSITIVITY ANALYSIS PERFORMANCE OBSERVATIONS 


Holding wow fo fo fs fof 
one = ole 


Convex LDHP — 


Concave 


3 
fpwvap to to to to fo 


Obsolete Convex 


Concave 


Concave 





The sensitivity analysis can be summarized as follows. 
Although the RL for the TRAD model displayed the most 
Sensitivity to changes in the parameter values analyzed, it 
had little effect on the performance of the TRAD model as 
compared to all other models analyzed. The UICP model 
performance displayed the most sensitivity to changes in the 


parameter values analyzed. 
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Mar. . OVERVIEW, CONCLUSION AND RECOMMENDATIONS 


A. OVERVIEW 

This thesis evaluated the effectiveness of the Navy's UICP 
economic retention model. The evaluation was performed by 
comparing several mathematical economic retention models with 
the Navy's retention model. There were two primary factors 
that motivated this thesis. First, the Navy does not 
currently apply economic retention theory when making 
retention decisions for the majority of the material managed 
by the Navy. Second, the excess inventory problem will 
continue to grow as the Navy's budget and fleet are further 
reduced. 

An analysis of the models was performed for a variety of 
demand scenarios in both steady state and declining demand 
Situations. The analysis was designed with two goals in mind. 
The first goal was to determine which model(s) were most 
effective in a demand environment similar to the Navy's 
stochastic demand environment. The second goal was to 
evaluate how the Navy's retention process performed with 
respect to the mathematical models. 

A simulation of the Navy's UICP demand process and the 


mathematical retention models was developed. The evaluation 


AS 


of the various models was based on the measures of 
effectiveness (MOE) of total cost (TC) over a specified period 
of simulation time and average customer wait time (ACWT) per 
requisition for all requisitions generated over a specified 
period of Simulation time. The research also examined model 
sensitivity to changes in various parameters common to the 
models. The parameters were chosen for the analysis because 
UICP uses estimates of the true rates and these estimates 
could vary considerably from the true rates. Results of the 
sensitivity analysis helped to determine the practicality of 


applying the models in the UICP environment. 


B. CONCLUSION 

The findings of this research showed that, of the models 
analyzed, there was not one economic retention model or 
retention quantity which yielded the lowest total cost and 
ACWT for all of the demand and retention scenarios analyzed. 
There were two factors which contribute to this. First, the 
optimal retention level varied significantly with demand 
scenario and management weighting of the MOEs of TC and ACWT. 
Second, all the models analyzed did not account for the 
stochastic nature of demand for Navy managed items. Bue 
based on the results of all analysis, the "net benefit" 
models, as a group, performed the best and generally performed 
better than the UICP retention model. Additionally, for most 


demand scenarios in both the Constant and Declining Demand 
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Analysis, the decision on which model to chose could typically 
Mmemaetermined by total cost alone. This was due to the fact 
that the difference in the models' ACWTs (measured in days) 
for each demand scenario, were generally small. 

The results of the Total Cost Analysis showed that there 
was a unique "optimal" retention level for a given demand 
scenario in a stochastic demand environment. It also showed 
that the “optimal" retention level varies significantly with 
changes in unit price, mean quarterly demand and variance of 
mean quarterly demand. 

The Constant Demand Analysis compared the models to the 
"Optimal" retention level determined in the Tttal Cost 
Analysis. In general, when considering both TC and ACWT the 
Mathematical models performed well in the Navy's stochastic 
demand environment with respect to the performance obtained 
from the "optimal" retention level. Additionally, there was 
typically at least one model which performed as well as the 
"Optimal" retention level with respect to TC alone. The NB 
and NB-MOD models consistently outperformed the UICP model 
when management emphasis was placed on total cost or mostly on 
meeal cost. 

Meerresults of the Declining Demand Analysis indicated 
meee the "net benefit" models, as a group, were the best 
performers over all scenarios and typically outperformed the 
UICP retention model. The average retention quantities of the 


best performers in the Declining Demand Analysis varied with 
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changes in the unit price, mean quarterly demand and the 
variance of mean quarterly demand ina pattern similar to that 
observed in the Total Cost Analysis for the "optimal" 
retention level. The declining demand pattern had little 
effect on overall model performance. 

The performance of the TRAD model dominated the 
performance of the other models across all analysis scenarios 
with respect to ACWT. But the performance of the NB, NB-MOD 
and UICP models was competitive with respect to ACWT Giemese 
of the Declining Demand Analysis scenarios. It 1S important 
to note that while there was generally a significant variation 
in ACWT in terms of percentage difference, in most cases the 
difference in terms of days was typically small. Thess 
observation applies to both the Constant and Declining Demand 
Analysis. 

The results from the sensitivity analysis showed that the 
performance of the "net benefit" models, aS a group, was 
robust with respect to changes in all the parameters analyzed. 
The UICP model performance showed the most sensitivity to 
parameter changes, especially with respect to the low demand 
scenarios. Although the RL for the TRAD model displayed the 
most sensitivity to changes in the parameter values analyzed, 
it had little effect on the performance of the TRAD model as 


compared to all other models analyzed. 
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C. RECOMMENDATIONS 

There are three areas related to this research which merit 
Meeener study. First, because all of the models' actual 
retention quantities are dependent upon the demand forecasting 
method, the effectiveness of a model is limited by the 
accuracy of the demand forecast. It would be interesting to 
see how performance would change if demand forecasts were 
adjusted for known changes in future demand (1.e. declining 
demand due to decommissioning of ships). Second, further 
modifications to the NB-MOD model could be made to improve the 
treatment of the stochastic nature of demand. Modifications 
could include changes in the holding cost savings and 
repurchase cost terms. The goal would be to develop a model 
which performed effectively across all demand scenarios. 
Third, the simulation developed for this thesis could be 
modified to include the Navy's repairable item demand process 
in the Forecasting, Levels and Supply/Demand Review procedures 


of the main program. 
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APPENDIX A. CONSTANT DEMAND ANALYSIS RESULTS 


HDHVHP 


CT (OPTIMAL [| TRAD | ~=NB_ | ~NB-MOD | NB-NPV | UICP | 


0.90] 0.86] BS 0.90 
75/25 p70} 0.95] G9] LT 
50/50 


HDHVLP 





| OPTIMAL | TRAD | NB | NB-MOD | NB-NPV | UICP | 
ACWT TOT 8B] SO] 498] 83]. 
MADM _% ACWT/% TC 


HDLVHP 


_ OPTIMAL | TRAD | ~=NB_ | NB-MOD | NB-NPV | UICP _ 
ACWT | 10.55] 5.87] 11.06] 9.921 10.94) 9D 
MADM _% ACWT/% TC 

75/25 | TT 92S] BB] 0BS] 0.70 


HDLVLP 


CT OPTIMAL | TRAD | NB | NB-MOD | NB-NPV | UICP | 
ACWT | 0987] 112] 9G] 1.23) 1.08 
YRS_RL 
MADM _% ACWT/% TC 

17 ee eT) ee a ae 




























pC OPTIMAL | TRAD | NB {| NB-MOD | NB-NPV | UICP | 
YRS RE | MT 13.88] 585] A] 5G] 8.00 
MADM _% ACWT/% TC 

| 5/75 | CO] BB] 8.90) | oles | 


LDLP 
PCS COPTIMAL | TRAD | NB {| NB-MOD | NB-NPV | UICP _| 
ACWT 2 OTT 20] 04] GT) 1S 


% ACWT/% ITC 
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Constant Demand Analysis Results Model Ranking by MADM Results 


ep | A | 


25% ACWT/75% TC | UICP* |NB-MOD* |OPTIMAL*| NB | NB-NPV | TRAD 
75% ACWT /25% TC | TRAD UICP | NB-MOD | OPTIMAL| NB | NB-NPV 
[50% ACWT/50% TC | TRAD UICP NB-MOD | OPTIMAL| NB | NB-NPV 


) 

: HDLVHP 

(25% ACWT / 75% TC [OPTIMAL" | NB-MOD* | NB | NB-NPV UICP TRAD 

| 75% ACWT/25% TC | TRAD UICP | NB-MOD | OPTIMAL | NB-NPV | NB | 
| 50% ACWT/50% TC | TRAD UICP | NB-MOD | OPTIMAL | NB-NPV | NB | 


HDHVLP 


| 25% ACWT/75% TC | TRAD | OPTIMAL} NB-MOD | UICP | NB | NB-NPV 
75% ACWT /25% TC | TRAD | OPTIMAL | NB-MOD UVICP | NB | NB-NPV 
50% ACWT/50% TC | TRAD PTIMAL | NB-MOD | UICP | NB _ | NB-NPV 


| UICP | 
| UICP 
|NB-MOD* 
| UICP | 
| UICP 
| OPTIMAL 
| OPTIMAL | 
HDLVLP 
25% ACWT/75% TC | TRAD* |NB-MOD* |OPTIMAL*| UICP | NB __ | NB-NPV 
75% ACWT /25% TC | TRAD | NB-MOD | OPTIMAL| UICP | NB | NB-NPV 
| 50% ACWT/50% TC | TRAD |NB-MOD | OPTIMAL| UICP | NB | NB-NPV 
i NB 
_NB-NPV _ 
| NB 
| NB 


LDHP 


25% ACWT / 75% TC | NB-MOD* OPTIMAL*|] NB-NPV | UICP TRAD 
175% ACWT /25% TC | TRAD ICP | NB | NB-MOD | OPTIMAL | NB-NPV 
50% ACWT/50% 1C | TRAD ICP | NB | NB-MOD | OPTIMAL | NB-NPV 


| LDLP 

25% ACWT/75% TC | NB-MOD | NB-NPV | OPTIMAL| TRAD | NB | UICP 
NB-NPV | OPTIMAL | TRAD UICP 
NB-NPV | OPTIMAL | TRAD UICP 


Note: * indicates models have same rank and are both ranked as 1. 


| 25% ACWT /75% TC | NB- 
75% ACWT / 25% TC | NB-MOD 
50% ACWT /50% TC | NB-MOD 


CO 
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APPENDIX B. DECLINING DEMAND ANALYSIS RESULTS 


Declining Demand Analysis Results: HDHVHP 





SSTEPDECREASES | ~=TRAD | ~=NB___| NB-MOD | NB-NPV |  UICP__| 
AVGYRS RL] 13.88] 55] 4] 94) 8.00 


25% ACWT / 75% TC 
75% ACWT / 25% TC 
50% ACWT/50%TC | 0.96} 0.88] 8B] BG] 0.91 
















CONVEXDECREASES | TRAD | NB | NB-MOD | NB-NPV | UICP | 


25% ACWT / 75% TC 
75% ACWT/25% TC | 0.99] 0.95] 0.95] 0.93] 0.98 
50% ACWT/50%TC | 0.98] 0.96 9G] S| 9B 


CONCAVE DECREASES | TRAD | NB | NB-MOD | NB-NPV | UICP | 


25% ACWT / 75% TC 
75% ACWT/ 25% TC | 0.97] 0.87] 0.90 8G] 9 
50% ACWT / 50% TC 


Model Ranking by MADM Results 
[i 28 eee eee 









ISTEPDECREASES || 
| 25% ACWT/75% TC | TRAD‘ | NB‘ | UICP* | NB-MOD | NB-NPV __ 

75% ACWT/25%TC | TRAD | UICP_ | NB-MOD | NB | NB-NPV __ 
| 50% ACWT/50% TC | TRAD | UICP | NB-MOD | NB | NB-NPV 












CONVEXDECREASES | | 
| 25% ACWT/75% TC | NB* | UICP* | NB-NPV* | NB-MOD | TRAD _ 
| 75% ACWT/25%TC | TRAD | UICP | NB | NB-MOD | NB-NPV__ 
50% ACWT/50% TC | UICP* | TRAD* | NB | NB-MOD | NB-NPV _ 
















CONCAVEDECREASES | £| #7» » (2 D2 ae 
25% ACWT/75% TC. | NB-NPV* | NB* | UICP* NB-MOD | TRAD 






75% ACWT / 25% TC TRAD UICP NB-MOD | NB | NB-NPV 
50% ACWT / 50% TC TRAD* UICP* NB-MOD | NB | NB-NPV 


Note: * indicates models have same rank and are both ranked as 1. 
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Declining Demand Analysis Results: HDLVHP 













STEPDECREASES se | ~=TRAD | ~=NB_ | NB-MOD | NB-NPV | UICP 
ACWT COT 0] .93| 0.90 | 
| 25% ACWT/75%TC | 0.96] 0.99] 0.98) 0.99 0.97 
TS% ACWT/25%TC | 0.99} 0.99) 0.99) 0.98 | 

| 97] 9} 0.99] 0.98] 0.98 | 


=| 









50% ACWT / 50% TC 













CONVEX DECREASES | TRAD | NB | NB-MOD | NB-NPV | UICP 
MADM 

0.98} 1.00/00} 1.00) 0.99 
| 75% ACWT/25%TC | 0.99] 1.00) 1.00) 0.99] 0.99 | 
50% ACWT/50%TC | 0.98) 1.00) 00] 00] 0.99 | 








CONCAVE DECREASES | TRAD | NB | NB-MOD | NB-NPV | UICP 














25% ACWT/75%TC | 0.90} 0.99 
PS%ACWT/25%TC [| 0.97] 0.99] 0.99 00] 0.99 
BO%ACWT/50%TC [| 0.93] 0.99] 0.99] 1.00 0.97 


Model Ranking by MADM Results 


a eel os se SS 


MeeEDECREASES =| 

25% ACWT/75%TC | NB*___| NB-NPV* | NB-MOD | UICP | TRAD | 
| 75% ACWT/25%TC | NB-MOD* | UICP* | —NB* | ~TRAD* | NB-NPV | 
50% ACWT/50% TC | NB* | NB-MOD* | NB-NPV | UICP | TRAD | 


avd, 










NB-NPV* UICP TRAD 
NB-NPV* UICP TRAD 















ee ae 
NB-MOD | UICP TRAD 
B-MOD | UICP TRAD 
NB-NPV B-MOD | _UICP TRAD 


Note: * indicates models have same rank and are both ranked as 1. 


50% ACWT / 50% TC 
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Declining Demand Analysis Results: HDHVLP 


STEPDECREASES | ~=TRAD {| ~=NB___ | NB-MOD | NB-NPV | UICP | 
| 7S% ACWT/25%TC | 066] 87} 1.00] 0.81] 5 
50% ACWT / 50% TC 









= 








CONVEX DECREASES | TRAD | NB | NB-MOD | NB-NPV | UICP 
MADM. 
50% ACWT/50% TC | 0.90] 0.91] 100] 84] OS 







CONCAVE DECREASES | TRAD | NB | NB-MOD | NB-NPV | UICP | 
ACWT CC CCC741 | 542 48:60) 
75% ACWT/25%TC | 0-69] 088] 00] 074] 0.62 | 









SSTEPDECREASES | | 

25% ACWT/75%TC | NB-MOD | NB | NB-NPV | TRAD | UICP 
| 75% ACWT/25% TC | NB-MOD | NB | NB-NPV | TRAD | UICP 
| 50% ACWT/50%TC | NB-MOD | NB | NB-NPV | TRAD | UICP | 












[CONVEXDECREASES || 
| 25% ACWT/75%TC | NB-MOD | NB | TRAD | NB-NPV | UICP 

75% ACWT/25%TC_ | NB-MOD | NB | TRAD: | NB-NPV | UICP 
| 50% ACWT/50%TC | NB-MOD | NB [| TRAD | NB-NPV | UICP 











CONCAVE DECREASES | | CC | eee 
| 25% ACWT/75%TC | NB-MOD | NB | NB-NPV | TRAD | UICP 
75% ACWT/25% TC | NB-MOD | NB | NB-NPV | TRAD | _UICP 

50% ACWT/50%TC | NB-MOD | NB | NB-NPV | TRAD | UICP 


Note: * indicates models have same rank and are both ranked as 1. 
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CONVEX DECREASES | TRAD | NB | NB-MOD | NB-NPV | UICP 
‘JACWT CT | 4G | 4G] 4G| 0.4 


IAVG YRS RL 13.88 12.95 15.16 10.10 


Declining Demand Analysis Results: HDLVLP 










ISTEPDECREASES _—| ~=TRAD | ~=NB__ | NB-MOD | NB-NPV | UICP 
CWT CC] OG OG |G] TE. 





25% ACWT / 75% TC 
IS%ACWT/25%TC | 0.97] 0.96} 0.95] 0.63] 1.00 
50% ACWT / 50% TC 
















MADM 
25% ACWT/75%TC | 0.99] 0.95] 0.94] 1.00] 0.96 
TS% ACWT/25%TC | 0.98] 0.98} 0.98) 1.00) 0.96 


50% ACWT / 50% 













CONCAVE DECREASES TRAD | NB _ | NB-MOD NB-NPV UICP 
OTAL COST 3375.91 3810.11 3918.66 3676.27 3078.21 







+ 






MADM 
25% ACWT/75%TC | 0.93] 0.86] 0.84] 0G] 1.00 
75% ACWT / 25% TC 
SO%ACWT/50%TC | 0.96} 0.90] 089] 048] 1.00 











MeebecREASES =| | CEC 
| 25% ACWT/75%TC | UICP | TRAD | NB | NB-MOD | NB-NPV | 
| 75% ACWT/25%TC | UICP | TRAD | NB _ | NB-MOD | NB-NPV | 
S0%ACWT/50%TC | UICP | TRAD | NB | NB-MOD | NB-NPV 



























MueeeeDECREASES | fT 
| 25% ACWT/75%TC | NB-NPV | TRAD | UICP | NB | NB-MOD | 

75% ACWT/25% TC | NB-NPV | TRAD | NB | NB-MOD | UICP 
50% ACWT/50% TC | NB-NPV | TRAD | NB | NB-MOD |  UICP 















CONCAVEDECREASES | | 
75% ACWT / 25% TC 
50% ACWT / 50% TC 


Note: * indicates models have same rank and are both ranked as 1. 
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Declining Demand Analysis Results: LDHP 






STEPDECREASES | SC TRAD | ~—ONB__ | NB-MOD | NB-NPV | UICP | 
_ 25% ACWT/75%TC | 0.97] 91 96] 90] 0.92 | 
| 75% ACWT/25%TC | 0.99] 0.74] 088] OTT 0.79 











(CONVEX DECREASES | =TRAD | ~=NB__ | NB-MOD | NB-NPV | UICP | 
MADM- 

| 75% ACWT/25%TC | 0.99} 0.93] 0.96] T0938 | 
50% ACWT / 50% TC 













(CONCAVE DECREASES | TRAD | NB | NB-MOD | NB-NPV | UICP 
| 25% ACWT/75%TC | 0.94] 0.94) 0.96] 0.93] 0.94 | 
50% ACWT/50%TC | 0.96} 0.88] 92] 8G] 90 


Model Ranking by MADM Results 


ISTEPDECREASES | CC | CC LULU eee 

25% ACWT/75%TC | TRAD | NB-MOD | UICP | =NB_ {| NB-NPV | 
| 75% ACWT/ 25% TC | TRAD | NB-MOD | UICP | NB | NB-NPV | 
| 50% ACWT/50%TC | TRAD | NB-MOD | UICP | NB | NB-NPV | 



















CONVEXDECREASES | | Ce 
| 25% ACWT/75%TC_ | NB-MOD* | TRAD* | NB |  UICP | NB-NPV__ 

75% ACWT/25% TC | TRAD | NB-MOD | UICP | =NB_ | NB-NPV_ | 
| 50% ACWT/50%TC | TRAD | NB-MOD | UICP | NB | NB-NPV | 












CONCAVE DECREASES | 
50% ACWT / 50% TC 


Note: * indicates models have same rank and are both ranked as 1. 
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EO 





Declining Demand Analysis Results: LDLP 






ISTEP DECREASES | ~=TRAD | ~=NB_ | NB-MOD | NB-NPV | UICP | 


25% ACWT / 75% TC 
IS%ACWT/25%TC | 0.69} 1.00] 1.00) 1.00) 0.51 | 
50% ACWT / 50% TC 



















CONVEX DECREASES | TRAD | NB | NB-MOD | NB-NPV | UICP | 
25% ACWT/75%TC | 0.93] 0.99} 0.99) 1.00) 081 














CONCAVE DECREASES | TRAD | NB | NB-MOD | NB-NPV | UICP | 


25% ACWT/75%TC | 0.84) 0.99) 0.99] 1.00 
75% ACWT / 25% TC 
SO%ACWT/50%TC | 0.73) 0.99) 0.99 0 . 


Model Ranking by MADM Results 
























STEPDECREASES || 
25% ACWT/75%TC | —NB* | NB-MOD* | NB-NPV* | TRAD | UICP 
75% ACWT/25%TC | —NB*___| NB-MOD* | NB-NPV* | TRAD | UICP 

50% ACWT/50% TC | _NB* | NB-MOD* | NB-NPV* | TRAD | UICP | 









MeEXDECREASES | | 

25% ACWT/75%TC | NB-NPV | =NB_ | NB-MOD | TRAD | UICP 
| _ 75% ACWT/25%TC | NB-NPV* | NB* | NB-MOD* | TRAD | UICP 
50% ACWT/50%TC | NB-NPV | NB | NB-MOD | TRAD |  UICP | 















CONCAVEDECREASES | | | a ees 
25% ACWT/75% TC | NB-NPV | NB | + -MOD 

75% ACWT/ 25% TC | NB-NPV* | NB* | NB-MOD* 
| 50% ACWT/50%TC | NB-NPV | NB | NB-MOD | TRAD | UICP | 


Note: * indicates models have same rank and are both ranked as 1. 
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APPENDIX C. SENSATIVITY ANALYSIS RESULTS 


Sensativity Analysis: HIGH DEMAND /CONVEX /STORAGE RATE 












RATE =.01 {Default setting for DDA} 
Es ele TRAD] NB] ~ONB-MOD NB-NPV| UICP 
ACWT 25.69 24.55 


| 7T5% ACWT/25%TC | 99] 0.95] 95] 930.98 | 
50% ACWT/50%TC {| 0.98} 096] 9G} S88 | 








RATE = .03 


po TRAD] NB] OUNB-MOD] —ONB-NPV| _CUICP | 
AVGYRSRE | 98 TT GTA] 462] 8.00 

TS%ACWT/25% TC | 0.99] 0.94) 093 0.92] 0.99 | 
| 50% ACWT/S50%TC | 0.98] GS S| 0.98 | 





RATE = .05 





| TRAD] CUNB] ONB-MOD[ ONB-NPV|_CUICP| 
29.27 
PT5% ACWT/25%TC | 0.99] 0.94] 7] TT 0.99 
50% ACWT/50%TC | 0.98] 9G] 9B] 4.98 | 








RATE = .07 


pT RAD] NB] ONB-MOD| —ONB-NPV| _CUICP | 
AVG YRS RL 
| 25% ACWT/75%TC | 0.96] 96] G96 | 0.97 
| 75% ACWT/25%TC | 0.95] 88093] 087] 0.99 
50% ACWT/50%TC | 0.96] 2S] 0.98 | 


Model Ranking by MADM Results 


wei | el eS 2 ee ee ee 


Rate = 0.01 {Default setting for DDA 


25% ACWT/75%1C | NB* | UICP® NB-NPV* | _NB-MOD TRAD 
| UICP | NB | NB-MOD NB-NPV 





75% ACWT / 25% TC TRAD 
50% ACWT/50% TC | UICP* | TRAD* | NB | NB-MOD NB-NPV 





Rate = 0.03 


25% ACWT / 75% TC UICP* | NB* | TRAD NB-NPV NB-MOD 
75% ACWT / 25% T TRAD* UICP* { NB | NB-MOD NB-NPV 
50% ACWT/50% TC | UICP* | TRAD* | NB | NB-MOD NB-NPV 


OO 


Rate = 0.05 


25% ACWT/75% TC | NB-MOD* | NB* |. TRAD* UICP* NB-NPV 


75% ACWT / 25% TC UICP* TRAD* NB-MOD | NB | NB-NPV 
50% ACWT / 50% TC VICP* TRAD* NB-MOD* {| NB ~~ | _ NB-NPV 
Rate = 0.07 

25% ACWT / 75% TC NB-MOD TRAD | NB | NB-NPV 


| _UICP | 
75% ACWT/25% TC | UICP | TRAD NB-MOD | NB ~~ | _ NB-NPV 
50% ACWT/50%1TC | UVICP | TRAD NB-MOD | NB  |_ NB-NPV 


Note: * indicates models have the same rank and are both ranked as 1. 





Sensativity Analysis: HIGH DEMAND /CONVEX /ORDER COST RATE 


RATE = 200 






Po TRAD NB] ONB-MOD| —ONB-NPV| CUICP | 

AVGYRSRAL | 13.8857] TT 8BL 8.00 | 

MADM 

| 75% ACWT/25%TC {|  —s0.99] 0.93] 95] 0.98 | 
SO%ACWT/50%TC | 0.98] 95] 9G] 4B 














RATE = 400 


a TRAD| ——NB| NB-MOD NB-NPV UICP 
25.94 25.74 24.55 







TOTAL COST 348303.18|  332428.02| 337065.42|  332190.22 
AVG YRS RL 






MADM 
TS%ACWT/25%TC | 0.99] 94] 9S] 92] 0.98 | 
[50% ACWT/50%TC | 0.98] 9G] G48 | 






RATE = 800 {Default setting for DDA} 


po TRAD NB] ONB-MOD| —ONB-NPV] CUICP 
MADM 

| 75% ACWT/25%TC | 0.99] 95 S| 93] 0.98 | 
|SO%ACWT/50%TC {| —s0.98f Gf 9G] S| 0.98 | 






= 











RATE = 1200 


po TRAD] NB] OUNB-MOD| —ONB-NPV]CUICP | 
AVGYRSRE | 13.88f 552] 77] 509] 8.00 


25% ACWT / 75% TC 
75% ACWT/25%TC | 0.99] 94] 95] 0.93] 0.98 
50% ACWT/50%TC | 0.98] 96] 9G] 9G] 0.98 | 





Model Ranking by MADM Results 


Rate = 200 


25% ACWT / 75% TC uIcP* | NB* | NB-MOD NB-NPV TRAD 
75% ACWT /25% TC TRAD UICP NB-MOD | NB | NB-NPV 
50% ACWT /50% TC UICP? TRAD’ NB-MOD |. NB | _NB-NPV 


Rate = 400 

25% ACWT/75%TC | NB*_ | —UICP*, | NB-MOD | NB-NPV [| TRAD 
75% ACWT/25%TC | TRAD |  UICP_ | ~NB-MOD | NB |_—NB-NPV__ 
| 50% ACWT/50%TC | UICP* | TRAD’ | NB | NB-MOD [ NB-NPV | 





Rate = 800 {Default setting for DDA 
| 25% ACWT/75%TC { NB* | UICP* | NB-NPV* | NB-MOD | TRAD 
| 75% ACWT/25%TC | TRAD | UICP | NB | NB-MOD_ | NB-NPV_ | 
50% ACWT/50% TC | UICP* | TRAD’ | NB | NB-MOD | NB-NPV | 






Rate = 1200 
25% ACWT / 75% TC UICP* NB-NPV* NB-MOD TRAD 


| NB* 
75% ACWT / 25% TC TRAD | UICP | NB-MOD { #4=.\NB___|_ NB-NPV 
50% ACWT /50% TC UICP* TRAD* | NB | NB-MOD NB-NPV 


Note: ° indicates models have the same rank and are both ranked as 1. 
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Sensativity Analysis: HIGH DEMAND /CONVEX /OBSOLESENCE RATE 


RATE = .06 


TRAD] —sSY*NB NB-MOD NB-NPV UICP 
14.97 18.46 19.85 17.99 
321842.45|  312977.96| 314537.71 312285.90|  313143.04 


25% ACWT/75%TC | 0.98] 0.95] 0.96] 94]. | 
75% ACWT/25%TC | 0.99] 0.86] 90] B28 
50% ACWT/50%TC | 99] 90f 93} BB]. 





RATE =.09 













po TRADT NB] ONB-MOD} ——ONB-NPV| _CUICP | 
| 0.99 | : | 0.95 | 







75% ACWT / 25% TC 


50% ACWT /50% TC 






ARATE .12 {Default setting used in DDA} 


| TRAD] NB] ~—ONB-MOD| —ONB-NPV] OO CUICP 
75% ACWT/25%TC |  —Sss.99] 95] SB] 0.98 | 
| 50% ACWT/50%TC | 0.98] 96] 9G] SE] 0.988 













A 
+ 







RATE = 15 


TRAD] NB] ONB-MOD] —ONB-NPV[OUICP | 


.25% ACWT/75%TC | 0.96] 098] 97] 97] 97 | 
| 75% ACWT/25%TC | 0.99] 93.95] 92] 98 
50% ACWT/50%TC}] _—o.g7f 96] 9G] OS] 988 | 





Model Ranking by MADM Results 


Rate = 0.06 


25% ACWT / 75% TC TRAD NB-MOD VICP | NB | NB-NPV 
75% ACWT / 25% TC TRAD NB-MOD VICP | NB | NB-NPV 
TRAD vICP_ | NB | —UNB-NPV 


Rate = 0.09 


25% ACWT/75% TC | NB* | UICP* TRAD NB-MOD NB-NPV 
VICP | =—NB_ | ~NB-MOD | _NB-NPV 
50% ACWT /50% TC TRAD | NB | __UICP NB-MOD NB-NPV 






Rate = 0.12 {Default setting used in DDA 
25% ACWT/75%TC | NB* | _ ~—UICP*—| =~NB-NPV* | NB-MOD | TRAD 
| 75% ACWT/25%TC | TRAD | UICP | NB | NB-MOD | NB-NPV | 
| 50% ACWT/50%TC | UICP* | TRAD* | NB | NB-MOD | NB-NPV_ | 







Rate = 0.15 

| 25% ACWT/75%TC |  =NB* {| UICP* =| NB-MOD* | NB-NPV_ | TRAD 
| 75% ACWT/25%7TC | TRAD | UICP | NB-MOD | NB | NB-NPV_ | 
| 50% ACWT/50% TC | UICP | TRAD | NB-MOD | NB | NB-NPV_ | 


Note: * indicates models have the same rank and are both ranked as 1. 
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| 
; 
| 
| 





Sensativity Analysis: HIGH DEMAND/ CONVEX/ SALVAGE RATE 


RATE =.01 





75% ACWT / 25% TC p95] 93] 0.98 | 
50% ACWT / 50% TC po GTS] 0.98 | 
RATE = .02 {Default setting for DDA} 













25.77 
TOTAL COST 349545.12| 334089.74 333267.81 
AVG YRS RL 


25% ACWT / 75% TC 
75% ACWT/25%1C | 0.99] 95] 95] 93] 0.98 | 
SO% ACWT/50%TC {  —=s_0.98] 96] 9G] S98 












RATE = .05 













Po Cid TRADY NB] ONB-MOD] —ONB-NPV] CUICP | 

25% ACWT / 75% TC 
| 75% ACWT/25%TC | 99] 95] 7 93] 0.99 
|50%ACWT/50%TC | 0.98] 7] 7] OS] 0.99 









RATE =.15 


po TRAD] NB] ONB-MOD{—ONB-NPV| CUICP | 
| 25% ACWT/75%TC | 0.99] 97] 98] 9G] 0.99 
| 75% ACWT/25%TC | 1.00] 92] 95] BO] 0.99 

50% ACWT/50%TC{  —-—s.99] 95] 0.9 0.92] 0.99 | 


Model Ranking by MADM Results 





Rate = 0.01 
25% ACWT / 75% T UICP* NB-NPV* NB-MOD TRAD 


Ci] NB 
75% ACWT/25% TC | TRAD | VICP | =NB | NB-MOD NB-NPV 
50% ACWT /50% TC UICP* TRAD’ | NB | NB-MOD NB-NPV 


Rate = 0.02 {Default setting for DDA 

|} 25% ACWT/75%TC | NB’ {| UICP* =| NB-NPV* | NB-MOD | TRAD | 
| 75% ACWT/25%TC | TRAD | UICP | NB | NB-MOD | NB-NPV_ | 
| 50% ACWT/50%TC | UICP* | TRAD* {| NB { NB-MOD | NB-NPV_ | 





Rate = 0.05 


25% ACWT/75% 1C UICP* | NB* | NB-MOD* | NB-NPV* TRAD 
75% ACWT/25% 1C | TRAD UICP* NB-MOD NB NB-NPV 





50% ACWT / 50% TC UICP TRAD NB-MOD 


Rate = 0.15 


25% ACWT / 75% TC UICP* TRAD* NB-MOD | NB {| _ NB-NPV 
75% ACWT/25%TC | TRAD | UICP NB-MOD {| NB |_ NB-NPV 
50% ACWT / 50% TC TRAD* UICP* NB-MOD | NB | NB-NPV 


Note: ° indicates models have the same rank and are both ranked as 1. 


95 


Sensativity Ananlysis: HIGH DEMAND/ CONCAVE/ STORAGE RATE 


RATE = .01 {Default setting for DDA} 











rr =) TRAD] —s—sNB/ _—ONB-MOD NB-NPV UICP 
ACWT 12.28 172 14.04] 15.05] 13.03 
TOTAL COST 231634.28| 208435.28| 213789.73|  207017.16| 217823.08 


| 25% ACWT/75%TC | 92] 0.95) 94 95] 0.95 | 
| 75% ACWT/25%TC | 0.97} B87] 90] BG] 9 
| 50% ACWT/50% TC] 95] 2] 0.95 | 











RATE = .03 

See TRAD| _——sSNB|__sONB-MOD NB-NPV UICP 
ACWT TQ] 1644] 1547] 16.96] 14.25 
TOVAL COST 222797.49| 214975.46| 230106.07 






AVGYRS RL] BT] ATSB] 8.00 | 
| 25% ACWT/75%TC |] 0.93] 95] 49S] 0.94 | 
|) 75% ACWT/25%TC | 0.98] BG] BD] BAT 94 
|50% ACWT/50%TC | oe 0 2] OT 0.94 | 

| CUICP | 



















RATE = .05 


| TRAD] NB] OUNB-MOD|_NB-NPV. UICP 
AVGYRS RE | 787] 453] BOT 4 24t8.00 
| 25% ACWT/75%TC | 0.94] 96] 9S] E984 | 

| 2] C96 | 


50% ACWT/50%TC}] 0.95} 0.92] 984 | 


RATE = .07 


po TRAD] NB] ONB-MOD| —ONB-NPV CUICP | 
ACWT TG] 1880] 1760] 19.58) 15.18 


| 25% ACWT/75%TC | 0.91} 95] 94 940.93 | 
50% ACWT/50%TC | 0.90] 90] 91] BD] 





Model Ranking by MADM Results 


Rate = 0.01 {Default setting 
25% ACWT/75% TC | NB-NPV* [| NB* | UICP* =| NB-MOD | TRAD | 
| 75% ACWT/25%TC | TRAD | UICP | NB-MOD | NB | NB-NPV_ | 
| 50% ACWT/50%TC | TRAD* | UICP* | NB-MOD | NB | NB-NPV_ | 










Rate = 0.03 

| 25% ACWT/75%TC | NB-NPV* [|  NB* =| NB-MOD | UICP | TRAD | 
| 75% ACWT/25%TC | TRAD | UICP | NB-MOD | NB | NB-NPV_ | 
| 50% ACWT/50%TC | TRAD | UICP {| NB-MOD | NB [| NB-NPV_ | 


Rate = 0.05 

| 25% ACWT/75%TC | NB* | NB-NPV* | NB-MOD | UICP | TRAD 
| 75% ACWT/25%TC | UICP | ~=TRAD | NB-MOD | NB | NB-NPV 
50% ACWT/50% TC |  UICP. | TRAD | NB-MOD | NB | NB-NPV_ | 






Rate = 0.07 


25% ACWT/75%TC | NB | NB-NPV NB-MOD UICP 






TRAD 





75% ACWT / 25% TC UICP NB-MOD TRAD | NB |_ NB-NPV 
50% ACWT / 50% TC UICP NB-MOD | NB | TRAD NB-NPV 





Note: * indicates models have the same rank and are both ranked as 1. 
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Sensativity Analysis: HIGH DEMAND/ CONCAVE/ ORDER COST RATE 


RATE = 200 


TRAD NB-MOD NB-NPV UVICP 


| CNB 
12.28} 15.55) 14.98} 16.04 | 
230883.87| 207366.24| 212826.75| 206063.49 
13.88 


25% ACWT/75%TC |] ss 0.92 0.94] 0.93] 0.984 | 
75% ACWT / 25% TC 
50% ACWT/50%TC{ 0.95] 89] BO] | 













RATE = 400 

 _ TRAD] NB|_NB-MOD NB-NPV UICP 
ACWT | 228] 1490] 14.85] 15.92] 13.03 | 
TOTAL COST 231114.77| 207692.54|  213235.53|  206532.50| 217199.76 
AVG YRS _RL 





MADM 
25% ACWT/75%TC | 0.92] 0.95 | 
| 50% ACWT/50%TC | 0.95] 91 | 





RATE = 800 {Default setting for DDA} 


TRAD __Net__NB Mop ene UICP 
12.28 13.03 
38 


25% ACWT/75%TC | 0.92] 0.95] 0.94] 95.95 
75% ACWT/25%TC |]  —s.97] B87] OBE]. 
50% ACWT/50%TC | 0.957 9 Z| 0.95 | 















RATE = 1200 


Po TRAD] NB] ONB-MOD/—ONB-NPV| CUICP 
ACWT | 2B] 1524] 3.89) 14.91] 13.03 | 
AVGYRS RET 13.88] 544] 93f 503] 8.00 
MADM 

L25% ACWT/75%TC | 0.92} 95 9S 9G 0.9 
|75% ACWT/25%TC | 0.97] 85] | BZ | 0.94 | 
50% ACWT/50%TC | 0.95} 90] 93] 9 


Model Ranking by MADM Results 


2. Ee eee aa eee 
Rate = 200 


25% ACWT / 75% TC UICP | NB | NB-NPV NB-MOD TRAD 
75% ACWT /25% TC TRAD UICP NB-MOD | NB | _ NB-NPV 
50% ACWT /50% TC TRAD* UICP* NB-MOD | NB | NB-NPV 

















Rate = 400 
Ba ACW /75% TC NB-NPV NB-MOD TRAD 


eee) UICR 2 | 
75% ACWT / 25% TC TRAD UICP* | NB | NB-MOD NB-NPV 
50% ACWT /50% TC TRAD* |  UICP* | NB | NB-MOD NB-NPV 






Rate = 800 {Default setting for DDA 

25% ACWT/75% TC | NB-NPV* | NB* | UICP*, =| NB-MOD | TRAD 
| 75% ACWT/25%TC [| TRAD | UICP {| NB-MOD | NB | NB-NPV_ | 
50% ACWT/50%TC {| TRAD* | UICP* | NB-MOD | NB | NB-NPV_ | 









Rate = 1200 


25% ACWT/75%1C | NB-NPV* UICP* NB-MOD*’ | NB. |. TRAD 
75% ACWT/25%1C |__ TRAD UICP NB-MOD NB-NPV. |. NB. 
50% ACWT/50%1C | TRAD UICP* NB-MOD NB-NPV | NB 


Note: * indicates models have the same rank and are both ranked as 1. 
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Sensativity Analysis: HIGH DEMAND/ CONCAVE/ OBSOLESENCE RATE 










RATE = .06 
TRAD NB] ONB-MOD[ —ONB-NPV]_CUICP | 
25% ACWT/75%TC | 0.94] 0.92] 92] 0.90} 0.911 | 






75% ACWT / 25% TC 
SO% ACWT/50%TC | 096] BS] BET TTB | 





RATE = .09 












po TRAD] NB] ONB-MOD| —ONB-NPV[ CUICP 
| 25% ACWT/75%TC | 0.93] 94] 93] 0.93] 0.93 | 

75% ACWT / 25% TC 
| 50% ACWT/50%TC {| 0.95] 9B] BE] 0.90 | 








RATE =.12 (Default setting for DDA} 


po TRAD NB ONB-MOD] —ONB-NPV] UIC | 

25% ACWT/75%TC | 0.92] 0.95] 0.94] 95] 0.95 
P75% ACWT/25%TC}] 0.97} 87] OBE 8 
| 50% ACWT/50% TC} 0.95] 91 92 9.95 


















RATE =.15 


| TRAD] NB] ONB-MOD| —ONB-NPV| CUICP 


25% ACWT/75%TC | 0.91} 96] 94] 9G] 0.93 


S 
75% ACWT/25%TC | 0.97] 8990] B87] 0.983 
50% ACWT/50%TC | 0.94] 92] 92] 92] 93 





Model Ranking by MADM Results 


Rate = 0.06 


25% ACWT / 75% TC TRAD NB-MOD | NB. |. UICP |__NB-NPV 
75% ACWT / 25% TC TRAD NB-MOD {| NB | UICP | NB-NPV 
50% ACWT /50% TC TRAD NB-MOD_|._NB. | _UICP NB-NPV 


Rate = 0.09 
| 25% ACWT/75%TC | NB | —OUICP, | NB-NPV_ | TRAD | ~NB-MOD | 
| 75% ACWT/25%TC | TRAD |  UICP | NB | NB-MOD | NB-NPV | 
| 50% ACWT/50%TC | TRAD | UICP | NB | NB-MOD | NB-NPV | 










: i@) 
25% ACWT /75% TC | NB-NPV* | NB* | UICP* NB-MOD TRAD 
75% ACWT / 25% TC TRAD UICP NB-MOD | NB | NB-NPV 
50% ACWT / 50% TC TRAD* UICP* NB-MOD | NB | NB-NPV 





Rate = 0.15 


25% ACWT/75% TC | NB-NPV*’ | NB* | NB-MOD UICP TRAD 
75% ACWT/25% TC | TRAD UICP NB-MOD | NB | NB-NPV 
50% ACWT/50% TC | TRAD UICP NB-MOD |. _NB | NB-NPV 


Note: * indicates models have the same rank and are both ranked as 1. 
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Sensativity Analysis: HIGH DEMAND/ CONCAVE/ SALVAGE RATE 


RATE =.01 


po TRAD NB{ONB-MOD{ONB-NPV[ OO CUICP 
ACWT | 222] 145] 14.03] 14.98] 13.03 | 
AVG YRSRE | 14GB] 4BT G92] 5.04) 8.00 


25% ACWT/75%TC | 0.92 0.95] 94 9S 9S 
TS% ACWT/25%TC | 0.971 87{ 90F BE] 8G 
50% ACWT/50%TC |  —s_—.94] 92 T0994 | 












= 























RATE = .02 {Default setting for DDA} 
TRAD| —sSO*NB | NB-MOD NB-NPV UICP 


aa 
ee TB 72] 408] 15.05] 1.03 | 
MADM 

25% ACWT/75%TC | 92] 95] 49S 0.95 | 
He A ere ee | 0 RL 
| 50% ACWT/50%TC | 0.95] 0.91 0.92] 0.91] 0.95 





RATE =.05 


TRAD| —_sCNB- NB-MOD NB-NPV UICP 
ACWT 


TOTAL COST 229654.07| 210586.55| 215664.68|  209210.58| _219896.04 
AVG YRS RL ease 523] e672) abt 8.00) 


MADM 


25% ACWT/75%TC | 0.93 96] 95] 96] 0. 
75% ACWT/25%TC {| 0.98] 88] 90] 87] 0.97 | 
S0% ACWT/S0%TC | 0.96] 92] 92] 9 





RATE =.15 


po TRAD] NB] ONB-MOD]—OUNB-NPV[_CUICP 
AVGYRSERL | AT 73] 22) 425] 8.00 
-—— as ee 
|25% ACWT/75%TC | 0.96] 96] 9G} 487 
| 75% ACWT/25%TC | 0.98] 8] 2] BI 0.99 | 

50% ACWT/50%TC | 0.97] 0.92] 94} BF} 0.98 | 


Model Ranking by MADM Results 


Rate = 0.01 


25% ACWT/75% TC | NB-NPV* | NB* | UICP* NB-MOD TRAD 
75% ACWT /25% TC TRAD UICP NB-MOD | NB | NB-NPV 
50% ACWT /50% TC UICP* TRAD* NB-MOD | NB | NB-NPV 














— 





Rate = 0.02 {Default setting for DDA 


25% ACWT/75%1C | NB-NPV’ | NB* | UICP® NB-MOD TRAD 


75% ACWT / 25% TC TRAD UICP NB-MOD | NB | NB-NPV 
50% ACWT /50% TC TRAD* UICP* NB-MOD | NB | NB-NPV 





Rate = 0.05 

| 25% ACWT/75% TC |  UICP* | =NB* | ~NB-NPV* | NB-MOD | TRAD 
| 75% ACWT/25%TC | TRAD | UICP {| NB-MOD | NB | NB-NPV | 
| 50% ACWT/50% TC | UICP* | TRAD* | NB-MOD | NB | NB-NPV 


Rate = 0.15 

|25% ACWT/75%TC {| UICP_ | TRAD | NB-MOD | NB | NB-NPV_ | 
| 75% ACWT/25%TC | UICP =| TRAD | NB-MOD | NB | NB-NPV_ | 
| 50% ACWT/50%TC | UICP | TRAD | NB-MOD | NB | NB-NPV 


Note: ° indicates models have the same rank and are both ranked as 1. 


Sensativity Analysis: LOW DEMAND /CONVEX /STORAGE RATE 


RATE =.01 {Default setting for DDA) 


oT TRADT NB] ONB-MOD] —ONB-NPV[ UCP 
AVG YRSRL | 13.88] BOT 54] 5 39f 8.00) 
75% ACWT/25%TC | 0.99] 9396 0.93 | 

| 0.95 | 


50% ACWT/50%TC | 0.98] 095] 97] 0.984 | 















RATE = .03 











po TRAD] NB] NB-MOD] ONB-NPV]UICP 
ACWT 8.36] 1565) 14.16] 16.27] 14.63 | 
AVG YRSRE | BT 5B2T 2] 499] 8.00 
p25% ACWT/75%TC | 0.99] 96 98] 9S 97 
75% ACWT/25%TC | 1.00{ 89] 96] B87] (0.93 | 
50% ACWT/50%TC]  —_0.99f 93] 97 TT 0.95 








RATE =.05 


po TRAD NB] ONB-MOD] —ONB-NPV| UICP | 
AVGYRSRL | 787] 493] 55B] 4 B5Y 8.00 | 
}25% ACWT/75%TC | 99] 96] 98} 95.99 | 
75% ACWT/25%TC | 1.00] 89] 95] BE 1.00 
| 50% ACWT/50%TC | 99] 92] 9G 0.99 





RATE = .07 












po TRADT NB] NB-MOD} ONB-NPV| UICP | 
ACWT CT 24] 19.05] 1853] 19.62] 16.46 | 
0.97] 0.96 | 0.99 

75% ACWT/25%TC | 0.96] 90 OT BB] 1.00 | 
50% ACWT/50%TC | 0.97] 93] 94] 92} 99 


Model Ranking by MADM Results 


Rate = 0.01 {Default setting for DDA 
25% ACWT / 75% TC NB-MOD TRAD | NB |  UICP | NB-NPV 
75% ACWT / 25% TC TRAD NB-MOD UICP | NB | NB-NPV 











50% ACWT / 50% TC TRAD NB-MOD uicP._|..NB|__NB-NPV 


Rate = 0.03 

| 25% ACWT/75%TC | TRAD | NB-MOD | UICP | NB | NB-NPV_ | 
75% ACWT/25%TC | TRAD | NB-MOD | UICP | NB | NB-NPV | 
| 50% ACWT/50% TC | TRAD | NB-MOD | UVICP {| NB {| NB-NPV_ | 


Rate = 0.05 
25% ACWT/75% TC | TRAD* |  UICP* NB-MOD | NB |} NB-NPV 
75% ACWT /25% TC UICP* TRAD* NB-MOD 
50% ACWT /50% TC UICP* TRAD* NB-MOD 
NB 


NB 

| eee 

Rate = 0.07 

[25% ACWT/75%TC |  UICP [| TRAD | NB-MOD | NB | NB-NPV 
75% ACWT/25%TC | UICP | TRAD | NB-MOD | NB 

| 50% ACWT/50% TC} UICP | TRAD | NB-MOD | NB 


Note: * indicates models have the same rank and are both ranked as 1. 
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Sensativity Analysis: LOW DEMAND /CONVEX /ORDER COST 


RATE = 200 


po TRAD NB] NB-MOD] —ONB-NPV/ UCP | 
AVGYRSRAE OT 13-88] 5B6] TT 5.03] 8.00 | 


25% ACWT / 75% TC 
TS%ACWT/25%TC {| 9992 9390 0.93 | 
SO%ACWT/50%TC}] _098f Ss 0.94f 95] 949 | 





RATE = 400 


Po TRAD] NB NB-MOD] ONB-NPV] CUICP | 
VG YRS RL po 13.88P 550] 4] 514] 8.00 | 
25% ACWT / 75% TC 

P7S%ACWT/25%TC {| 0.99] 92] 96] T0983 

[50% ACWT/S50%TC | 0.98f 947] 949 








: 





RATE = 800 {Default setting for DDA} 


po TRAD] NB] OUNB-MOD] —ONB-NPV] OO UICP 
AVGYRSRE | 13.88} BO} SAY 5398.00 
| 75% ACWT/25%TC | 0.99] 0.93] 96] 91 0.93 | 
| 50% ACWT/50%TC{ 0.98} 95] 7 4S 





RATE = 1200 


po TRAD] NB] NB-MOD] —ONB-NPV] CUICP | 
CWT CCS] 13-75{ 13.18] 14.00 13.79 
AVGYRSRE | 13.88f 0377] 5B] 8.00 | 


MADM 
25% ACWT / 75% TC 
ee 0192) 


75% ACWT/25%TC {| 0.99] 0 93{ 0.986 | 
50% ACWT/50%TC {| 0.98} 95] 7] 9 


Model Ranking by MADM Results 


Rate = 200 


25% ACWT/75% TC | TRAD* | NB-MOD* | NB* | UICP* NB-NPV* 
75% ACWT/25% TC | TRAD | NB-MOD UICP | NB | NB-NPV 
50% ACWT / 50% TC TRAD NB-MOD UICP | NB | NB-NPV 


| 0.95 | 


Rate = 400 


25% ACWT/75% TC | NB-MOD TRAD.|._NB. | UIcP NB-NPV 
75% ACWT/25% TC | TRAD NB-MOD uicP. | NB | NB-NPV 
50% ACWT/50%1C | TRAD NB-MOD uicP. | NB | NB-NPV 


Rate = 800 {Default setting for DDA 


25% ACWT /75% TC | NB-MOD* TRAD‘ | NB | UICP NB-NPV 
75% ACWT / 25% TC TRAD NB-MOD UICP | NB | NB-NPV 





50% ACWT/50%1TC | TRAD | NB-MOD UICP.|.NB. | NB-NPV 


Rate = 1200 


25% ACWT/75% TC | NB-MOD* {| TRAD* | NB* =| NB-NPV UICP 
75% ACWT / 25% TC TRAD NB-MOD | NB | __ UICP NB-NPV 
50% ACWT / 50% TC TRAD NB-MOD | NB | __UICP NB-NPV 


Note: * indicates models have the same rank and are both ranked as 1. 
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Sensativity Analysis: LOW DEMAND /CONVEX /OBSOLESENCE RATE 


RATE = .06 


TRAD| NB] _NB-MOD NB-NPV UICP 


TOTAL COST 32578.42| 31942.98{  32004.91| 3183459|  31757.68 
AVG YRS RL p 1856] 8.00] BS] 713] 8.00 


25% ACWT/75%TC |  —_—o.98| 7] 7] G0. 
75% ACWT/25%TC | 99] 993 BBL 88 | 
50% ACWT/50%TC] 99} 0949S] 2092 | 





RATE = .09 












po TRAD NB] NB-MOD| —ONB-NPV| UIC 
ACWT GAT 29] 10.48] 1.64] 10.84 | 
AVG YRS RL | SBS] IT SIT GB] 8.00 | 
|25% ACWT/75%TC | 0.98] 96] 97] 9G] C97 | 
75% ACWT/25%TC | 0.99 BO] 94] 87] 0.92 | 

50% ACWT/50%TC |]  —s.99f 93] 96} TY 0.984 












RATE’= .12 {Default setting for DDA} 


po TRADT NB] ONB-MOD] —ONB-NPV|_CUICP | 
po BOT SAT 539] 8.00 | 

25% ACWT /75% TC 
75% ACWT/25%TC | 0.99 93] 9G] OT 0.93 | 
| 50% ACWT/50%TC] 0.98} 95} 7] A095 





RATE =.15 


po TRAD] NB] ONB-MOD] —ONB-NPV|_CUICP | 
ACWT | 409] 15.96] 15.27] 16.92] 15.05) 


}25% ACWT/75%TC | 0.97] 97] 9B] 9G] 97 | 
| 75% ACWT/25%TC | 0.99} 9 A879 
50% ACWT/50%TC | 0.98} 94] 96] 92] 986 | 





Model Ranking by MADM Results 


Rate = 0.06 


25% ACWT / 75% TC TRAD NB-MOD | NB | UICP NB-NPV 
75% ACWT / 25% TC TRAD NB-MOD | NB | UVICP NB-NPV 
50% ACWT /50% TC TRAD NB-MOD {| NB _ | VICP NB-NPV 


Rate = 0.09 
25% ACWT/75%TC | TRAD | NB-MOD | UICP | NB | NB-NPV_ | 
75% ACWT/25%TC | TRAD | NB-MOD | UICP | NB | NB-NPV | 
50% ACWT/50%TC | TRAD | NB-MOD | UICP | NB | NB-NPV | 


Rate = 0.12 

| 25% ACWT/75% TC | NB-MOD* {| TRAD* | NB | UICP | NB-NPV__ 
| 75% ACWT/25%TC | TRAD | NB-MOD | UICP | NB | NB-NPV 
| 50% ACWT/50%TC [ TRAD | NB-MOD [| UICP | NB | NB-NPV | 


Rate = 0.15 

[25% ACWT/75% TC | NB-MOD |  UICP {| TRAD | NB | NB-NPV | 
75% ACWT/25%TC | TRAD [| UICP | NB-MOD | NB | NB-NPV_ | 

(50% ACWT/50% TC | TRAD | UICP  [ NB-MOD [| NB {| NB-NPV_ | 


Note: * indicates models have the same rank and are both ranked as 1. 
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Sensativity Analysis: LOW DEMAND /CONVEX /SALVAGE RATE 


RATE = .01 

















po TRAD] NB] NB-MOD| —ONB-NPV| CUICP | 
AVGYRSRE | 14GB] BA] SOT 5 AS] 8.00 


MADM 


25% ACWT / 75% TC 
75% ACWT/25%TC | 0.99} 92 95 99 
50% ACWT/50%TC |  —Ss.98f 95] FT 94] | 





RATE = .02 {Default setting for DDA} 


po TRAD NB] ONB-MOD| ——ONB-NPV| CUICP 
AVGYRSRE | 13.88f 5 BO] SA] 539] 8.00 
MADM 

25% ACWT / 75% TC 
P7S%ACWT/25%TC | 99 0.93] 9G} T0983 
[50% ACWT/50%TC | 0.98] 5G 97] 940.9 





RATE = .05 














po TRAD] NB] ONB-MOD] —ONB-NPV] OUICP 
AVGYRSRLE | 12.05] GE] AO] 523] 8.00 


25% ACWT/75%TC | 0.98] 97} 99] 97] 0.97 | 
PSAACWT/25%TC} = 0.99} 27] 0.94 | 
50% ACWT/50%TC]  —ss0.99].95f 98] 94] 0.986 | 


= 
> 
O 
= 





RATE =.15 


po TRAD] CNB] ONB-MOD| —ONB-NPV]OUICP | 


25% ACWT/75%TC | 0.99] 0.99] 00f 96} 0.99 | 
75% ACWT/25%TC | 0.99] 96 00] 0.90 | 0.98 
SO% ACWT/50%TC | 0.99} 97] 00 0.93 





Model Ranking by MADM Results 


Rate = 0.01 

| 25% ACWT/75% TC | NB-MOD* | NB* | ~TRAD* | UICP* | —NB-NPV__| 
| 75% ACWT/25%TC | TRAD | NB-MOD |  UICP | NB | NB-NPV 

| 50% ACWT/50%TC | TRAD | NB-MOD | UICP | NB | NB-NPV | 


Rate = 0.02 


25% ACWT / 75% TC | NB-MOD* | NB | UICP NB-NPV 
75% ACWT / 25% TC TRAD NB-MOD UICP | NB | NB-NPV 
50% ACWT/50%TC | TRAD | NB-MOD UICP | NB | NB-NPV 


Rate = 0.05 
| 25% ACWT/75% TC | NB-MOD | TRAD | UICP | NB { NB-NPV_ | 
75% ACWT/25%TC | TRAD | NB-MOD |  UICP | NB | NB-NPV | 


Rate = 0.15 


a NB 

25% ACWT/75% TC | NB-MOD | TRAD | __ UICP | NB | NB-NPV 
75% ACWT /25% TC | NB-MOD TRAD | JUICP | NB | NB-NPV 
50% ACWT/50% TC | NB-MOD TRAD UICP | NB | NB-NPV 


Note: * indicates models have the same rank and are both ranked as 1. 
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Sensativity Analysis: LOW DEMAND/ CONCAVE/ STORAGE RATE 








RATE = .01 {Default setting for DDA} 
res | TRAD| —NB|__—NB-MOD NB-NPV{_UICP 
ACWT TS] 0] 54] 3B} C7 | 





AVGYRSRE | 138. 8BY29T 720) 5.83] 8.00 
L25% ACWT/75%TC | 094) 94} 96} 9B] 0.94 
75% ACWT/25%TC | 0.98] 82] B89] BOT 0.85 | 
50% ACWT/S50%TC |] 0.96] BB] 2} BET 0.90 







RATE = .03 


el TRAD] —NB{ _NB-MOD NB-NPV UICP 
ACWT OTT] G27] 769] 
TOTAL COST 24232.85 24205.94 

po BITS 78t BBY 540] 8.00 


25% ACWT/75%TC | 0.97] 95] .98| 4G 
75% ACWT/25%TC | 0.99] 8G] 9S BZ] 0.93 | 
50% ACWT/50% TC] 0.98 90 GBB]. | 





RATE =.05 


gee. TRAD} NB! NB-MOD NB-NPV| _UICP | 
ACWT | 5G] 8G] 740] 4G] SO. 
TOTAL COST 25761.12 

87{ 5.35] 15] 504] 8.00 | 


25% ACWT/75%TC | 0.98 95} 9G] 940.98 | 
75% ACWT/25%TC | .99t 85] BB] 0.99 | 
50% ACWT/50%TC {| 0.98] 90] 93] 8B] C99 | 





RATE = .07 


po TRAD] NB] ONB-MOD] —ONB-NPV]UICP | 
CWT 795 871) | eo 
AVGYRSRL | ES] 98] 573] 472] 8.00 
|} 25% ACWT/75%TC | 0.96] 0.96] 9G] 0.9 0.97 
| 75% ACWT/25%TC |  —.92| 87] BG] BE] 0.99 | 
| 50% ACWT/50% TC} 0.94} TF 93]. | 0.98 


Model Ranking by MADM Results 


Rate = 0.01 Default setting for DDA 
25% ACWT / 75% TC NB-MOD TRAD | UICP | NB |. NB-NPV 
75% ACWT/25%1C |__ TRAD NB-MOD |. UICP.|._NB.|__NB-NPV 








50% ACWT /50% TC TRAD NB-MOD |  UICP | NB ___|_NB-NPV 


Rate = 0.03 

| 25% ACWT/75%TC | NB-MOD | TRAD | UICP | = NB | NB-NPV 
75% ACWT/25%TC | TRAD | NB-MOD |  UICP | NB | NB-NPV __ 
| 50% ACWT/50% TC | TRAD | NB-MOD | UicP [| NB | NB-NPV__ 
Rate = 0.05 

| 25% ACWT/75%TC | UICP* | TRAD* | NB-MOD | NB | NB-NPV | 
| 75% ACWT/25%TC | UICP* | TRAD* | NB-MOD [| NB | NB-NPV_ | 
| 50% ACWT/50%TC | UICP_ | TRAD {| NB-MOD {| NB | NB-NPV | 
Rate = 0.15 

| TRAD | NB-MOD {| NB | NB-NPV | 


| UICP - 
75% ACWT/25% TC | UICP | TRAD NB-MOD | NB | NB-NPV 
50% ACWT/50% TC | UICP | TRAD NB-MOD | NB _ | 


Note: * indicates models have the same rank and are both ranked as 1. 
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Sensativity Analysis: LOW DEMAND /CONCAVE /ORDER COST 


RATE = 200 


po TRAD] NB] OUNB-MOD] —ONB-NPV] UIC | 
ACWT S| 7] 57] BO] i786 | 


25% ACWT/75%TC | 094} 0.92] 95] 092] 0.94 | 
75% ACWT / 25% TC 
SO% ACWT/50%7C | 0.96] BS] TBA] 0.90 | 








= 

















RATE = 400 


po TRAD] NB] UNB-MOD| —ONB-NPV]_OUICP 
AVGYRSRE | 13.88] 5 BBY BO] SSO] 8.00 | 


25% ACWT/75%TC {| 0.94] 93] 6920.94 | 
75% ACWT / 25% TC 0.89 
50% ACWT/50%TC | 0.96] 0.87 | O93] 0.85] 0.90 





RATE = 80) {Default setting for DDA} 


po TRAD] NB ONB-MOD] —ONB-NPV| UIC 
eC] (tC COT SATB] S76 | 
ee 

25% ACWT/75% TC | 0.94] 4G 93] 0.94 
| 50% ACWT/50%TC | 0.96] Bf 92] BE] 0.90 













RATE = 1200 












po TRAD] NB NB-MOD| —ONB-NPV] CUICP | 
ACWT |S] 5735 BBT 21] SE 
AVGYRSRE | 13.88] OT 752] 6.08] 8.00 | 
p25% ACWT/75% TC} 0.95] 95] ET 940.94 | 
|75% ACWT/25% TC] 0.98] BET OT BT 8 
SO% ACWT/S0%TC{ ——.97] 91 9387] 0.90 


Model Ranking by MADM Results 


Rate = 200 


25% ACWT /75% 1TC | NB-MOD TRAD UICP | NB | NB-NPV 
75% ACWT /25% TC TRAD NB-MOD UICP | NB | NB-NPV 
50% ACWT / 50% TC TRAD NB-MOD UICP | NB | NB-NPV 











Rate = 400 

NB-MOD_ | TRAD |  UICP | NB | =NB-NPV__ 
| 75% ACWT/25%TC | TRAD | NB-MOD [| UICP | NB | NB-NPV _ 
| 50% ACWT/50%TC | TRAD | NB-MOD | UICP { NB | NB-NPV_| 


Rate = 800 {Default setting for DDA 


25% ACWT / 75% TC NB-MOD | TRAD | VICP | NB | __NB-NPV 


75% ACWT / 25% TC TRAD NB-MOD UICP | NB | NB-NPV 
50% ACWT /50% TC TRAD NB-MOD UICcP. | NB | NB-NPV 





Rate = 1200 

| 25% ACWT/75% TC | NB-MOD | NB | TRAD | UICP | NB-NPV __ 
| 75% ACWT/25%TC | TRAD | NB-MOD | NB |  UICP | NB-NPV | 
| 50% ACWT/50% TC | TRAD | NB-MOD | NB | UICP | NB-NPV | 


Note: ” indicates models have the same rank and are both ranked as 1. 


105 


Sensativity Analysis: LOW DEMAND /CONCAVE /OBSOLESENCE RATE 













RATE = .06 

TRAD] _NB| _NB-MOD NB-NPV UICP 
TOTAL COST 21881.18| _20559.13| 2078988] 2041857] 2056938 
AVG YRS RL 18.56 


25% ACWT/75%TC | 0.95] 0.93] 0.92] 0.93 | 
75% ACWT/25%TC | 0.98] BT] BOT BOT 0.79 | 
50% ACWT /50% TC 


RATE = .09 


TRAD] NB] ONB-MOD] —ONB-NPV[CUICP 
ACWT CT 462 
AVGYRSRL | 5.89 728f B27] EB] 8.00 
|} 25% ACWT/75% TC] 0.95 OS] 9S S| 0.95 | 
75% ACWT/25%TC | 0.98] 85] 87 BO] 0.87 
50% ACWT/50%TC}] 97] 90 B79 
















RATE =.12 (Default setting for DDA} 

Ps TRAD] ——NB/_—sCNB-MOD NB-NPV|UICP | 
ACWT | SO} 544] 38] 576 | 
TOTAL COST 25046.76 23241.83 23400.07 23180.42 23542.96 


Pp 25% ACWT/75%TC | 0.94] 0.94] 0.96] 093]. 
| 75% ACWT/25%TC {| 0.98} 82 B89] BO] 0.85 | 
| 50% ACWT/50%TC | 0.96] BB] 92{ BETO. 


RATE =.15 















po TRAD] NBYONB-MOD] —ONB-NPV| CUICP 
AVGYRSRE | 12.30] 555] BT 52T{ 8.00 | 
| 75% ACWT/25%TC | 0.98] 8G] 9B] BA] 95 | 

50% ACWT / 50% TC po OTT ST BG] 0.986 | 





Model Ranking by MADM Results 


Rate = 0.06 


25% ACWT /75% TC TRAD NB-NPV | NB |. UICP | NB-MOD 
75% ACWT /25% TC TRAD. | NB | NB-NPV NB-MOD |. UICP | 
50% ACWT / 50% TC TRAD | NB | NB-NPV NB-MOD | UICP | 


Rate = 0.09 
25% ACWT / 75% TC TRAD* | NB* — | NB-MOD* NB-NPV 


| UICP* 
75% ACWT / 25% 1C TRAD NB-MOD UICP | NB | NB-NPV 
50% ACWT /50% TC UICP NB-MOD | NB | NB-NPV 






Rate = 0.12 {Default setting for DDA 


25% ACWT /75% TC NB-MOD TRAD VICP | NB___|_NB-NPV 







75% ACWT / 25% TC TRAD NB-MOD | UICP | NB | NB-NPV 

50% ACWT /50% TC TRAD NB-MOD | JUICP | NB | NB-NPV 
Rate = 0.15 

25% ACWT /75% TC | NB-MOD* UICP* | NB | TRAD NB-NPV 


75% ACWT / 25% TC TRAD UICP NB-MOD | NB |_NB-NPV 
TRAD* | vicP* | NB-MOD {| NB | NB-NPV 


Note: * indicates models have the same rank and are both ranked as 1. 
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Sensativity Analysis: LOW DEMAND /CONCAVE /SALVAGE RATE 


RATE = .01 


TRAD|  —NB. NB-MOD NB-NPV UICP 
p43} 608] A] 29] 576 | 
25140.99|  23192.11| _2331459|  23090.34| 23445.31 


25% ACWT/75%TC | 0.94{ 94] 9G] S| 0.94 | 
75% ACWT / 25% TC 
50% ACWT/S50%TC | 0.96] 8B] AFB] 8D 





RATE = .02 {Default setting for DDA} 


eo TRAD] NB] NB-MOD NB-NPV 


ACWT 1) RD C8 S|) 9 


TOTAL COST 25046.76| _-29241.83| _23400.07| _-23180.42| 2354296 
AVG YRS RL | SEE) ES a 


25% ACWT/75%TC | 0.94 0.94] 96] 93] 0.984 | 
P75% ACWT/25%TC | 0.98} 82] 89] BO] BS 
[50% ACWT/50%TC | ss o.96{ 88 92} BE] 0.90 





RATE =.05 






_ TRAD] ~— NB] ~—ONB-MOD] —ONB-NPV{_CUICP | 
ACWT TS] B] B15 76 | 
AVGYRSRE | 12.05] 4] 706] 5G] 8.00 | 

25% ACWT/75%TC | 0.96] 95] 97 S| 09S 
| 7I5% ACWT/25%TC | ss o.99| 85 94] 84] 89 
50% ACWT/50%TC | 0.97} 90 BD] 09D 












RATE =.15 






po TRAD] NB] ONB-MOD| —ONB-NPV] OUICP 
ACWT | A} BB] SOT 684] 5.7 
}25% ACWT/75%TC | 0.98} 97 OT 4] 0.98 | 
| 75% ACWT/25%TC | 0.98] 90 1.00 B49 
50% ACWT/50%TC] _—s 0.98] 0.93] 00] 89] 0.97 | 


Model Ranking by MADM Results 


Rate = 0.01 


25% ACWT /75% TC | NB-MOD UICP TRAD | NB |. NB-NPV 
75% ACWT / 25% TC TRAD NB-MOD UICP |  =NB | NB-NPV 
50% ACWT/50%TC | TRAD | NB-MOD UICP | NB | NB-NPV 
















Rate = 0.02 {Default setting for DDA 


25% ACWT /75% TC | NB-MOD TRAD UICP | NB | NB-NPV 
75% ACWT/25% TC | TRAD | NB-MOD UICP | NB | NB-NPV 
50% ACWT /50% TC TRAD NB-MOD UICP | NB | NB-NPV. 






Rate = 0.05 

| 25% ACWT/75%TC | NB-MOD | TRAD | UICP | NB-NPV_ | NB | 
| 75% ACWT/25%TC | TRAD | NB-MOD | UICP | NB | NB-NPV_ | 
| 50% ACWT/50%TC | TRAD | NB-MOD | uiCP | NB | NB-NPV_ | 


Rate = 0.15 


25% ACWT /75% TC | NB-MOD UICP | NB | 
75% ACWT/25%TC | NB-MOD | TRAD | UICP | NB | NB-NPV 
50% ACWT /50% TC | NB-MOD TRAD UICP | NB | 


Note: ° indicates models have the same rank and are both ranked as 1. 
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APPENDIX D. SIMULATION CODE 


This appendix contains the following pascal code for the 


Simulation: 
NAME TYPE OF CODE PAGE 
- UICP_Simulator main program 109. 
- toolbox unit 144. 
- unirand unit ae 
- PDUnit unit 5 3. 
- PQueue. unit 6 Se 
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program UICP Simuiator (input, ouput); 


{SM 54000,0,0) {{r+} (SN+,E+) (SG+) 


uses dos, crt, toolbox, unirand, PDUnit, pqueue; 


type quarterArray=array [1..100] of real; 
weeklyArray=array [1..1300] of real; 
repArray=array [1..750]) of real; 
qtrintArray=array [1..100] of integer; 
changeRealArry = array [1..5]} of real; 
changelIntArry = array [1..5] of integer; 
pas2field=string([15]; 


deseriptType=string[40}; 


const COEFF1l=1.386; 
POWER 1=0.746; 
COEFF2=3.869; 
POWER2=1.378; 
MAXPLT=14.0; 
MINPLT=2.0; 
ERROR=1 . OO000000000000E-0010; 
YRSERR=8; 
MINERR=5; 


var wklyObserv:weeklyArray; 
observ, frest, mad, EOOQArry, ROLevelArry,ERRArry, 
SSADDBO, SSADD, SSSMA ,meanDmdArry ,varDmdArry, investQtr,qtrSMA:quarterArray; 
stepIndArry, trndiIndArry,mkCodeArry:qtriIntArray; 
observType,distrType, output Type, seedtype,wkDataType,qtrDataType, 
PDDataType, repStatType, ERRType, anal IndType: char; 
numberRep,i, numberOfReps,numberOfQtrs,numberOfWks,markCode, initInv,simCount: rnteger; 
meanDemand, varDemand: real; 
noInt,trendOn,StepOn,nmbrSteps, nmbrTrends, TWUS, orderCount: integer; 
s,seedIndex,numOtr: integer; 
currSeed: longint; 
Inputfrle,outputfile:text; 
noReal, f1xERR: real; 
stringval:pd82field; 
stop: boolean; 
Startstep, startrnd, endtrnd: changelIntArry; 
stepmult, trendcoeff, trendpower: changeRealArry; 
hourl,minutel,secondl,hdSecl, hour2,minute2,second2 ,hdSec2:word; 
out PileName: string; 
OSHeap, BOHeap: PriorityQueueType; 
ADDBO, ADD, SMA, Invest: real; 


S1mADDBO, simADD, simSMA, simInvest, simOrdercCount:real; 


IGN, 





ciADDBO,c1ADD,ciSMA,cilnvest,ciOrderCount:real; 


ciDisposals,cidisposalCount,ciEndOH,c1iEndOS: real; 


varADDBO,varADD,varSMA, varInvest,varOrdercount: real; 


varDrsposals,varD:sposalCount , varEndOH, varEndOS:1eal; 


disposalCount,disPosals,endOH,endoS: integer; 


simDisposalCount,simDisposals,s1mEndOH,srmEndOS:real; 


runDesecript :descriptType; 


totCost ,holdTC,orderTC,shortTC, salvTR: real; 


totCostArry,holdTCArry, order: TCArry, shortTCAr ry, salvTRArry:quarterA? ray; 


simTotCost,simHoldTC,simOrderTC, simShortTC, srmSalvTR::eal; 


varTotCost ,varHoldTC,varOrderTC, varShortTC, varSalvTR: real; 


c1TotCost,ciHoldaTC, c1OrderTC, «iShortTC, c1SalvTR: real; 


procedure Frontscreen; 


begin 


clrser; 
writeln; 
writeln; 


writeln; 


writeln; 


writeln; 
writeln; 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 


writeln 


Delay (1500); 


classe: 


end; 


. 
P 
" 
. 
‘ 
. 
. 
“ 
. 
. 


{For 


procedure runtype {var 


var 


done: 


1,maxStart:integer; 


demandInFile: 


begin 


boolean; 


var 


var 


var 


Var 


var 


var 


val 


var 


UOOUO TOO ORO OURO tit GOD NOON OO ORG) ic 
, 


: UICP LEVELS FORECASTING aed 
SIMULATOR poner 
: FOR CONSUMABLES Se 


G Ge Gy Robrllara Lr, se 
. D. C. Miller LCDR,SC 


a, 


POOONOMOTOO OOO OOS OCO OOO SD OODG OTN s 
, 


1500 ms} 


distrType,outputType,wkDataType,qtrDataType, 

PDDataType, repStatType, ERRType, anal IndType: char; 
numberOfOtrs,numberOfWks,numberOfReps, seedIndex: integer; 
meanDemand, varDemand: real; 
numYrsOH,numYrsERR: real; 
inputfile,outputfile: text; 
frest,mad: quarterArray; 
seeds: seedArryType; 


outFileName: string; 


runDescript:descriptType) ; 


string; 


a0 


writeln; 

Writeln (" ‘*** THIS SCREEN W1LL ALLOW SELECTION OF RUN TYPE OPTIONS ***'); 
done: =FALSE; 

writeln; 

writeln; writeln; 

write ('Ente1 the number of replications (from ! to 750) to be run: ‘° ); 
numberOfReps:=Get_Intege1(1,750); 

writeln; 

write('Ente: Run Description: '); 

readin (runDescript) ; 

writeln; 

writeln('Quarterly observations will be generated based on your selectron of drstiibution'); 
writeln('(Porsson or Normal) and seed selectron.'); 

writeln; 

repeat 


wirteln ('Random Numbe1 Generator Seed Selection: '); 


writeln; 
writeln (' 1] - Default array - unique seed for each replication'); 
writeln (' 2 - Select seeds - max numbe1 of replications rs 100'); 
writeln; 


fomee {*Choice: '}; 
seedtype:=readkey; 
writeln (seedcdtype) ; 
writeln; 
case seedtype of 
el: begin 
done:=TRUE; 


maxStart :=20001]1-numberOfReps; 


- 

— 
ae 

a 
- 


write(*Enter Random Seed Start Index ( 1 to ',maxStart:2, 
seedIndex:= Get_Integer(l,maxStart); 
end; 
a2. begin 
done: =TRUE; 
1f numberOfReps » 100 then numberOfReps:=100; 
for 1 := 1 to numberOfReps do begin 
write ('Enter Seed value for replication ',i,° : '); 
seeds([i):=Get_LongInt (1,2147483646); 
writeln; 
end; {for} 
end 
end 
until done=TRUE; 
er scr; 
writeln {° 7o** RUN SELECTION OPTIONS CONTINUED ****'); 
writeln; writeln; 
write (‘Enter the number of simulation quarters: ‘' }); 
numberOfOtrs:=Get_Integer(1,100); 
numberOfWks:=13*NumberOfOtrs; 
writeln; 
done: =FALSE; 


repeat 


dS 





writeln (‘Type of Distribution: ")-; 


writeln; 

writeln ({‘ I - Normal‘); 
writeln {' 2 - Poisson'); 
writeln; 


write: f° Chorce: “yi 
distrType: sreadkey; 
writeln (distrType); 
writeln; 


case distrType of 


TY 3S ferzxe(Qia 
done: =TRUE; 
write (‘Enter quarterly mean demand: '‘'). 


meanDemand: =Get_Real (0.0001 ,999999.0); 
writeln; 
write (‘Enter demand variance: ‘); 
varDemand: =Get_Real (0.0001,999999. 0); 
writeln 
end; 

‘2°: begin 


done: =TRUE; 


write (‘Enter quarterly mean demand: ‘'); 


meanDemand: =Get_Real (0.000],999999.0); 
var Demand: =meanDemand; 
wiiteln; 
end 
end 
until done=TRUE; 
frest [1] :=meanDemand; 
mad{1]:=COEFF1*exp(POWER1]*In(frest{[{1])); 
done: =FALSE; 


eleiseia: 


writeln {' **44 RUN SELECTION OPTIONS CONTINUED ****'); 


writeln; 


repeat 


writeln (‘Initial Inventory and Outstanding Reorders 


Selection: 


writeln; 

writeln (' 0O - Default: Initial Inv = EOQ + Safety stock'); 

writeln (' 1 - User specified Initial Inv, No Outstanding Reorders'); 
writeln; 

write ('Choice: '); 


anal IndType: =readkey; 

writeln {analIndType); 

writeln; 

case analIndType of 
*O': done: =TRUE; 


‘1': begin 


write('Enter initial inventory in years of annual demand 


numY2sOH:=Get_Real(0.0,1000.0); 
done: =TRUE; 


end; 


1 A 


a ey 





end; {case} 
unt:] done=TRUE; 
done:=FALSE; 


Cluser : 


done: =FALSE; 
cliser; 
wiiteln (' **** RUN SELECTION OPTIONS CONTINUED °****) 
writeln; 
repeat 
writeln (‘Type of Economic Retention Model Selection: ‘'}; 
writeln; 


wiitteln {*° 0 - No economic retention model used'); 


Wrateln (° 1 - Navy UICP-B20'); 

writeln (' 2 - Net Benefit Model'); 

writeln (' 3 - Modified Net Benefit Mocel'); 
wiiteln (' 4 ~- NPV Net Benefit Model"); 


writeln {°° 5S - Tradition Retention Mocdel'); 
writeln {' 6 - Fixed Retention Requirement (in years)'); 
writeln; 
Winter ( Cholce: ‘); 
ERRType: =reacikey; 
writeln (ERRType); 
writeln; 
case ERRType of 
aoe. Os Aone:=TRUE; 
‘o': begin 
wiite(‘Enter retention requirement in years : '‘'); 
numYrsERR:=Get_Real (0.0,1000.0); 
done: =TRUE; 
end; 
end; {case) 


until] done=TRUE; 


q FALSE; 
c l; 
wi:teln (' oo] RUN SELECTION OPTIONS CONTINUED *****)); 


writeln; writeln; 
repeat 
writeln; 
Mrictein (°Sencdi Output to: '); 
writeln; 
writeln {* 1 - Screen’); 
writeln ({' 2 - File‘); 
writeln; 
write ('Choice: ‘); 
output Type: =readkey; 
writeln (outputType) ; 
case outputType of 
mle: begin 
done: =TRUE; 


assign(output file, ‘*con'); 


ia ae 


end; 


ho 


begin 
done: =TRUE; 
repeat 

writeln; 

write (‘Enter Path and Filename: ‘); 

readin (outFileName) ; 

writeln; 

writeln (‘Path and FileName entered: ‘,outFileName); 

writeln; 

Write ("Is this correct? (Y o10N): 3); 
until Get_Answer; 
assignitoutput file, out FileName) ; 

end; 
end; 
until done=TRUE; 
wkDataType: ='0°'; 
wiiteln; 
write{'Include Weekly SDR Data? {Y or N): ‘); 
if Get_Answer then wkDataType:='1'; 
qtrDataType:='0'; 
writeln; 
write('Include Quarterly SDR Data? (Y or N): '); 
if Get_Answer then qtiDataType:='1'; 
PDDataType:='0'; 
writeln; 
write{'Include Quarterly demand, forecast and PD82/8b6 Data? (Y 01 N): ‘); 
if Get_Answer then PDDataType:='1'; 
repStatType:='0'; 
writeln; 
write('Include Replication Statistics? (Y or N): °'); 
if Get_Answer then repStatType:='1'; 


end; 


procedure RunAgalin (var outputfirle:text;var runDescript:descriptType; 
var outputType, ERRType: chai ; 
var stop: boolean; 
var numYrsERR:1eal; 


var outPFileName:string); 


var demandInFile: string; 
donel: boolean; 
begin 
stop: =FALSE; 
elvsea: 
writeln {° #*4* RE-RUN SIMULATION OPTIONS SCREEN *****); 
writeln; 
writeln('Re-running the srmulation will marntarn the same run-type parameters, but will'); 
writeln('allow the user to change the destination (output) file and vary NIIN'); 


writeln('and model parameters.'); 
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wiiteln; 
write(*Dc you wish te re-run the simulatica)); fY or N 
1f Get_Answer then begin 
writeln; 
wiirte('Change Run Description? (Y 01 N): 
if Get_Answer then begin 
wiiteln; 
write (‘Enter Run Description: '); 
readin (runDescript); 
end; 
writeln; 
done] :=FALSE; 
write('Change Economic Retention Model? (¥Y 01 N): ‘); 
if Get_Answer then begin 
writeln; 
writeln; 
done]: =FALSE; 
writeln; 
repeat 
writeln (‘Type of Economic Retention Model Selection. '); 
writeln; 
writeln (* 0 - No economic retention model used'); 
writeln (*° 1 - Navy UICP-B20'); 
writeln (' 2 - Net Benefit Model’); 
writeln (*' 3 - Modified Net Benefit Model'); 
writeln (' 4 - NPV Net Benefit Model'); 
writeln (° 5 - Tradition Retention Model’); 
writeln (' 6 - Fixed Retention Requirement (in years)'); 
writeln; 
mrate (*Choice: ‘); 
ERRType: =readkey; 
writeln {(ERRType); 
writeln; 
case ERRType of 
Seen 50s done] :=TRUEs 
Bie begin 
write(‘Enter retention requirement in years : °'); 
numYrsERR:=Get_Real (0.0,1000.0); 
done]: =TRUE; 
end; 
end; {case} 
unt:1 donel=TRUE; 
elrscr; 
writeln (’ ##** RUN SELECTION OPTIONS CONTINUED ‘****'); 
writeln; writeln; 
end; {if} 
tf outputType='2' then begin 
writeln; 
write('Change Output File? (Y 01 Nj): ‘'); 
if Get_Answei then begin 


repeat 


dla Wes 





wirteln; 
write (‘Enter Output Path and Filename: '); 


readin (outFileName) ; 


writeln; 

writeln (‘Path and FileName entered: ',outFileName) ; 
writeln; 

write {(7ls this comect 25 (Y tor Nii) ji: 


until Get_Answer; 
assign{output file,outFileName) ; 
end; 
end; 
end else begin 
stop: =TRUE; 
end; 
clisci: 


end; 


function GetMarkCode {t,oldMark:integer; frest, unitPrice:real):integei; 


begin 
if t=] then begin 
if frest < 0.25 then getMarkCode:=0; 
if {frest »-= 0.25) and {(frest <« 2.0) then begin 
if (unrtPrice -= 300.00) then begin 
get Mar kCode: =3; 
end else begin 
getMarkCode:=1; 
end; 
end; 
if frest >= 2.0 then begin 
if (unitPrace‘*frest) »= 6000.0 then begin 
getMarkCode: =4; 
end else begin 
getMarkCode:=2 
end; 
end; 
end else begin 
getMarkCode: =oldMark; 
if oldMark = 0 then begin 
if frest >= 0.5 then begin 
if {unitPrice -= 300.00) then begin 
getMarkCode: =3; 
end else begin 
getMarkCode:=1; 
end; 
end; 
1f frest »=3 then begin 
rf (unitPrice‘frest) »= 600.0 then begin 
getMarkCode: =4; 


end else begin 


ue 





getMarkCode: =2 
end; 
end; 
end; 
if (oldMark=1) 0: (oldMaik=3) then begin 
Peeticst -=3 then begin 
1f (unitPrice*frest) -= 000.0 then begin 
getMarkCode: =4; 
end else begin 
getMarkCode:=2 
end; 
end else if unitPrice <= 200 then begin 
getMarkCode:=1; 
end else 1f unitPrice »= 400 then begin 
getMarkCode: =3; 


end; 


It 
°o 


if frest -= 0.25 then getMarkCode: 
end; 
1f (oldMark=2) o: (oldMark=4) then begin 
af frest <= 1.0 then begin 
if (unitPrice += 300.00) then begin 
getMarkCode: =3; 
end else begin 
getMarkCode:=1; 
end; 
end else if (unitPrice*frest) »= 800.00 then begin 
getMarkCode: =4; 
end else if (unitPrice*frest) <= 400.00 then begin 
getMarkCode: =2; 
end; 
if frest -= 0.25 then getMarkCode:=0; 
end; 
end 


ena; 


procedure InitializeArrays (var observ, EOQArry, ROLevel, SSADDBO, SSADD, 

SSSMA, ERRArry: quarterArray; 

var stepIndArry, trndIndArry,mkCodeArry: qtirintArray; 

numberOfOtrs,numberOfWks,numbe1Rep: integer; 

meanDemand: real; 

var wklyObserv:weeklyArray; 

var meanDmdArry,varDmdArry: quarterArray; 

var totCostArry,holdTCArry,orderTCArry, 
short TCArry,salvTRArry,investQOtr, 
qtrSMA:quarterArray} ; 


var t:integer; 
begin 


for t:=l to numberOfOtrs do begin 


observ ([(t]:=0.0; 


jUAEy, 


meanDmdAi ry[t]:=0.0; 


Vv 


arlmdAaArry[t]:=0.0; 


EOOAR ry (tt ]=0..0; 


ROLevel [t]:=0.0; 


S 
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SADDBO[t] :=0.0; 


SADD[t] :=0.0; 


SSSMA [t] :=0.0; 


stepIndArry([t]) :=0; 


trndindaArry([t] :=0; 


mkCodeArry [t]:=0; 


if numberRenp = 


totCostArry([t]:=0.0; 


holaTeArry [oc] 2—0).0: 


OVNaeiw TCA Tey [ital =O). 0); 


shortTCArry [t]:=0.0; 


salvTRArry [Et ]/2=0.0; 
ERRArry (t]:=0.0; 
investOtr[t]:=0; 


qtrSMA(t]<:=0; 


end; ate} 


end; 


for 


1 then begin 


t:=l1 to (numberOfWks) do begin 


wklyObserv([(t]:=0.0; 


end; 


end; 


procedure LoadObserv (var observ, frest,mad:quarterArray; 


var wklyObserv:weeklyArray; 


var meanDmdArry,varDmdArry:quarterArray; 


observType,distrType:char; 


numberOfQOtrs,numberOfWks, repNum,simCount: integer; 


var trendiInd,stepiInd,nmbrSteps, nmbrTrends: integer; 


meanDemand, varDemand: real; 


Val 


var 


var 


var 


Varioonehare 


1, t, Min, observWeek;: 
randnorm, currMeanDmd, 


demandInFile:string; 


begin 


input file:text; 
seeds: seedArryType; 
startstep, startrnd, 


stepmult, trendcoeff, 


integer; 


endtrnd: changelntArry; 


trendpower: changeRealArry) ; 


initTrendMean, coeffVar, qtrObserv:real; 


if (repNum = 1) and (simCount = 1) then begin 


for i:=l1 to 5S do begin 


startstep[1] :=0; 
stepmult[1]:=0.0; 


end; 


startrnd[{i]:=0; endtrnd[i] :=0; 


trendcoett(7 ):20-0- 


trendpower[i]:=0.0; 
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nmb:i Steps: =0; 
nmbrT1 ends: =0; 
end; {1f) 
currMeanDmd: =meanlemand; 
coef fVal :=sqrt (varPemand) /meanDemand; 
fo. t:=0 to inumbe1Of0Ot:s) do begin 
1f {t=0) and (1epNum = 1!) and (simCount=1} then begin 
Gos=*Y'; 
wiiteln; 
wiite(’Do you wish to vary mean demand rate ovel time? {Y o1 N): ‘' 
1f Get_Answe1 then begin 
Some one > 
stepInd:=0; 
tirendind: =0; 
elirscr; 
writeln; 
wiiteln (' **#* Mean Demand Variants *** ins 
writeln; 
writeln (‘You have the option to vary mean demand rate ovei time. If the noimal'); 
writeln (‘distiibution was selected, variance will also change to maintain your'); 
writeln (‘original variance to mean ratio. You may choose between step change‘); 
writeln (‘or trend or any combination of the events. If more than one event is'); 
writeln (‘chosen to occur at the same time, step changes will occur first.'); 
writeln (‘A maximum of 5 occurances of each event is allowed.'); 
writeln; 
So.- Y'; 
wiite {'Do you still wish to vary mean demand rate ove time? {Y or N): ‘'); 
if Get_Answer then begin 
so. = Ni"; 
elxsci; 
writeln(' #** Step Changes Screen ***'); 
writeln; 
write ('Do you wish to have step increases o1 decreases? (Y or N): '); 
if Get_Answer then stepInd:=1]; 
if stepInd=1] then begin 
writeln; 
write(’Ente: the number of steps changes desired (max 5): ‘'); 
nmbiSteps:=Get_Integer (1,5); 
writeln; 
writeln(‘The step function 1s of the form: Mean(t}) = A * Mean(t-1).'); 
writeln(‘You must specify the value of “A” for each step.'); 
min: =l; 
for 1:=] to nmbrSteps do begin 
writeln; 
Mice ime oten °,1,°%"): 
writeln; 
write (‘Step quarter: ‘); 
startstep(i]:=Get_Integer (min,numberOfOtrs) ; 
writeln; 
write {*Step Multiplier {A): ‘); 


stepmult [1] :=Get_Real (0.00001 ,9999.0); 


IA, 





writeln; 
Minz=staitstep| i | 
end; 
end; 
SCT 
writeln(' Trend Setting sseLeenmaunanns 
writeln; 
write (‘Do you wish to have trends? (Y or N):'); 
1f Get_Answer then trendInd:=1]; 
if trendiInd=1 then begin 
writeln; 
writet'Ente1 the number of trend periods desired imax 5): '): 
mmbitrends: =Get_Integer(1,5); 
writeln;: 
writeln('The trend function 1s of the form:'); 
wiriteln(' Mean(t) = InitTiendMean * ( 1 +A * £(0) ** B)*) 
writeln('where t(0) is reset to "l* at the beginning of each trend pertod'); 
writeln('and InitTrendMean 1s the Mean at the beginning of the trend period.'); 
writeln('Parameters A and B must be specified for each trend period.'); 
min:=1; 
for 1:=l to nmbrtrends do begin 
writeln; 
peer ig (Mp esytel F528 8 Ol) G 
writeln; 
WEIte [Start quarter: sn, 
startrnd[1):=Get_Integer({min,numberOfOtrs) ; 
writeln; 
write (‘End quarter: ‘}; 
endtrnd[i):=Get_Integer(startrnd(1],numberOfOtrs) ; 
writeln; 
write (‘Trend coefficent (A): ‘'); 
trendcoeff[i]:=Get_Real (-9999.0,9999.0); 
writeln; 
write ('Trend power (B): ‘); 
trendpower[1):=Get_Real (-9999.0,9999.0); 
writeln; 


min:=endtrnd[i]¢l; 


end; {for} 
end; {if trend=1) 
end; {1f getans) 
end; {if getans) 


end else 1f t + 0 then begin 
if SS='Y" then begin 
meanDmdArry [t ] : =meanDemand; 
1f distrType='1l' then begin 
varDmdArry [t]) :=varDemand; 
end else begin 
varDmcdArry [t ] :=currMeanDmd; 
end; 
end else begin 


if steplnd = 1 then begin 


ZG 


for 1:=l1 to nmbrSteps do begin 
if t = startstep[i] then currMeanDmd:=stepmult (1) ‘cur:MeanDmd; 
end; 
end; 
if trendiInd = 1 then begin 
for 1:=] to nmbrTrends do begin 
if t = startrnd[i)] then initTrendMean: =currMeanlma: 
Poa Scartrnndii}) and (t «= endtrnd(1]'} then begin 
currMeanDmd: =1nitTrendMean’ (l+trendcoeff[1]* 
{exp {trendpower [1] *In({t-startrnd[1j+1)))); 
if currMeanDmd « 0.0 then cur:rMeanDmd:=0.0; 
end; 
end; 
end; 
meanDmdArry [t] : =cur1 MeanDmd; 
if distiType='!1' then begin 
varDmdArry [(t ]:=sqr {coef fVar*currMeanDma) ; 
end else begin 
varDmdArry [t ] :=currMeanDmd; 
end; 
end; 
1f distrType='1' then begin 
randnorm: =GetNormal; 
qt rObserv: =round(meanDmdArry [t}+{randnorm‘sqrt (varDmcArry [t]))); 
if qtrObserv < 0.0 then qtrObserv:=0.0; 
for 1:=1 to round{qtrObserv) do begin 
observWeek: =Get Uni formlint (13); 
wklyObserv[(t-1)*13+observWeek] := 
wk lyObserv [(t-1) *l13+observWeek] +1; 
end; 
end rit) 
else if distrType='2' then begin 
qt rObserv: =Get Poisson (meanDmdArry [t]); 
for 1:=l1 to round{qtrObserv) do begin 
observWeek: =Get Uni formInt (13); 
wklyObserv { (t-1) *13+observWeek]:= 
wklyObserv [ (t-1) *13+observWeek } +1; 
end; 
end; {else} 
observ([t]:=qt rObserv; 
end; {else,1rf] 
end; {for] 


elrscr: 







end; 


procedure Forecast {var observ, frest, mad:quarterArray; 
var steplIndArry, trndIndArry,mkCodeArry: qtrintArray; 


numberOfQtrs,repNum: integer; unitPrice:real); 


const ALPHA=0.1!; 


a2 0 


STEPBOUND 1 =3.0:; 
STEPBOUND2=2. 0; 


var upper, lower, sum, sampleMean, sampleStdDev, stdDevToMean: real; 
upInd, downInd, steplnd, trendind, trenduUp, 
trendDn, t, 1, j, W, S, table:integer: 
kendTest, lowDemand: boolean; 
begin 
writeln('Running Replication # ',1epNum) ; 
mkCodeArry[1]):=sgetMarkCode (1,0, frest{1]),unitP11ce); 


Usd Neale =O AGownlncts =r 

for t:=2 to numberOfOtis do begin 

lowDemand: =FALSE; 

trendind: =0; 

stepInd: =0; 

1f ((mkCodeArry[t-1] = 0) or 

if lowDemand then begin 
upper: =STEPBOUND] 'frest [t-1]; 
lower: =0.0; 


end else begin 


{mkCodeArry[t-1] = 


{Compute quarterly forecast) 


1} or {mkCodeArry([t-1] 


upper:=frest [t-1])+1.25*mad[t-1] *STEPBOUND2; 


lower:=frest [t-1]-1.25*mad[t-1] *STEPBOUND2;: 


end; 


if (lowDemand and (observ[t-1) < 5)) or 


({observ[t-1]) - upper) and 
uplI nei: =0; 


downlnd: =0; 


(observ[t-1] 


>= 


lower)) then begin 


frest[t]:=ALPHA*observ[t-1}]+{1-ALPHA) *frest{t-1]; 


mad([t]:=ALPHA* (abs (observ [t-1]-frest [t-1])})+{1-ALPHA) ‘mad[t-1]; 


end else begin 


1f ({observ[t-1] > upper) 
((observ[t-1] « 


if t>4 then begin 


frest [t]:=(observ[t-4]+observ[t-3] +observ[t-2]+observ[t-1]})/4; 


end else if t = 4 then begin 


and (upInd=])) or 


lower) and (downInd=1)) then begin 


frest[t]:=(observ[t-3)+observ[t-2]+observ[t-1])/3; 


end else 1f t = 3 then begin 


frest[t]):=(observ[t-2]+observ(t-1])/2; 


end; 


rf frest(t) = 0.0 then mad[(t):=0.0 


else mad[t):=COEFF]*exp(POWER1*In(frest(t))); 


steplnd:=1; 
uplna: =0; 
down 1nd: =0; 
end else begin 
if ((observ[t-1] » upper) and 
uplnd:=1; 
Prest [eG] = trese [t—) | 
mad[t}:=mad[t-1]; 


end else begin 


(upInd=0)} then begin 
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=3)) then lowDemand: =TRUE; 


1f ((observ(t-1)] lower) and {(downInd=0); then begin 
Gout nes =. ; 
frest(t):=frest(t-1]: 
mad(t]:=mad[(t-1]); 
enci; 
end; 
end; 
end; 
1f (t-4) and (stepInd=0) then begin {Conduct Kendall Trend Test } 
sum: =0.0; 
if t -= 8 then begin 
forsi'=) to t-1 clo begin 
sum: =sum+observ [1]; 
end; 
sampleMean: =sum/ (t-1); 
sum: =0.0; 
for 1:=] to t-1 do begin 
sum: =sum+sqr (observ ([1]-sampleMean); 
end; 
sampleStdDev:=sqrt (sum/ (t-2)); 
end else begin 
for 1:=t-8 to t-1 do begin 
sum: =sum+obseirv([i]; 
end: 
sampleMean: =sum/8; 
sum: =0.0; 
for 1:=t-# to t-1 do begin 
sum: =sum+sqr (observ [i] -sampleMean) ; 
end; 
sampleStdDev:=sqrt (sum/7); 
end; 
if sampleMean » 0.0 then begin 
stdDevToMean: =sampleStdDev/sampleMean 
end else begin 
stdDevToMean: =939999.0 
end; 
kendTest :=false; 
if (sampleMean -= 3.0) and (stdDevToMean <= 1.75) then begin 
kendTest: =true; 
if stdDevToMean » 1.0 then begin 
table:=3; 
end else begin 
table:=2; 
end; 
end; 
if ({SampleMean ~= 1.0) and (sampleMean <« 3.0)) and 
(stdDevToMean «= 1.75) then begin 
kendTest :=true; 
if stdDevToMean »~ 1.25 then begin 
table: =3; 


end else begin 


LZ 


table:=2; 
end; 
end; 
if ((sampleMean -= 0.125) and (samp]eMean enc and 
(stdDevToMean <= 2.00) then begin 


kendTest :=true; 


table:=2; 

end; 

1f kendTest=true then begin {Conduct Kendall S-Test for Trend) 
W: =8; 


if (sampleMean -= 3.0) and {(sampleMean <« 9.0) then begin 
rf (stdDevToMean - 0.30) then W:=5; 
end; 
if (sampleMean ~= 9.0) and (sampleMean ~« 20.0) then begin 
1f (stdDevToMean - 0.93) then W:=6; 
if (stdDevToMean < 0.28) then W:=4; 
end; 
if {sampleMean -= 20.0) then begin 
if (stdDevToMean < 0.53) then W:=6; 
if (stdDevToMean < 0.28) then W:=4; 
end; 
1€ W > (t-1) then W:=( (t=) )} div 2) 2: 
S:=0; 
for 1:=(t-W) to (t-2) do begin {Compute Kendall S-Statistic} 
for j:=(1+1) to (t-1} do begin 
if obseiv{1i} < obseiv(j]) then S:=S+l1; 
1f obseiv{i] » observ{j] then S:=S-1; 
end; 
end; {for} 
1f table = 2 then begin 
if W = 4 then begin 
trendUp:=4; trendDn:=-4; 
end; 
if W = 6 then begin 
trendUp:=9; trendDn:=-9; 
end; 
if W = 8 then begin 
tiendUp:=13; trendDn:=-13; 
end; 
end else begin 
if W = 4 then begin 
trendUp:=6; trendDn:=-6; 
end; 
if W = © then begin 
trendUp:=11; trendDn:=-11; 
end; 
if W= 8 then begin 
trendUp:=lo; trendDbn:=-lo; 
end; 
end; {if} 


trendiInd:=0; 
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rf S += trendUp then trendind:=1; 
if S -= trendDn then trendind:=1; 
rf trendInd =1 then begin 
sum:=0.0; 
for r:={t-4) to (t-1) dc begrn 
sum: =sum+observ(1]; 
end; 
Prest (t)z=sum/4; 
if frest{t} = 0.0 then mad[t}:=0.0 
else mad(t]}:=COEFF] texp(POWERI]*In(frest{t}) i; 
end: {if} 
end; {1f] 
end; {:f) 
mkCodeArry([t}:=getMarkCode (t,mkCodeArry([t-1},frest[t],unitPrice)-: 
steplndArry(t]:=stepiInd; 
trndiIndAaArry [(t ] :=trendind; 
end; {for} 


end; 


procedure LoaclLevels (var frest, mad, observ, EOQArry, ROLevelArry, 
SSADDBO, SSADD, SSSMA:quarterArray; 
var mkCodeArray:qtriIntArray; 
var numberOfOtrs: integer; 
prbBrkPt:rnteger; meanDemand: real; 


PDDat aType: char); 
vat A023B,BRLDC,B011A,B019A,B023C,B023D,B073,M, PPV,B019,B021,BRLDCU: real; 


PD82strl: string(24]; 
PD82str2, PD82str3, PD¥2str4, PDS2strS, PNS2stro, PH&2str7, 
PDS2str8: string(255); 


PD86strl: string[24]J; 
PD8ostr2, PD8estr3, PD8estr4, PDSostrS, PD8estroe, PD8o6str7, 
PD86str8: string(255}; 


PD86str9: string[60}; 


infile,outfile:text; 
LTVar:real; 


t:integer; 


begin 
for t:=1 to numberOfOtrs do begin 
Gorexy (1,3); 
write('Quarter # °,t); 
mssion (rtnfile,’pd&21n.f11’); 
reset {inftile); 
read(infile,PD82strl, PD82str2, PD82str3, PD82str4, PD82st:5, PD82stro, 
PD82str7, PD82strs8); 
close (rnfrle); 


Begespe=troest (t}; {current quarterly forecast} 


1s. 





A023B:=meanDemand; 
if t-4 then begin 
A023B:=(observ(t-4]+observ(t-3)+observ(t-<|]+observ(t-1}) °4; 
end else 1f t = 4 then begin 
A023B:=(observ[t-3]+observ[t-2)+observ(t-1})/3; 
end else 1f t = 3 then begin 
A023B:=(observ([t-2]+observ[t-1])/2; 
end; 
1f AQ23B <= 0.0 then A023B:=1.0; 
strTemp: =copy (PD82st12,46,15); BOl1A:=StringToReal (StrTemp) ; 
B023C:=B011A*B023D; 
PPV: =B023C; 
delete (PD82str2,1,15); 
insert (NumToSti ing (A023B),PD82st12,1); 
delete (PD82str2,121,15); 
insert (NumToString(BO023D) ,PD82str2,121); 
Aelete (PDR2str2,106,15); 
inseit (NumToString(B023C) ,PD82str2,106); 
delere (PD82str5,91,15); 
insert (NumToString(PPV) ,PD82str5, 91); 
M:=mkCodeArry (tJ); {current mark code} 
delete (PD82str4,241,15); 
insert (NumToString(M),PD82str4,241); 
if (mkCodeArry(t] = 2) 0: {mkCodeArry[t]=4) then begin 
LTVar:=1.57*BO11A; 
BO19A: =BOLLA* (sqr (mad[t})*1.57)+(sqr(frest[t}))*LTVar; 
end else begin 
1f abs(B023C)< ERROR then B023C:=0.0; 
1f BO23C=0.0 then begin 


BOI9A: 


tl 
o 
Oo 


end else begin 
BO19A;: =COEFF2 *exp( POWER2 * 1n(B023C) ) 
end; 
end; 
delete (PD82str2,76,15); 
insert (NumToString(BO19A),PD82str2,76) ; 
if mkCodeArry[t] = 0 then begin 
BRLDC: =3; 
end else begin 
1f prbBrkPt = 0 then begin 
BRLDC: =5; 
end else begin 
if BO23C <« prbBrkPt then begin 


BRLDC: =4; 


end else begin 


BRLDC: 


a 
end; 
end; 
end; 
delete (PD82st12,16,15); 


insert (NumToString(BRLDC) ,PD82st12,16); 
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assign (outfile,'pdsZin. fil‘); 
rewiite (outfile); 
writein(outfile,PD82strl, PD82stiz, PD82str3, PH8zstr4. Ph¥2str5, PHB2stro, 
PO82sti7, PDd2str$); 
close (outfile); 
SwapVect ors; 
exec ('dA:\ulcp\PPDY2KRO.exe',‘d:\urep pds82in.fi1 pdB2out.fil ' j; 
SwapVector S; 
if DosErroi -- 6 then begin 
writeln; 
Sound (220): 
Aelay (300); 
NoSound; 
writeln ("Dos error #', DosError); 
HrtToCont; 
end; 
assign (infile, ‘pdsg2out.f11"'); 
reset (infile); 


read(infile,PD82st1:1, PD82str2, PD82str3, PD82str4, PD82strS, PD82stro, 





PD82st17, PD82str8); 
close (1infirle); 
strTemp: =copy (PD¥2st17,1960,15); BO19:=StringToReal (StrTemp); 
ROLeve lArry[t}:=B019; 
strTemp: =copy (PDB2str7,226,15); BO21:=StringToReal (StrTemp) ; 
EOOAr?y {t} : =B021; 
strTemp: =copy (PD82st17,121,15); BRLUCU:=StringToReal (StrTemp) ; 
if PDDataType = 'l‘ then begin 
InitPD8oFile; 
SwapVectors; 
exec ('dl: \uicp\PPD86KR4.exe',‘'d:\uicp pd8oin. fil pd8Boout. fil ° j; 
SwapVectors; 
if DosError <> 0 then begin 
writeln; 
Sound (220); 
delay (300); 
NoSound; 
writeln (‘Dos error #', DosError); 
HitToCont; 
end; 
assign (infile,''pd8oo0ut.fil'); 
reset (infile); 
read(infile,PD86str1, PD8ostr2, PD8ostr3, PD8estr:4, PDBestrS, PhH8oestro, 
PD8osti7, PD8ostr8, PH8ostr9); 
close (infile); 
st1Temp: =copy (PD8ostr8,1l66,15); SSADDBO[t]:=StringToReal (StrTemp) ; 
strTemp: =copy (PD8ostr#,181,15); SSADD[t]:=StringToReal (StrTemp) ; 
strTemp: =copy (PD8ostr8,1960,15); SSSMA[t} :=StringToReal (StrTemp) ; 
end; 


end; 





end; 
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Procedure ComputeERR{var ROLevelArry, EOQArry, frest,Mad, ERRArry:quarterArray; 
var mkCodeArry:qtrintArray; 
var qtr,OHCur: ,disPosals,disposalCount,qtrDispose: integer; 
ERRType: chat; 
unitPrice,orderCost,holdFrac,shortCost,salvRate, 


PLT, obsolRate,discRate,numY:sERR,milEssent:real); 


var WII,WI,ERR.TZe10o,dummy: real; 


pStockOut ,Z, LTD, LTVai , sigmaLTD, pafZ,probShort ,expShort:real; 


fDb] PrimeOfT, FPrimeOfT,Tn,Tnl]l,1,k,P,F,R,QO,C,Ps,M,t,delta:real: 


begin 
case ERRType of 

‘O': begin {no disposal} 
ERR: =OHCurr; 
ERRArry [qtr]:=ERRArry[qtr}) + 0; 

end; { case 0} 

‘l': begin {ulcp} 
Wie 24 eee St [(qGirai: 
W1l:= YRSERR * WII; 
if WI - MINERR then ERR: =WI 
else ERR: =MINERR; 
ERRArry[qtr):=ERRArry[qtr} + YRSERR; 


end; {case 1} 


'2': begin (net ben} 
if (frest[qtr] <> 0) and {EOQArry[qtr]) <>+0) then begin 
TZero:={{unitPrice - (unitPrice * salvRate)} + 
{orderCost / EOQArry[qtr]))) / (unitPrice * holdFrac)) + 
(EQOArry (qtr) / {8 * tresticqtr i; 
ERR<=TZero = 4" Crest iccri: 
ERRArry [qtr}:=ERRArry[qt1] + TZero; 
end {1} 
else begin 
ERR:=1; 
ERRArry [qtr}:=ERRArry[qtr} + 0; 
end; {else} 
end; {case 2} 
'3': begin {mod nb} 
pStockOut: =(holdFrac*unit Price) / 
({holdFrac*unitPrice)+(shortCost*milEssent)); 
Li=Z inv (potockOut) > 
ERD = fics tle tile e Live 
if (mkCodeArry[qtr} = 2) or (mkCodeArry[qtr]=4) then begin 
LiVarz=1],57* PLT; 
sigmaLTD: =sqrt (PLT* (sqr (mad[qtr])*1.57)+(sqr{frest [qtr})) 
*LTVar); 
end else begin 
if abs(LTD)< ERROR then LTD:=0.0; 
1f LTD=0.0 then begin 


sigmaLTD: =0.0 
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end else begin 
sigmaLTD:=sqrt (COEFF2 ‘exp (POWER2*in(LTD: | : 
end; 
end; 
pafZ:=ZPdf (Z); 
probShort: =utNormal (Z); 
expShort:=(LTD - ROLevelArry[qtr]) ‘probShort + sigmaLTD‘pcdfZ:;: 
lf (frest[(qti)] «- 0) and (EOQArry[qtr]) <« -0) then begin 
TZer10:=(unitPiice* (l-salvRate)}/(unitPrice*holdFrac) + 
EOCAIyiatr}/(2°4*frest[qtr]) + 
{orde:Cost+sho1tCost ‘expShort) / 
(EQOQArry [qtr] ‘unitPricetholdFrac) ; 
BRR = ieenOeg  frest (jt) > 
ERRA1 ry [qt:):=ERRArry[qtr] + TZero; 
end 1 2f) 
else begin 
ERR: =1; 
ERRArry [qt1r]:=ERRArry[qtr] + 0; 
end; {else} 
end; {case 3} 


fan2 begin 


Pe=infRate; k:=discRate; P:=unitPrice: F:=holdFrac; 
Cs=E00Arry(qtr): R:=4‘frest [qtr]; C:=ordercCost; delta:=90000; 
Ps:=unitPrice‘salvRate; M: =4*frest [qtr] ‘numY2rsOH; 


if (frest[qtr] < +» 0) and (EOQArry[qtr] - - 0) then begin 

ts=((unrtPrice - (unitPrice ‘ salvRate) + 
(orderCost / EOQQArry[qtr])) / (unitPrice * holdFrac)) + 
(EQOQArry[qtr] / (8 * frestiqtir])); 

JUS le & 

dummy : = (exp ( ((i-k)‘Q)/R)-1); 

fDb1PrimeOfT: =1; 

while (delta »0.01) and (dummy «<» 0) and (Tn » ERROR) 


and (abs(fDb]PrimeOfT) > ERROR) do begin 


fPrimeOfT:=((P*F*R)/(2*k)-(P*F*t*R) /2) texp(-k*t)+ 
((P*F*Q) /2+(P4Q* (1-k)} 4+C* (1-k)})}/ 
jexmut(i-k} *O) /R)=1)) expt {i-k) *t) 
-Ps*R-(P*F*R) / (2*k); 


fDb1lPrimeOfT:=P*F*R*((k*t-1)/2) *exp(-k‘t) + 
((P*F*Q* (i-k)) /2+(P*Q* sqr (1-k)+C* sqrv(i-k))/ 
texp(((1—k) *O) /R)-1))* exp ((i—k) *t); 

Tnl:=Tn-fPrimeOfT/fDb1PrimeOfT; 

delta:=abs(Tn1-Tn); 

G2=Tnl; 

Tn:=Tnl; 

dummy :={exp(({1-k)*Q)/R)-1); 


end; {while} 


if Tnl » ERROR then begin 


ERR:=Tn1*4*frest [qtr]; 


IAS, 


ERRArry [qtr]:=ERRArry([qtr] + Tnl; 
end {1f} 
else begin 
ERR:=1; 
ERRArry [qt2z]:=ERRArry(qtr] + 0; 
end; {else} 
end {1f) 
else begin 
ERR-=1- 
ERRArry [qtr):=ERRArry[qtr} + 0; 
end; {else} 
ena; {case 4} 
'S* ss begrh {trad} 
TZero:=In((salvRate* (discRate+obsolRate)+storRate'’ (1-obsolRate) * 
{(l+discRate)) /(discRate+obsolRate+storRate* 
(1-obsolRate) ‘(l+discRate)})/1n((1l-obsolRate) / 
{l+discRate)); 
ERR: =TZero*4*frest [qtr]; 
ERRArry (qt: ]:=ERRArry[qtr}) + TZero; 
end; {case 5} 
'6': begin { fixed) 
ERR: =numYrsERR*4* frest [qtr]; 
ERRArry [qtr]:=ERRArry[qtr} + numYrsERR; 
end; {case 6} 
end; {all cases} 
1f ERR » MINERR then ERR:=ERR 
else ERR: =MINERR; 
if OHCur:r - ERR then begin 
disposalCount:=disposalCount + 1; 
disPosals:=disPosals + (OHCurr - round(ERR)); 
qtrDispose:=OHCur: - round(ERR); 
OHCurr: =round (ERR) ; 
end; {if} 


end; {computeERR} 


procedure SDR(var OSHeap, BOHeap: Prior:tyQueueType; 
var wklyObserv:weeklyArray; 
var EOOArry, ROLevelArry, observe, frest, ERRArry:quarterArray; 
var numberOfOtrs,initInv,orderCount: integer; 
var disPosals,disposalCount:integer; 
meanDemand, rat ioPLTSTDMU, unit Price,orderCost ,holdFrac: real; 
shortCost,salvRate,PLT,obsolRate,discRate:real; 
var numYrsERR,numYrsOH: real; 
mi lEssent:real; 
var TWUS,endOH,endOS: integer; 
var ADDBO,ADD, SMA, Invest:real; 
wkDataType,qtrDataType, outputType, ERRType, anal Indtype: char; 
var totCost,holdTC,orderTC, shortTC, salvTR: real; 
var totCostArry,holdaTCArry,orderTCArry,shortTCArry, 


salvTRArry,investOtr,qtrSMA:quarterArray) ; 
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var wrlyBO,wklyOS:dataiecord; 
amt BO, amt kRecv, receipt ,wkhlyDemand, date, inrtOrders,1nitOS:integel 
wk ,qti,$12¢0S,s12eBO,qtiDispose,numbe10S,day: integer; 
landnorm, 1 andPLT,wklyInvest,qtrinvest,1epInvest,reorderTime:ieal: 
flagl,flag2: boolean; 
BOF111,4meaTot ,OSTot ,OSCu1: , BOTot , BOCurr, OHcur: , OHPrev, IPcurr, IPPrev: intege : 
cumBO, cumRO, cumHC, cumSR,orderlnterval:real; 


startint, intLength:real; 


begin 
Set Seed (seedArry {numberRep] ) ; 
OSCur1 : =0; 
OSTot :=OSCur? ; 
rnitInv:=round(numY1sOH’ frest [1] °4); 
InitializePriorityQueue (OSHeap); InitializeP1iorityQueue(BOHeap) ; 
if anallndType = '0' then begin 
InttInv:=round(EOOArry[1} + ROLevelArry[1]-frest[1]*PLT); 
numberOS: =round(PLT/ (EOQOArry[1])/frest[1})); 
1f numberOS «> 0 then begin 
for 1:= 1 to numberOS do begin 
wklyOS. Oty:=1ound (EOOArry[1]); 
if (PLT - (i-1) * (EOOArry[1])/frest({1]))) - 0 then begin 
day:=round(PLT - (1-1) * (EOOArry[1}/frest[(1})) * 13 + 1; 
wklyOS.Week:=day; 
Insert PriorityQueue {(OSHeap,wklyOS); 
OSTot:= OSTot + wklyOS.Oty; 
OSCurr:= OSCurr + wklyOS.Oty; 
end; {if} 
end; {for} 
end; {1f) 
end; {1f analind) 
1f (qtrDataType = ‘1') or (wkDataType = ‘1') then begin 


writeln(output file); 


writein{outputfile,'SDR Data Mioestaul lst Minhyee Opel Migw?, 2 indie ell (ele) fopvelsspga 9 .Okyetisr)) 4 
writeln(output file, ‘-----------~-------- nr en re rr er rr rr ren ren ener nen rr eres ncne= et 
end; 


OHCurrs=tnitinv; 
OHPrev:=OHCurr; 
BoCurr: =0; 
replnvest:= 0.0; 
BORIIN:= 0; 
iwusz= 0; 
ADDBO:= 0; 

Pbins= 0; 
cumBO: =0. 
cumRO:=0. 
cumHC: =0. 
cumSR:=0. 


dmdTot :=0; 
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disPosals:=0; 
disposalCount:=0; 
I PCurr : =OHCur?r+O0SCurr; 


VPPRev = Peurts 


foi qtr:= 1 to numberOfOtirs do begin 
if wkDataType = '1l' then begin 
writeln(output file); 
writeln(outputfile,'OTR WK REC DEM BO OS OH LPO ORDENT SOS ebook 
end; 
qtilInvest:= 0.0; 
qt rDispose: =0; 
wklyInvest:= 0.0; 
if (analIndType = '1l') and (qtr = 1) then 
Comput eERR (ROLevelArry, EOOArry, frest ,Mad, ERRArry ,mkCodeArry,qtr, 
OHCur: ,disPosals,disposalCount,qtrDispose, ERRType, 
unit Price,orderCost,holdFrac,shortCost,salvRate, PLT, 
obsolRate,discRate,numY1 sERR,m11lEssent); 
if (analIndType = '1') and (qtr «+ 1) then ERRArry [qtr] :=0; 
if ({(qti+l) mod 2) = 0) and (anallIndType = ‘0') then 
Comput eERR (ROLevelArry, EOOQArry, frest,Mad, ERRArry,mkCodeArry,qtr, 
OHCurr,disPosals,disposalCount,qtrDispose, ERRType, 
unit Price, orde:Cost ,holdFrac, shortCost,salvRate,PLT, 
obsolRate,discRate,numYrsERR,m1lEssent); 
rf (({qtr+1) mod 2) <-~+ 0) and (analiIndType = '0') then 
ERRArry (qtr) :=ERRArry (qtr-1]; 
for wk:= I to 13 do begin 
wk lyDemand: =round (wk lyObserv [date] ); 
amaTot:= dmdiTot + wklyDemand; 
receipt: =0; 
amt Recv: =0; 
amt BO: =0; 
wklyBO.Oty:=0; 
wk lyBO.Week: =date; 
wklyOS.QOty:=0; 


flagl:=FALSE; flag2:=FALSE; 


if not (EmptyPriorityQueue(OSHeap)) then begin {receive] 
repeat 
lf CurrWeek(OSHeap) = date then begin 

amt Recv: =Ext ractOty (OSHeap) ; 

receipt:=amtRecv; 

OSCurr:= OSCur:r - amtRecv; 

while (amtRecv - 0) and not (EmptyPriorityQueue(BOHeap)} do begin 

if CurrQty(BOHeap) <= amtRecv then begin 

amt BO: =CurrOty (BOHeap) ; 
amtRecv:= amtRecv - amtBO; 


BOCurr:= BOCurr - amtBO; 


BOFi1l1:= BOFi11 + amtBO; 


TWUS:= TWUS + (amtBo* (date - ExtractWeek (BOHeap}) } }; 


Loe 


TWUS'); 


end else begin 
BOHeap.HeapArray([1j).Oty:= BoOHeap.HeapArray({!|-vty - amtRec 
TWUS:= TWUS + famtRecv‘ ijdate - BOHeap.HeapArray[1] .Week))- 


BOCurr:= BOCurr - amtRecv; 


BOF111:= BOF:i1]1 + amtRecv; 
amtRecv:= 0; 
end; (sue); 


end; {while} 
OHPrev: =OHCurr; 
OHCurr:=OHCurr + amtRecv; 
end; 
1f EmptyPriorityQueue(OSHeap) then flag2:= TRUE 
else 1f currWeek(OSHeap) «= date then flagl:=TRUE; 
untrl flag] o: flag2; 


end: fif recelve} 


if wklyDemand - 0 then begin {1ssue} 
if wklyDemand ~ OHCur. then begin 
wklyBO.Oty:= wklyDemand - OHCurr; 
OHCurr:=0; 
Insert PriorityQueue (BOHeap,wk1lyBO) ; 
BOTot :=BOTot + wklyBO.Qty; 
BOCurr:=BOCurr + wklyBO.Qty; 
end {if} 
else OHCurr:= OHCurr - wklyDemand; 


end; {if 3rssue} 


PPPreyv:=]PCurr; {order} 
{ 1f wk=13 then begin } {for quarterly SDR) 
Preurr:= OHCur: + OSCurr - BOCurr; 
if IPCurr «<= ROLevelArry[qtr} then begin 
wklyOS.Qty:=round(ROLevelArry{gtr] + EOQArry[qtrj}) + BOCurr - 
‘OHCurr + OSCurr); 
randnorm:=- Normal; 
randPLT: =abs (PLT+ (randnorm*rat 10PLTSTDMU* PLT) ) ; 
1f randPLT + MAXPLT then begin 
randPLT: =MAXPLT; 
end else if randPLT <« MINPLT then begin 
randaPLT: =MINPLT 
end; 
wklyOS.Week:=date + round(randPLT‘’13) + 1; 
Insert PriorityQueue (OSHeap,wklyOS} ; 
OSTot:= OSTot + wklyOS.Qty; 
OSCurr:= OSCurr + wklyOS.Qty; 
orderCount:= orderCount + I; 
end; {if} 
ieends} {for quarterly SDR} 


if wkDataType 'l* then begin 


i 


writeln(output file,qtr:3,date:5,receipt:6,wklyDemand:6,BOCurr:6, 
oSCurr:6,OHCurr:6, [PCurr:6,orderCount:6,O0STot :6, BOTot:6, TWUS: 6) ; 


1f (outputType = '1') and ((wk mod 13) = 0) then begin 


5 


Hverocont: 
wiiteln(output file): 
end; (1) 


end; 


{test code only} 
{ $1Z2eBO:=51zeP1 1021 1 tyQueue (BOHeap) ; 


$12e0S:=S1zePriorityQueue (OSHeap) ; 


writeln{outputfile,'BO OQ Wk: ',currWeek(BOHeap):3,' BO O OTY: ‘,currOQty(BOHeap) :3, 


Sz: ',s1zeBO:3, ' OS OQ Week: ‘',currWeek(OSHeap):3,' OS 0 OTY: ' 
Currotcy (OSHeap) 23, "Sze. > sslzeony: 


writeln(outputfile); } 


REce late =r 
date:=date+1; 
wklyInvest:= wklyInvest + OSCurr + OHCurr; 
cumBO:=cumBO + ( (wkTWUS/52) ‘shortCost) ‘exp(-discRate/52‘date) ; } 
cumBO:=cumBO + ( (BOCurr/52) *shortCost) texp(-discRate/52‘*date); 
cumHC:=cumHC + (OHCurr* (holdFrac*unit Price) /52) *exp(-discRate/52‘*date) ; 
if wklyOS.Qty » O then 
cumRO:=cumRO + {unitPrice*wklyOS.Qty + orderCost) texm(-discRate/52*tdate) ; 
end; {for week] 
qtriInvest:= wklyInvest/13; 
investOtr[qtr] :=investOtr{(qtrj+qtrinvest; 
repInvest:= repInvest + qtriInvest; 
cumSR:=cumSR + (unitPrice’salvRate‘qtrDispose‘texpm(-discRate? (qtr-1)}/4)); 
totCostArry [qtr] :=totCostArry[qtr] + cumBO+cumRO+cumHC-cumSR; 
holdaTcArry [qtr]: =holdTCArry [qtr] + cumHC; 
orderTCArry [qtr] :=orderTCArry [qtr] + cumRO; 
shortTCArry [qtr]):=shortTCArry [qtr] + cumBo; 
salvTRArry [qtr] :=salvTRArry[qt:] + cumSR; 
1f BOFi11 <~ 0 then ADDBO: =7* (TWUS/BOFi11}; 
if amdTot <«- 0 then begin 
ADD: =7* (TWUS/dmdaTot } ; 
SMA:=1] - BOTot /dmdaTot ; 
qt rSMA [qtr] :=qtrSMA [qtr] +SMA; 
end; {if) 
if qtrDataType = ‘1' then begin 
1f (qtr=l) or (((qtr-1} mod 20 } = O}then begin 


writeln(output file); 


writeln(outputfile, 'OTR DMD OH IP os BO ADDBO ADD SMA 


end; 
1f (qtr »- 1) and (wkDataType = '1') and not(((qtr-1) mod 20) = 0) 
then begin 
writeln(output file); 
writeln(outputfile, ‘QTR DMD OH IP OSs BO ADDBO ADD 
end; 
end; 
if qtirDataType = 'l' then 


writeln(output file,qtr:3,observ[qgtr]):6:0,OHCurr:6,I1PCurr:6, 


oSCurr:6, BOCurr:6,ADDBO: 7:2, ADD: 72:2, SMA: 7:2, 9Cr Invest :9:2, 


eA 


INVEST 


INVEST 


DISP 


DISP 


ERR"); 


ERR‘); 


qtrDispose:6, ERRArzry [qtr] :9:2) ; 
rf (outputType = ‘'1') and (qtrDataType ='1') and 
not (qtr = 1) and {((lqtr-1) mod 20) = 0) then begin 
HrtToCont; 
writeln(output file); 
ere; {if} 
eee) (£02 qtr} 
Invest:= 1epInvest /numberOfOti s; 
endOH: =OHCurr; 
endos: =O0SCurr; 
totCost : =cumBO+cumRO+cumHC-cumsSR ; 
holATC: =cumHC; 
orderTC: =cumRO; 
short TC: =cumBo; 
salvTR: =cumSR; 


gotoXY(1,13); 
end; (sdr} 


procedure PrintHeader (prbBrkPt , seedIndex: integer ; 
salvRate,numYrsOH,ratioPLTSTDMU,meanDemand,varDemand: real ; 
var outputfrle:text; 
outputType,distrType, ERRType, anal IndType: char; 
out FileName: string; runDescript:descriptType; 
nmbrSteps,nmbrTrends: integer; 
stepMult,trendCoeff,trendPower: changeRealArry; 


startStep,starTrnd, endTrnd: changelIntArry) ; 


var 1:integer; 
errUsed, distrUsed, analUsed: string[7]; 
infrle:text; 
Year,Month,Day, Dayofweek:word; 


C028 : string[1]); 


A023B,B010,B011A,B020,B023C,B023D,B055,B057, B058,B061 ,B073,C008C, DO25E, 


MSLOD, SCR, TD,TSDRS, VO15R,V022,V101A,V102,V1034,V295: real; 


PD¥2stri: string[24]; 
PD82str2, PD¥2str3, PD82str4, PD82ser5, PD82stro, PDB2st17, 
- PO82str8: string[255]; 


begin 
distrUsed:=' Normal’; 
if distrType = ‘'2' then distrUsed:=*Poisson'; 
Syrused:=" UICP '; 
case ERRType of 
fo. : errUsed:=° None‘; 
*2': errUsed:=‘Net Ben'; 
‘3‘: errUsed:=' Mod NB‘; 


‘4*: errUsed:=' NPV NB‘; 


Bays 





Seep Uses: Traces 
"os er rUsed:=" Rye 
end; {case} 
analUsed:='Default'; 
case analIndType of 
‘l': analUsed:='UserSpc'; 
end; {case} 


1f outputType = '2' then begin 


wiiteln{output file, ' *** * outFrleName,' ***'); 


GetDate (Year ,Month,Day,Dayofweek) ; 
writeln(outputfile,' Date: ',Month,'-',Day,'-',Year); 
end; 
writeln(output file); 
writeln(outputfile,' Model: UICP - WILSON EOO '); 
writeln(output file); 
writeln(outputfile,' Description: ',runDescript); 
wirteln(output file); 
writeln{(outputfile,' Initzral simulation settings '); 


writeln(output file); 


writeln(outputfile,' Number: of quarters to simulate: 
writeln(output file,' Numbe: of replications of simulation to run: 
wirteln(outputfrle,' Random number generator seed type: 


1f seedType = 'l' then 

writeln{(outputfrle,' Random number seed start tndex: 
writeln(output frle,' Economic Retention Model: 
if ERRType = ‘o' then 

wrirteln(outputfile,' Number years economic retention used: 
writeln(outputfile.' Initial Inventory Type: 


if anallIndType = '1l' then 


writeln(outputfile,' Numbe1: years initral rnventory: 
writeln(outputfrle,' Type of demand distribution: 
wiiteln(outputfile,' Mean Demand: 
writeln(outputfile,' Vai Demand: 
writeln(outputfile,' Numbe1 of steps: 


if nmbrSteps -0 then begin 


for 1:=] to nmbrSteps do begin 


' mumberOfOtrs:5 ); 
' number OfRens:5); 


', seedtype); 


', seedindex:6); 


, errUsed); 


" numYrsERR:6:2); 


, analUsed); 


* NUuMl?sSOH: 6.2); 
, AistrUsed); 
' meanDemand:6:2); 
' varDemand:6:2); 


' mnmbrSteps:5); 


writeln{outputfile,' Sten? “,i:2," Sten Otrs “ Sstaresteniii- somal 
end; 
end; {rf} 
writeln(outputfile,' Number of trends: ‘ mmbrTrends:5); 


1f nmbrTrends +0 then begin 


for i:=] to nmbrTrends do begin 


writeln(outputfrle,' Trend:',i:2,° Start Qtr: *,starTirnd[:|]24,*° Ston Cty] “ enalinals | -47 


* Coeff: *, trendCoefflij:6:3,° Power: (trenaPowe: | 1] 675). 


end; 

end; {1if) 

writeln(output file); 

if outputType = ‘1l' then begin 
Hirt Tocane: 
clrser: 


end; 
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“ Stepmult (a}2oc3)s 


writelnioutput file.‘ Initial parameter settings 7 
assign (infile, ‘pag2in.f11'); 
reset {infirle); 
readiinfile, PD82sti1, PD82sti2, PD82sti3, PD¥2str4. PD¥2st:5, PD#2sti6, 
PD¥2st17, PD82st18&); 
close (infile); 
CO2Z8-=copy (PD82st:1,5.,1):; 
strTemp:=copy {PD82st:2.46,15); BOLLA:=StringToReal (StirTemp) ; 
strTemp:=copy({PD82st:2,91,15); BO20:=StringToReal (StrTemp) ; 
stiTemp: =copy (PD82sti2,121,15); BO23D:=StringToReal (Str Temp) ; 
strTemp: =copy (PD82st:2,181,15); BOSS:=St:1ingToReal (StrTemp); 
strTemp: =copy(PD82st:2,211,1!5); BOS7:=StringToReal (StrTemp) ; 
strTemp: =copy (PD82st:12,2260,15); BOS8:=StringToReal (StrTemp); 
strTemp: =copy (PD82st1r3,1,15); BOol:=StringToReal (StiTemp) ; 
strTemp: =copy (PD82str3,31,15); BO73:=StringToReal (St1Temp) ; 
strTemp: =copy (PD82st:3,760,15); CO08C:=<StringToReal (StrTemp) ; 
strTemp:=copy (PD8¥2str3,121,15); DO2ZSE:=StringToReal (StrTemp) ; 
strTemp: =copy (PD82str5,31,15); MSLOD:=StringToReal (StrTemp) ; 
strTemp:=copy (PD82st15,181,15); SCR: =StringToReal (StiTemp) ; 
strTemp: =copy (PN82st15,211,15); TO:=StringToReal (StrTemp) ; 
strTemp:=copy (PD82st15,226,15); TSDRS:=StringToReal (StrTemp) ; 
strTemp: =copy (PD82st15,241,15); VOI5SR:=StringToReal (StrTemp) ; 
strTemp: =copy (PD82st:16,16,15); V022:=StiingToReal (StrTemp); 
strTemp: =copy (PD8¥2str6,100,15); VIlO1A:=StringToReal (StrTemp) ; 
strTemp: =copy (PD82stro,121,15); V1l02:=StringToReal (StrTemp) ; 
strTemp: =copy (PD82stro,1360,15); V1034:=StringToReal (StrTemp}; 


strTemp: =copy (PD82str6,1660,15); V295:=StringToReal (StrTemp) ; 


writeln (output file); 


writeln {outputfile,' Prob Break = “",PrbBrkPt:38, * Min Risk.” = 277022: 9:2)); 
writeln (outputfile,' Shelf Life : Y AOA 5 © Max Risk :;: ‘,V102:8:2); 
writeln (outputfile,' Reqn Size :; *,B0O73:8:0, ' Ord Cost : ‘,V015R:8:2); 
writeln (outputfile,' Unite Price ©. “7 B055:8:2, : MSLOD 3 0 WEIL) 0)p G9 2) 
writeln (output file,‘ Salv Rate : ‘,salvRate:8:2, ; Proc Meth : ‘,D025E:8:0); 
writeln (outputfile,'‘ Preocin. ET ws “[BOLIA=3s:2, ~< Shortage <: ',V1034:8:2); 
writeln (output file,‘ Essential 9s 9,CO08C: 8:2, * R/O Low ee BOZO Be 2): 
writeln (outputfile,‘ Mfg Set-Up : ',B058:8:2, . RVORCOnNS tats me cS ore ee 
writeln (output file,‘ Obsol Rate : °*,B057:8:2, : Stor Rate : ',SCR:8:2); 
writeln (output file,’ Disc Rate —= “,B061:8:2, i Time Pref : ',V101A:8:2); 
writeln (output file,‘ Time SURS + *,TSDRS:8:2, egay OL war, Lose 05 
writeln (outputfile,' Init Yrs OH: *,numYrsOH:8:2, * PLT STD/MU: ‘,ratioPLTSTDMU: 8:0); 


ce de ne re 8 8 SS FS ATES ERTS ESSEC LORE ETERS EES ERE ECE ER OCR OER AEH OT) : 
if outputType = ‘1l‘ then begin 
Hit ToCont; 
e:rscr: 
end; 


end; {printheader} 


procedure DisplayPDOutput (var observ, frest, mad, EOQArry. ROLevelAiry, 


Sy 





SSADDBO, SSADD, SSSMA:quarterArtay; 
var stepIndAaArry, trndiIndArry,mkCodeArry:qtrIntArray; 
numberOfOtrs,initinv,repNum: integer; 


output Type:char) ; 
VRUE IEG iqhe(=te(aie 5 


begin 
wrrteln (outputfile); 
writeln(outputfile,'Replrcatiron Number ',repNum); 
writeln(output file); 
writeln(outputfrle,'PD8&2/86 Data'); 
writeln(outputfile, '-------------------------+------------ 5-5 + 2 = - = - oN 5 
for t:=l to numberOfOtrs do begin 
if (t=1) or (((t-1) mod 20) = O)then begin 
if (outputType='1') and (t-1) then HitToCont; 
wirteln(output file); 
writeln (outputfile, OTR OBS FRCST MAD Q R/O ADDBO ADD SMA MK ST TR’); 
end; 
writeln (output frle,ts3,observ([t]:6:0, frest(t]<8s2 magith<ss2, 
E0OAri y[t]):6:0,ROLevelArry[t]J:6:0, 
SSADDBO[t]:8:2,SSADD[t]:8:2,SSSMA[t]:0:2,mkCodeArry[t]:3, 
stepIndArry([t)]:3,trndIndArry [tJ:3); 
end; 
writeln (output frle); 
if outputType= '1l' then HitToCont; 


end; 


procedure DirsplayRepStats {var ADDBO, ADD, SMA, Invest, totCost:real; 
var orderCount ,d:sPosals,disposalCount ,endOH,endOS: integer; 
outputType:char ); 
begin 
rf numberRep = 1 then begin 
writeln(output file); 
writeln (output file, 1 Os 008 08 fee tie oie tie amet aie AeA BRAC SUR EONS ORR RAS EO eRe nee ne cence eee ee 
writeln(outputfrle,*Replitcation Final Statistics’); 
writeln(outputfile,* Num ADDBO ADD SMA Ords Invest EndOH EndoS DspCt TotDsp TotCost’); 
Snel ene: 
writeln(outputfile,numberRep: 4,ADDBO:7:2,ADD:6:2,SMA:6:2,orderCount:6, Invest:8:2, 
endoOH:6,end0S:6,disposalCount:6,disPosals:7,totcost:14:2); 
tf numberRep = numberOfReps then 
writeln (output file, §** Feet eee ete eee eee ee ee eee ee ee 
rf outputType ='l' then begin 
delay (1500); 
clrscr: 
end; 


end; 


procedure DoStats(var currMean,currVar,sampleReal:real; 


var samplelInt:integer; 
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var confiInt:real; 


numbe: Rep: integer) ; 
var sample,oldMean,oldVar:real; 


begin 

if sampleReal = -9999.0 then sample:=samplelInt 

else sample:=sampleReal ; 

oldMean:=currMean; 

oldVar:=curr Val; 

1f numberRep<=1] then currMean:=sample 

else CurrMean:= (((numberRep - 1) * oldMean)+sample: /numbe1Rep; 

1f numberKep< =2 then currVar:= 0.0 

else currVar:= (((numberRep-2) 'oldVar) + ( (number Rep- 1) *SOR(oldMean} | - 
(number Rep‘ SOR {currMean))+SOR(sample))/(numberRep-1); 

1f mumberRep - 0 then confint:= 1.96 * SORT(currVar/numberRep) 

else contint:=0.0; 


end; {dostats} 


procedure DisplaySrmStats (var simADDBO,simADD,simSMA,sirmInvest,simTotCost, 
simOrderCount,simDisposals,simD1sposalCount, 
simEndOH,simEndoS,c1iADDBO,c1ADD,ciSMA,cilnvest, 
ciTotCost,c101rderCount,ciDisposals, 
cidisposalCount,ciEndOH,ciEndoS:1eal; 
output Type: char; 
hourl,minutel,second!],hdSecl,hour2,minute2, 
second2, hdSec2: word); 
vai upADDBO,upADD,upSMA, upInvest ,upOrderCount ,upDisposals,upDisposalCount, 
lwADDBO, 1wADD, lwSMA, lwInvest , lwOrderCount, lwDisposals,]lwDisposalCount, 


lwEndOH, lwEndoS, upEndOH, upEndoSs: real; 


begin 

{ uPADDBO:=simADDBO+ciADDBO; lwADDBO:=s1mADDBO-c1ADDBO; 
upADD: =simADD+c1ADD; 1lwADD:=simADD-ciADD; 
upSMA:=simSMA+ciSMA; lwSMA:=simSMA-c1SMA; 
upEndoOH : =simEndOH+ci EndOH; lwEndOH: =simEndOH-c1 EndOH; 
upEndoS: =simEnd0OS+c1Endos; lwEndOS: =simEndOS-ciEndos; 
upInvest:=simInvest+cilnvest; lwInvest:=simInvest-cilnvest; 
upOrderCount :=simOrderCount+ciOrderCount; 
lwOrderCount : =srmOrderCount -ciOrderCount ; 
upDisposals:=simDisposals+cidisposals; 
lwDisposals:=simD1isposals-ciDisposals; 
upDisposalcount:=simDisposalCount+c1iDisposalCount; 
lwDisposalCount:=simD1sposalCount-ciD1sposalCount; 
1f IwADDBO < 0.0 then IwADDBO:=0.0; 
if lwADD - 0.0 then lwADD:=0.0; 
if lwSMA - 0.0 then lwSMA:=0.0; 
if iwInvest < 0.0 then Iwinvest:=0.0; 


1f IlwOrderCount < 0.0 then lwOrderCount:=0.0; 


IES 2) 





if lwDisposals « Ji.) then lwDisposals:=0.0; 

if lwDisposalCount 0.0 then lwDisposalCount:=0.u; 
if lwEndOH - 0.0 then lwEndOH: =0.0; 

if lwEndOS - 0.0 then lwEndOS:=0.0;} 

writeln{(output file); 

wiiteln (output file, tte ee ee ee eee ee eee ee ee ee eee aie eet en 


writelni(outputfile,’Simulation Final Statistics’); 


writeln(outputfile,'Final Means and Confidence Interval (95%)'); 


writeln(outputfile,' Mean liane 
writeln(outputfile,' ADDBO SMA OBO GLAD Melton is 
writeln(outputfile,' ADD ", SIMADD2 12:72, C1 AUD: 1222) > 
writeln(outputfile,' SMA ', S1mSMA:12:2,¢c1SMA:12:2); 
writelnioutputfile,' ORDERCOUNT » -srmOrdercount sl 222 ncrordencount: lc aie. 
writeln(output file.' INVEST ', simInvest:l12:2,cilInvest:12:2); 
writeln(outputfrle,' ENDING OH DP SIMENGOH 2 s2rcMeEnao rezk 
writeln(outputfrle,' ENDING OS ", SImENdOS: 12:2,c1EndOS: 12:2); 


* 


writeln(outputfile,' DISPOSAL COUNT ',simDisposalCount:12:2,¢1D1sposalCount:12:2); 
writeln(outputfile,' DISPOSALS ',srmDisposals:12:2,c1Disposals: 12:2); 
writeln(output file,' TOTAL COST ‘ Simlot€ost sl2:2, crlotGoct-12-2)¢ 
writeln(output file) 8 808i e eee ee eater eee ee ee en ne ae eee 
writeln(output file); 

writeln(outputfrle,'Sim Start Time * nourl, = minutel,* =" second. . +o naseci }- 
writeln(outputfile,'Sxrm End Time ‘, hour2,°s" minute?, °:*;sécona2, > ,ndsec2); 


1f outputType ='1' then HitToCont; 


end; {displaysimstat} 


Procedure DisplayQtrArrysi(var totCostArry,holdTCArry,orderTCArry,ERRArry, 
shortTCaArry, salvTRArry: quarterArray; 
numbe1OfOtrs:integer); 

VAR qtr:integer; 
begin 

writeln (output file, * 8808645 88 8 Pm ata tial t amit R88 eee E SRR SO RO 8 SOR e SO eee Re ee oe ee 

writeln{outputfile); 

writeln(outputfile,' Quarter cummulative costs and years ERR for graphing'); 

writeln{output file); 

writeln{output file, ' QTR TOTAL HOLD ORDER SHORT SALVAGE ERR‘); 

for qtr := 1 to number1OfQtrs do 

writeln(outputfile,qtr:4,totCostarry [qtr):12:2,holdTCArry[qtr] :12:2, 
orderTCArry (qtr) 21222, snortTecarry (atrlei2- 2, 
salvTRArry [qtr] :12:2, ERRArry [qtr]:10:2); 

writeln(output file); 

writeln(outputfile,' Quarter: SMA and Invest for steady state graphing'); 

writeln(output file); 

writeln(outputfile, ' QTR SMA Invest'); 

for qtr := 1 to numberOfQtrs do 

writeln(output file,qtr:4,qtrSMA[qtr]:12:2,investOtr[qtr):12:2); 


end; {displayqtrarray} 
begin (main) 


textcolor(14); 


stom: =FALSE; 
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simCount :=0 
currSeer): =() 
noReal :=-999%,.0; 
Homme s=0 + 
randSeedArry (seedArry}; 
Frontscreen; 
Runtype {distrType, outputType,wkDataType, qtrDataType, PUDataType, 
repStat Type, ERRType, anallIndType,numberO0fOtrs,numberOfWks,number1OfReps, 
seedindex,meanDemand, varDemand, numY1 sOH,numYrsERR, input file,output frie, 
frest,mad, seeds, out FileName, runDescript); 
repeat 
rewrite (outputfile); 
simCount :=simCount +1; 
GetTime( hourl,minutel,secondl ,hdSec1); 


1 to numberOfReps do begin 


for numberRep : 
1f seedType = ‘i' then begin 
if numberRep = 1 then begin 
fo: s:= 1 to seediIndex do cu:1Seed: =GetNext Seed (currSeed) ; 
SetSeed (currSeed) ; 
end {2f} 
else begin 
currSeed: =GetNext Seed {(currSeed) ; 
Set Seed (currSeed) ; 
end; {lelse} 
end ({2f} 
else SetSeed(seeds[numberRep] ); 


InitializeArrays ({observ,EOQArry,ROLevelArry,SSADDBO,SSADD,SSSMA,ERRArry, 
stepIndArry, trndindArry,mkCodeArry,numberOfOtrs, 
numberOfWks,numberRep,meanDemand, 
wklyObserv,meanDmdArry,varDmdArry,totCostArry, 
holdTCArry,orderTCArry, shortTCArry,salvTRArry, 
investOtr,qtrSMA) ; 

LoadObserv (observ, frcst,mad,wklyObserv,meanDmdArry ,varDmdArry, 
observType,distrType, numberOfOtrs,numberOfWks,numberRep, 
simCount,trendOn,stepOn,nmbrSteps, nmbrTrends, 
meanDemand, varDemand, input file,seeds,startstep, 
startrnd, endtrnd,stepmult, trendcoeff, trencipower); 

1f numberRep = 1 then begin 

rf simCount=1 then InitPD82File (prbBrkPt,numYrsERR,salvRate, 
numYrsOH, rat ioPLTSTDMU, storRate, 
obsolRate,discRate,infRate,milEssent); 
PD82Edit (prbBrkPt,unitPrice, PLT, orderCost ,holdFrac, 
shortCost,salvRate,numYrsOH, rat :oPLTSTDMU, numYrsERR, 
storRate,obsolRate,discRate,infRate,milEssent) ; 

end; 

if numberRep=1 then PrintHeader (prbBrkPt,seedindex, salvRate,numYrsOH, 

rat 1oPLTSTDMU, meanDemand, varDemand, 

output frlie,outputType,distrType, 

ERRType, anal IndType, out FileName, runDescript, 
nmbrSteps,nmbrTrends, stepMult, 


trendCoeff,trencdPower,startStep, 
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starTrna, endTrnd) 
Forecast jiobserv,frest.mad, stepinedArry, trndindArry. 
mkCodeArry ,numberOfOtis,numberRepn,unitPrice}) ; 
LoadiLevels (frest, mad, observ, EOOQArry, ROLevelArry, SSADDBO, SSADD, SSSMA. 


mkCodeArry, numberOfOtrs, prbBrkPt, meanDemand, PDDataType) ; 


if PDDataType='1' then DisplayPDOutput (observ, frest, mad, EOQArry, 
ROLevelArry, SSADDBO, SSADD, 
SSSMA, stepIndArry, trndindArry, 
mkCodeArry, number OfQOtrs,initinv, 
numberRep, outputTyne) ; 

SDR (OSHeap, BOHeap, wk lyObserv, EOOArry, ROLeve lArry,observ,frest, 
ERRArry,numberOfOtis,initInv,o1derCount,disPosals,disposalCount, 
meanDemand, rat :oPLTSTDMU, unit Price,orderCost ,holdFrac, shortCost, 
salvRate,PLT,obsolRate,discRate,numYrsERR,numYrsOH,m:1Essent , TWUS, enciOH, endOS, ADDBO, 
AUD,SMA, Invest ,wkDataType, qtrDataType, outputType, ERRType, analIndType,totCost, 
holaTC, orderTC, shortTC, salvTR,totCostArry,holdTCArry, 
orderTCArry,shortTCArry, salvTRArry, investQtr,qt1SMA) ; 

if repStatType = ‘l' then DisplayRepStats (ADDBO,ADD,SMA,Invest,totCost, 

orderCount,disPosals, 
AisposalCount,endoH, 
endoS ,outputType) ; 

1f numberRep = 1 then begin 

simADDBO:=0.0; simADD:=0.0; simSMA:=0.0; simInvest:=0.0; 
simOrderCount:=0.0; simDisposals:=0.0; srmDrsposalCount: =0.0; 
simEndOH:=0.0; simEndOS:=0.0; simTotCost:=0.0; simHoldaTC:=0.0; 
simOrde1rTC:=0.0; simShortTC:=0.0; simSalvTR:=0.0; 

endy tte) 

DoStats{s1mADDBO, varADDBO,ADDBO,noInt ,c1ADDBO, numberRep) ; 

DoStats(simADD, varADD,ADD,noInt,ciADD,numberRep) ; 

DoStats(simSMA,va1rSMA,SMA,noInt,ciSMA,numberRep) ; 

DoStats(simInvest ,varInvest, Invest ,noInt,cilnvest,numberRep) ; 

DoStats(simOrderCount ,varOrderCount , noReal,OrderCount ,c1OrderCount , 

numberRep) ; 

DoStats(simDisposals,varDisposals,noReal ,Disposals,ciDisposals,numberRep) ; 

DoStats({s1mEndOH, var EndOH, noReal , endOH, ci EndOH, numberRep) ; 

DoStats(simD1sposalCount,varDrsposalCount ,noReal,disPosalCount, 

ciDisposalCount ,numberRep) ; 

DoStats(s:mEndOS, varEndOS,noReal ,endOS,ciEndOS,numberRep) ; 

DoStats(simTotCost,varTotCost,totCost,noInt,ciTotCost ,numberRep) ; 

DoStats{simHoldaTC, varHoldTC, holdTC,nolInt ,ciHoldTC, numberRep) ; 

DoStats(simOrderTC, varOrderTC, orderTC, noInt,ciOrderTC, numberRep) ; 

DoStats(simShortTC,varShortTC, shortTC, noInt,ciShortTC,numberRep) ; 

DoStats(simSalvTR,varSalvTR,salvTR,nolInt,ciSalvTR,numberRep) ; 

end; {for) 
for i:= 1 to numberOfQtrs do begin 

qtrSMA[i]}:=qtrSMA[1i])/numberOfReps; 

investOtr[1]:=investOtr[i}/numberOfReps; 

end; {for) 
Get Time (hour2,minute2, second2,hdSec2) ; 


DisplaySimStats{simADDBO,simADI),simSMA,simInvest,simTotCost,simOrderCount, 
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simDisposais,s:mD:i1sposalCount ,simEndOH, simEndOS,c1ADDBC: 
CiADD,ci:SMA,ciInvest,c:TotCost,c:0OrderCount,ciDisposals, 
c:DrsposalCount ,c1EndOH,c1EndOS,output Type,hour!, 


minutel,second!] ,hdSec] ,hourz,minute2, seconde, hdSecc, ; 














numQt:r := 1 to numberOfOtrs do begin 

—totCostArry [numOtr) :=totCostArry [numOtr) /numberOfReps;: 
holdaTcArry [numotr ] :=holdTCarry [numOtr] /numberOf Reps: 
order TCArry[numOti):=s01der1TCArry[numOt: ) /numberOfReps; 
short TCArry [numOt:}:=shortTCArry [numOtr} /numberOfReps; 
~ salvTRArry [numOtr):=salvTRArry (numOtr) /numberOfReps; 
ERRArry [numOt1 ) :=ERRArry [numOtr) /numberOfReps; 

end; (for) 

‘DisplayOtiArrys(totCostarry ,holaTCArry ,orderTCArry , ERRALTY , 
; short TCArry, salvTRArry,numberOfOtrs) ; 
close {output file); 

RunAgain (out put file, runDeseript,outputTynpe, ERRType, stop, 
numY: SERR,outFiieName) ; 

1 stop; 


feolor (15); 


{main program UICP-Simulator} 
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Unit TOOLBOX; 


a ee el 


‘This Untt piovides a toolbox of useful functions functions and 


*procedures for data input. ‘ 


ee Oe ene i eC CuCaLE CC Oa Oe Ce OO OOOO AOA Ie SOOO OOD oO] 


Interface 
Uses CRT; 


type pdsb2field=string[15]; 


» 


var strTemp:pds2field; 


function Get_Answe1r: boolean; 

procedure HitToCont ; 

function Get_Intege:r (low,high: integer) :integer; 
function Get_Real (low, hrgh:real):real; 

function NumToString (var value:real) :pd82field; 
function StringToReal (var S:pd82field):real; 


function Get_LongInt (low,hrgh:longint):longint; 


Implementation 


functron Get_Answer; {Returns a Boolean result for a yes/no query} 


val Char_In:Char; 


Correct: Boolean; 


begin 
Correct :=False; 
repeat 
Char_In:=ReadKey; 
write (Char_In); 
case Char_In of 
"Y< 'v* > bean 
Wreteeailig if CGS js 
Get_Answer:=True; 
Correct :=True 
end; 
“N* on’ = beorn 
writelm ("6*%); 
Get_Answer:=False; 
Correct :=True 
end; 


else begrn 
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writelin; 
Sound(220} ; 
delay (300); 
NoSound; 
writeln ('** Un-recognizable answer ‘**'); 
wrivcein (*Enter Y or N,'); 
wrifeln ('Re-enter your answer: ') 
end 
end; (case} 
until Correct; 


end; {Get_Answer} 














Procedure HitToCont; 
var dummy: char; 


begin 

writeln; 

write (' Hae any key to continue |...*); 
‘dummy : =readkey; 


end; 


{Gets an integer input between low and high, prompts until one is received) 


function Get_Integer (low,high:integer):integer; 


var numberString: string[10]; 


error, numberValue: integer; 


begin 
repeat 
readin {numberString) ; 
val (numberString, numberValue, error); 
1f error <> 0 then begin 
writeln; 
Sound (220); 
delay (300); 
NoSound; 
write ('*** Invalid number, enter an integer: '}) 
end else if (numberValue- low) or {numberValue-high) then begin 
writeln; 
Sound (220); 
delay {300); 
NoSound; 
writeln ('*** Invalid Range - value must be a positive integer’); 
write {'between ',low,’ and ',high,’ Enter number: '); 
error:=); 
end; 
until error=0; 


Get_Integer:=numberValue; 
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end; {function) 


{Gets an longint input between low and high, prompts until one 1s received} 


function Get_LongInt (low,high: longint):longint; 


var numberString: string[10]; 
error: integer; 


numberValue: longint; 


begin 
repeat 
readln (numberString) ; 
val (numberString. numberValue, error); 
if error «+ 0 then begin 
writeln; 
Sound (220) ; 
delay (300); 
NoSound; 
write ('*** Invalid number, enter an integer: '} 
end else if (mumberValue< low) or (numberValue-high) then begin 
writeln; 
Sound (220) ; 
delay (300); 
NoSound; 
writeln (‘*** Invalid Range - value must be a positive integer'); 
write (‘between ‘,low,' and ‘,high,' Ente: number: '‘'); 
jen cypress 
end; 
unt:1] error=0; 
Get_LongInt:=numberValue; 


end; {function} 


{Gets a real value between low and high, prompts until one is received} 


function Get_Real (low,high:real):real; 


var Number_String: string; 
Error: integer; 


Number _Value:real; 


begin 
repeat 

readin (Number_String}; 

val (Number_String, Number_Value, Error); 

if Error <> 0 then begin 
Sound(220); 
delay (300); 
NoSound; 
writeln (‘**You must enter a valid real number** '); 

end else if (Number_Value< low) or (Number_Value-high) then begin 


writeln; 
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Sound (zz) = 
delay (300) 


~NoSound; 


elioh. =): 


]}:=Number_Value; 


r_ Real} 


yr it=l to 16 do 
Sir) = ° * then S[i):='0' 


Meise 1f S{i) = *.*' then delete (S,i,1); 


yn StringToReal (var S:pd82field):1eal; 


R2: real; 


@eer1o12); 


‘oReal :=R1+(R2/100000000) ; 


{Unit Toolbox} 


write (‘between ',low:0:1,' and ',high:0:1,' 


NumToString (var value:real) :pd82field; 


Enter 


writeln ('*** Invalid Range - value must be a reai value'}; 


number: 
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ees 


unit unirand; 


inte1 face 

type seedArryType = array [1..1000} of longint; 
var seeds, seedArry: seedArryType; 

procedure SetSeed (seed: longint}; 

function Get Seed: longint; 

function GetNextSeed (lastSeed: longint):longint; 
function RandomUni form: real; 

procedure randseedArry (var seedArry:seedArryType) ; 
function GetPoisson(vai meanDemand:real):integer; 
function GetNormal:1real; 

function GetGeometric(p:1eal):integer; 

function GetNegBini(p:real;s:integer) : integer; 
function GetUniformInt (high: integer): integer; 
function ZInv (p:real):real; 


function ZPdf (Z:real):real; 


function utNormal (Z:1eal):real; 


implementation 

var a:longint; 

procedure SetSeed (seed: longint); 
begin 


a:=seeqd 


end; (procedure} 


function GetSeed:longint; 
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begin 
Getbeed:=a 


end; {proceduie) 


function RandomUni form: real; 


const B2E15: longint=32768; 
B2els: Jongint=05536; 
Modlus: longint=2147483047; 
Mult]: longint=24112; 


Maltz: longint =26143; 
vai H115,Hi31,Low15, Lowprd, Ovflow,Z1:longint; 


begin 

Z1:=a; 

HvlS2*Z1 div B2Elo; 

Dowprd:=(Zi - HilS * B2El6) * Multi; 

Low15:=Lowprd div B2El6; 

Hi31]:=H115 * Multl + Low15; 

Ovflow:=Hi31l div B2E15; 

Z1:=(({Lowpid - LowlS * B2El6) - Modlus) + 
(Hi31 - Ovflow * B2E15) * B2Elo) + Ovflow; 

Poeza< O then Z1:= Z1 + Modlus; 

Seee= 92) Giv B2Elo; 

Lowprd:= (Z1 - H115 * B2El6) * Mult2; 

Low15:=Lowprd div B2El6; 

Meee 115 * Mult2 + Lowl5; 

Ovilow:= Hi31 div B2E15; 

Zizs=(((Lowprd - LowlS * B2El6) - Modlus) + 
(H131 - Ovflow * B2E15) ‘ B2El6) + Ovflow; 

1f Z1 - O then Z1:= Z1 + Modlus; 

a:=Z1; 

Pemcomunitorm:= (2 * (Z1 div 256) + 1) / 160777210.0; 


end; 
function GetNextSeed (lastSeed: longint):longint; 


const M:extended=2147483647.0; 
a:extended=715.0; 
b:extended=1058.0; 


c:extended=1385.0; 
var Z:extended; 


begin 
Z:=lastSeed; 
rf lastSeed=0 then begin 
@e=1973272912.0; 


GetNextSeed: =round (Z) ; 
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end else begin 
Fag r4)) /f [ule 
ZS — NOUNG (UP o)eeeMs 
Z:=(B*2) / M; 
Z:=(Z-r1o0und(Z-0.5)) *M; 
Z:=(C*Z) / M;: 
Z:=(Z-round(Z-0.5})) *M; 
Get Next Seed: =10und (Z) ; 
end; 


end; {getnext seed) 


function GetPoisson({va1 meanDemand:real):integer; 


vat alpha,beta, Ul:real; 


lz:integer; 


repeat 
| ae 
alpha:= exp (-meanDemand) ; 
U1: =RandomUni form; 
beta:=beta‘UI1; 

until beta<alpha; 

Get Poisson:=1 


end; 
function GetNoirmal:real; 
vat U1,U2,V1,V2,W, Y:real; 


begin 

repeat 
U1:=RandomUni form; 
U2: =RandomUni form; 
Vise2*ar=15 2 V2s—2°U2-1; 
W:=sqr(V1l)+sqr (V2); 

MBE a) Wee 0 

Y:=sqrt ((-2¢In(W)) /W); 

GetNormal:=V1*Y; 


end; 
function GetGeometric(p:real):integet; 


var U:real; 


r:integer; 


ISO 





form; 











feerics=\; 


yet NegBin(p:real;s:integer):intege: ; 


‘ormint:=round( {high-1) ‘RandomUni form) +1; 


Zinv (p:real):real; 


tr (-2*Ini(p)); 
Pete 1 5517+0.802853*t+0.010328* sqr(t))/ 


(1+1.4327884t+0.1892609'* sqr(t)+0.001308*exp(3*in(t)))); 


ZPdf (Z:real):real; 


989 "exp(-(sqr(Z)/2)):; 


ie) 
Normal (Z:real):real; 
tArry= array [0..3] of real; 


,OsubJ: const antArry; 


J, sumOsubJ,R1X,erfX,X:real; 


9795523953175; 


ake 


PsubJ [1 )t=21. 27926 le le234 feo: 

PsubJ [2] :=6. 9963834006 191955; 

PsubJ (3) :=-0.0356009843701815385; 

OsubJ [0] :=215.0588758698oOl2; 

OsubJ [1] :=91.1¢04905404514901; 

QsubJ (2) :=15.082797630407787; 

QOsubJ[3]:=1.0000000000000; 

sumPsubJ:=0.0; 

sumQsubJ: =0. 0; 

A eavAi/ E¥gteic (741) 5 

1 X = 0.0 then X:=#0.000001; 

if X -« 0.0 then x:xabs(X); 

for 7:= 0) to sido begin 
sumPsubJ:=sumPsubJ + PsubJ[j)] * exp((2°3) ‘1n{X)); 
sumQsubJ:=sumOsubJ + OsubJ[j] * exp((2*j)*In{X)); 

end; (for} 

R1X:=sumPsubJ/sumOsubJ; 

erfX:=X*R1X; 

1f Z »= 0 then utNormal:=1 - ((l+er£X)/2) 

else utNormal:=({l+erfX) /2; 


end; 


end. {Unit Unirand} 


IE3y2 





unit PDUnit; 


Inter face 


uses dos, crt, toolbox; 


Val prbBrkPt :integer; 
unit P11ce,PLT,ordercCost ,hnoldFrac,shortUcvst: real; 
numYrsERK,salvRate,numY:sOH,1atioPLTSTDMU :real; 


storRate,obsolRate,discRate,infRate,milEssent:1eal; 


procedure Init PD82File (var prbBrkPt:inteqer: 
var numYrsERR, salvRate,. numYrsOH, ratioPLTSTDMU, storRate, 


obsolRate,discRate,infRate,milEssent: real); 


procedure PD82Edit (var prbBrkPt:integer; 
var unitPrice,PLT,orderCost,holdFrac, 
shortCost,salvRate,numYrsOH, 
rat 1oPLTSTDMU, numYrsERR,storRate,obsolRate, 


discRate, nfRate,m1rlEssent:real); 
procedure InitPD#oFile; 
Implementation 


procedure Init PD#2Firle (var prbBrkPt: integel; 
var numYrsERR,salvRate, numYrsOH, ratioPLTSTDMU,storRate, 


obsolRate,discRate,infRate,milEssent: real); 


var AAC,AL,B067A,B0607G,C028,DRLI,DO31C,D125N, ERRI,F024,HODI,MARLI,PVPI,RII,RO, 
YR7POC,Y006A,Y006B,EOOIND, PVUI : char; 


Dio, FILLER : string [2]: 


A023B,BRLDC,B010,B011A,B012F,B019A,B020,B023C,B023D, BO23F,B023H,BG, BOSS, 
B05S5A,B057,B058,B0S8A,B061 ,B070,B073,B093,B280,C008C, DOPTC, DTC, DO25E, 
F009,HOD,HO141,HO0142,H01I43,HO144,H0145,HO0146,H0147,HOI48,H0149,HOI4I0, 
HO1411,H01412,H01413,HO1414,HO1415,H0O14160,HO1417,HO1418,HO1419,H01420, 
ILR, IMECY,M,MOQOOAD, MSLOAD,MSLOD,NRFIDRT,OSQ, PDO, PPV,QDH, RFIDRT, RIYAYABY, 
RSV,RT,SCR,SSOH, TD, TSDRS, VO15R,V016,V022,V039,V041R,V042R,V043R,V044, 
V101A,V102,V1034,V108,V295,LILT, LILY, PCR3,01B,02B,RMNAST,SER, YDR,MNOOAD, 
APSR, ARCI, BOO, BRLCI , BRLDCU, BRLO, BRPLO, BRO, B014A,B019,B019B,B021,BO021A, 


ERR ,MONDO ,OOCI, POC, PPVBNDO, PZO,RCI,RLCI,RPLCI,ROC1,VPSR : real; 
POB2strl: string{24]; 
PD82sti2, PD¥2str13, PD82str4, PD82stxrS, PHO¥2sti6, PD82sti7, 


PDSZ2str8: string[{255]; 


outfrle:text: 


leat, 





{initialization values} 


AAC: ='N‘';: Ab:='N'; BOo7A:s='N'; BObO7G:='N'; CO28:='0'; DRLI:s='N'; DO31E:s=' |; 


DIZO0:="O0o0"; DIZSNr=" *; ERRIse*N*; FO24:=' *; HOD]:=' "> MARL] s=*Y"> 
PVPITs="Y*; Rils="N’; ROs=*N*'? YRIPOCs= * "s YOOGAs="N"'; YOO0cE:="H*-; 
EQOINDS="N*'; -PVOD2=. "7 FICEER = 

Ue ches il. (Os {system requisition average} 

BRLDC: =5.0; {basic reorder level distribution code} 
BOGE =OR0- {contract prod lead time} 

BO11A:=8.0; {contract proc lead time} 

BOVE —0n 0 

BO]9A;=20..0; {non cred group proc variance} 

BOZO =r {system reorder level low limrt qty} 
BO23D:=1.0; {gross sys demand end of lead time} 

BO23C: =B011A*B023D; {gross sys demand during lead time} 


° 


BO2Z3F:=0.0; BO23H:=0.0; BG:=0.0; 


B055:=100.00; {unit price} 

BOSSA: =0.0; 

BOS7:=0.12; obsolRate:=B057; {obsolescence rate} 
BOS58:=600.0; {manufac set-up costs) 


BOS58A:=0.0; 

BO61:=1.0; {drscount rate} 

BO70:=0.0; 

BO73:s=).0; {expected units per requisition) 

BO093:=0.0; B280:=0.0; 

Co08C:=0.5; {average item essentiality} 
DOFTC:=0.0- Dire :—0)..0; 

DO25E:=0.0; {procurement method} 

FO09:=0.0; HOD:=0.0; H0O141:=0.0; HO142:=0.0; HO143:=0.0; HO144:=0.0; 
HO145:=0.0; HO146:=0.0; HO147:=0.0; HO148:=0.0; HO149:=0.0: HO1410:=0.0; 


HO1411:=0.0; HO1412:=0.0; HO1413:=0.0; HO1414:=0.0; H01415:=0.0; H01416:=0.0; 


HO1417:=0.0; HO1418:=0.0; HO1419:=0.0; H01420:=0.0; ILR:=0.0; IMECY:=0.0; 


M:=1.0; {mark code} 

MOQOAD: =6.0; {max order qty attritron qtrs demand} 
MSLOQOAD: =99.0; {max number safety level gqtrs attrition} 
MSLOD:=2020- {max number of safety level qtrs demand} 


NRFIDRT: =0.0; 


CSO7=0e0r {non-parametric order stat qtrs} 
PHO: =8 50° {past qtrs demand} 

PPV: =BO023D*BO11A; {proc problem var (mean) } 

QDH: =0.0; {quarters demand history) 


RFIDRT:=0.0; RIYAYABY:=0.0; 


RSV; =0.. 0; {requisition size variance} 
RT: =0.0; 
SCR:=0.01; storRate:=SCR; {storage cost rate} 


SSOH:=0.0; 


TH: =93'001 70. {today's date} 
TSDRS: =0.08; {time between SDR's in qtrs} 
VO15R:=850.00; {mark code 1 and 2 order costs) 


V016:=850.00; 


154 


VO22:=U.1]; {min rrsk} 





Vos%:=0.0; 

V041R:=850.00; {low value annua! demand order cost) 
wO042Rh:=1920.00; {negotiated procurement ordei cost} 
V043R:=1790.00; {advertised procurement order costs) 
V044:=8000.00; {max unpriced order cost} 

VI01A:=0.07; drscRate:=VIOIA; {procurement interest rate} 
hao 22=0.35; {max risk) 

V1034:=1000.00; (shortage cost} 

vlogs =0.1; (repair time preference rate) 

¥295:=1.0; {reorder level constraint) 


LI“LT:=0.0; LILY:=0.0; PCR3:=0.0; O1B:=0.0; Q2B:=0.0; RMNAST:=0.0; SER:=0.0; 
FDR 2=0.0; 

MNOQOOAD: =1.0; {min order qty attrition qtrs demand} 

APSR: =0.0; ARC!I:=0.0; BOO:=0.0; BRLC1:=0.0; BRLDCU:=0.0; BRLO:=0.0; 
Pee 0 0, BRO:=0.0; BOI4A:=0.0; BO19:=0.0; BO19B:=0.0; BO21:=0.0; 
BO21A:=0.0; ERR:=0.0; MONDO:=0.0; OOCI:=0.0; POC:=0.0; PPVBNDO:=0.0; 


PZO:=0.0; RCI:=0.0; RLCI]:=#0.0; RPLCI:=0.0; ROCI:=0.0; VPSR:=0.0; 


Erber ket :=0; 
salvRate:=0.02; 
1atioPLTSTDMU: =0.5; 
infRate:=0.0; 


milEssent:=COO¥8C; 


AAC+ AL+ BOO7A+ BO67G+ CO28+ DRLI+ DO31C+ D120+ DI25N+ ERR1+ FO024+ 


pdg2stri: 
HOD1+ MARLI+ PVPI+ RII+ RO+ YR7POC+ YOO6A+ YOO6B+ EOOIND+ PVUI+ 


FILLER; 


PD82str2:= NumToString(A023B)+ NumToString(BRLDC)+ NumToString(B010)+ 
NumToString(BO11A)+ NumToString(BOl12F)+ NumToString(B019A) + 
NumToStr1ng{(B020)+ NumToSt ring (B023C)+ NumToString(B023D) + 
NumToSt ring (BO23F)+ NumToString(B023H)+ NumToString(BG) + 
NumToSt:ing{(B055)+ NumToString({B0O55A)+ NumToString(B057)+ 
NumToString(B058)+ NumToSt ring (BO58A); 

PD82str3:= NumToString(B06l1)+ NumToString({B070)+ NumToString (B073) + 
NumToSt ring (B093)+ NumToString(B280)+ NumToString(C008C) + 
NumToString(DOPTC)+ NumToString(DTC)+ NumToString(DO025E) + 
NumToSt ring (F009)+ NumToString(HOD)+ NumToString(H0141)+ 
NumToSt ring (H0142)+ NumToString(H0143)+ NumToSt ring (H0144)+ 
NumToSt ring(H0145)+ NumToString(H0146); 

PD82str4:= NumToString(H0147)+ NumToString(H0148)+ NumToString(H0149) + 
NumToSt ring(H01410)+ NumToString(H01411)+ NumToString (H01412)+ 
NumToString(H01413)+ NumToString(H01414)+ NumToString(H01415)+ 
NumToStr ing (H01416)+ NumToString(H01417)+ NumToString(H01418)+ 
NumToStr1ng(H01419)+ NumToSt: ing (H01420)+ NumToSt: ing(ILR)+ 
NumToSti ing{1MECY)+ NumToString(M); 

PD82str5:= NumToString(MOQQAD)+ NumToSt: ing(MSLOAD)+ NumToString (MSLOD) + 
NumToString(NRFIDRT)+ NumToString(OSQ)+ NumToString(PDQ) + 
NumToString(PPV)+ NumToString(QDH)+ NumToString(RF1DRT) + 


NumToString (RIYAYABY)+ NumToString(RSV)+ NumToString(RT) + 


ikSye 





NumToString(SCRi+ NumToSt1ingiSSOH)+ NumToString({TD) + 
NumToString(TSDRS)+ NumToString(V0O15R); 

PO82stro:= NumToString(V01l6)+ NumToString{V022)+ NumToString({V039) + 
NumToSt ring {V041R)+ NumToStiing(V042R)+ NumToString(V043R)+ 
NumToString(V044)+ NumToString(VI01A)+ NumToString{V102)+ 
NumToString(V1034)+ NumToString(V108)+ NumToString(V295) + 
NumToString(LILT)+ NumToString{LILY)+ NumToString(PCR3)+ 
NumToString{(Q1B)+ NumToString{Q2B) ; 

PD82str7:= NumToString{RMNAST)+ NumToString{(SER)+ NumToSt21ing(YDR) + 
NumToSt1ring (MNOQAD) + NumToString{APSR)+ NumToString{ARCI)+ 
NumToSt11ng{BOO)+ NumToString{BRLC1)+ NumToString (BRLDCU) + 
NumToSt ring {BRLO)+ NumToSt1iing({BRPLO)+ NumToSt:z 1ng (BRO) + 
NumToStiing(B014A)+ NumToString{B019)+ NumToSt1i1ng{B019B) + 
NumToString(B021)+ NumToString{BO21A); 

PD&S2str8:= NumToString(ERR)+ NumToString(MONDO)+ NumToSt ring (OQOCI)+ 
NumToString{POC)+ NumToString(PPVBNDO)+ NumToString {PZO) + 
NumToString{RCl)+ NumToString{RLC1l)+ NumToString(RPLCI) + 
NumToString{ROC1)+ NumToString{VPSR); 

assign (outfile, 'pd82in.fil'); 

rewrite (outfile); 

writeln{(outfile,PD82strl, PD82sti2, PD82str3, PD82str4, PD82str15, PD82stro, 

PD82str7, PD82strs); 


close (outfile); 


end; 


procedure PD82Edit (var prbBrkPt: integer; 
var unitPrice,PLT,orde1Cost,holdFrac, 
shortCost,salvRate, numYrsOH, ratioPLTSTDMU, 
numYrsERK,storRate,obsolRate,discRate,infRate, 


milEssent:real); 
var C028 : string[1); 


A023B,B011A,B020,B023C,B023D,B055,B057,B058,B001,B073,C008C,D025E, 


MSLOD,SCR,TD, TSORS,VO15R,V022,V101A,V102,V1034,V295: real; 


Poe2strl: string i2aie 

PD82str2, PDB2str3, PD82str4, PD82str5, PD82stro, PD82str7, 
PD82str8: string[255]; 

editChoice:char; 

done: boolean; 


infile,outfirle:text; 
begin 
{retrieve selected default varrables from file to edit} 
assign (infile,’*pdaS2inJf11")-; 


reset (infile); 


read{infile,PD82strl, PD82str2, PD82str3, PD82str4, PD82str5, PDS2stro, 


Ee) 


PD82sti7, 


close 


(infile); 


PD82st1¥) ; 


CO2Z8:=copy (PD82st1i1,5,1): 


strTemp: =copy (PD82st12,46,15); 


sti Temp: =copy (PD¥2st12,91,15); 


strTemp: 
strTemp: 
sti Temp: 
strTemp: 
strTemp: 
strTemp: 
st1Temp: 
strTemp: 
strTemp: 
strTemp: 
str:Temp: 
st1Temp: 
strTemp: 
strTemp: 
strTemp: 
strTemp: 
strTemp: 


st1Temp: 


unitP11ice:=B055; 


holdFrac:= BO5S7 + VIOLA + SCR; 


xcopy (P82st:2,121,15); 
=copy (PDB2str2,181,15); 
=copy (PD82str2,211,15); 
=copy (PDB2str2,226, 15); 


=copy (PD82st13,1,15); 


eeony (PUS2str3,31,15); 
=copy (PD82st13,76,15); 
mony (PD¥2st:3,121,15); 
=copy (PD82st15,31,15); 
=copy (PD¥2str15,181,15); 
eeopy (PDB2str15,211,15); 
=copy (PD82str5,226,15); 
=copy (PDB2str5,241,15); 
=copy (PD82st16,16,15); 
=copy (PD82st16,106,15); 
=copy (PD8B2stro6,121,15); 
=copy (PD8Zstr6,136,15); 


=copy (PD82st16,1660,15); V295:=StringToReal (St1rTemp) ; 


PLT:= BOLIA; 


done: =FALSE; 


lepeat 

clrsc. ; 

wiiteln{ 
writeln; 
writeln; 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 


writeln; 


writeln {' 


write {' 


editChoice:=upcase(readkey); 


(* 


( 


a3 


_ 
. 


< rm RR GG 


oe ae 


zr Q 39 wm Oo na DoD YS 


Prob Break 
Shelf Life 
Reqn Size 
Unit Price 
Salv. Rate 
Procur LT 
Essential 
Mfg Set-Up : 
Obsol Rate 
Disc Rate 
Time SDRS 
imat) Yrs OH: 


Num Yrs ERR: 


Hit ENTER to 


or letter of 


writeln(edstChoice); 


orderCost :=V015R; 


THIS SCREEN ALLOWS 


* mumYrsERR:8:2,° 


B06] :=StringToReal (StrTemp) ; 


PPL DOBYrKPEZ Be. “ 
aie Be 
ELON Os (any OY 
7B055:822;, : 
,salvRate:8:2, 
PBOLTIAS+2,. 5% 
PASTOR UTS GH SP ee 
,B0S8:8:2Z, : 
,B05728:2, : 
PENG] U2, : 
FLSOURS: 822, * 


»numYrsOH: 8:2, ' 


M. 


N. 


= < C¢ 4 HH 


milEssent :=C0O08C; 


xX. 


Zoe 


BO11A:=Str ingToReal (StrTemp! ; 
BO20:=Sti1ingToReal (StrTemp} ; 
BO23D:=Sti1ingToReal (StrTemp) ; 
BOSS: =Str1 ingToReal (StrTemp}; 
B0S7:=StringToReal (StrTemp) ; 


B058:=Str ingToReal (StirTemp); 


BO73:=StringToReal (StrTemp) ; 
CO08C: =St ringToReal (StrTemp) ; 
DO2SE:=StringToReal (St1rTemp) ; 
MSLOD: =St1ingToReal (StrTemp) ; 
SCR: =StringToReal (StrTemp) ; 
TD: =StringToReal (St1Temp) ; 
TSDRS: =StringToReal (StrTemp) ; 
VO15R:=StringToReal (StrTemp) ; 
V022:=StringToReal (St: Temp) ; 
V101A:=St1 ingToReal (St2rTemp) ; 
V102:=StringToReal (StrTemp) ; 


V1034:=StringToReal (StrTemp) ; 


shortCost;=V1034; 


EDITING OF DEFAULT NIIN INPUT PARAMETERS 


Min Risk 


Max Risk 


Ord Cost 


MSLOD 


QO. Proc Meth 


Shortage 


R/O Low 


R/O Constr: 


Stor Rate 


See a) as 


7 WDPACIET SSG Fd) G 
",V102:8:2): 
SS VOlSRO8: 2); 
Sey yf) 0) 0151: Re) ee 
' DO2ZSE:8:0); 
,V1034:8:2); 
tf BOZ0G 82); 
DO MPADE Bile PAN ¢ 


Oke (3 24)) s 


Time Pref a OLA ea) 5 
Today DT © AAU UIDR SEALY 1 
PEs Os MU) ~ 7 rat boPETSTOMU: B22); 
Inflation Rate: ',infRate:5:3); 


accept current values '); 


field to change: 


ee A 


LE 


msase editChoice of 


PAU 


Pie i 


| FU 


aha 


vGe 


Uo te 


begin 
writeln; 
write ('Ente: new Probability Break Point: ')- 
PibBrkPt :=Ger_Inteqe: (0,20); 
end; 
begin 
writeln; 
write (‘Enter new Shelf Life code: '); 
readin (C028); 
delete (PD82st11,5,1); 
Insert (C028,PD82st11,5); 
end; 
begin 
wiiteln; 
writeln ('** Information Only - Model assumes requisition size of one. 
HitToCont; 
end; 
begin 
writeln; 
write ('Enter new Unit Price: '); 
B055: =Get_Real (0.0,999999.0); 
delete (PD82str2,181,15); 
insert (NumToString(B055) ,PD82str2,181); 
unitPrice:=B055; 
end; 
begin 
w1iiteln; 
wiite (‘Entel new Salvage Rate, fraction of unlit cost: ‘'); 
salvRate:=Get_Real (0.0,1.0); 
end; 
begin 
writeln; 
write ('Enter new Procurement Leadtime Forecast: '); 
BO11A:=Get_Real(0.0,40.0); 
BO23C:=B011A*B023D; 
Gelete (PD82sti:2,46,15); 
insert (NumToString(BO11A) ,PD82st12,46); 
delete (PD&2str2,106,15); 
insert (NumToString(B023C), PD82st1r2,106); 
PLT: =B011A; 
end; 
begin 
writeln; 
write (‘Enter new Average Item Essentiality: ‘'); 
CO08C: =Get_Real (0.0,999999.0); 
milEssent :=C008C; 
delete (PD82st13,760,15); 


Insert (NumToStiing(C008C) ,PD82st13,76); 


end; 


begin 


INS 


are) Ie) 


writeln; 
write (‘Enter new Manufacture: Set-up Cost: ‘'); 
BO58:=Get_Real (0.0,999999,0); 
deiete (Phyzstr2,226,15); 
insert (NumToSt:ring(B0O58) ,PD&¥2st:r2,226); 
end; 
wie: begin 
wiiteln; 
wiite {‘Ente1: new Obsolescence Rate: '‘'); 
BOS7:=<Get_Real (0.0,999999.0); 
obsol Rate: =B057; 
delete (PD#2str2,211,15); 


inseit (NumToString(B057),PD82str2,211); 


end; 
‘J’ : begin 
wiiteln; 
a te ('Enter: new Dirscount Rate: ‘); 
B :=Get_Real (0.0,999999.0); 


delete (PD82st13,1,15); 
insert (NumToString(B0ol),PD82str3,1); 
end; 
*K' 3 begin 
writeln; 
write ('Enter new Time Between SDRs: '‘'); 
TSDRS: =Get_Real (0.0,999999.0); 
delete (PD82str5,2260,15); 
insert (NumToString(TSDRS),PD82str5,226); 
end; 
‘L' : begin 
writeln; 
write (‘Enter number of years demand of initral inventory: ‘}; 
numY1rsOH: =Get_Real (0.0,200.0); 
end; 
‘M’': begin 
writeln; 
write (‘Enter new Minimum Risk: ‘); 
V022:=Get_Real (0.0,1.0); 
delete (PD82stro,16,15); 
insert (NumToString(V022) ,PD82stro,16); 
end; 
‘N' ; begin 
writeln; 
write (‘Ente new Maximum Risk: '); 
V102: =Get_Real(0.0,1.0); 
delete (PD82stro6,121,15); 
insert (NumToString(V102) ,PD82stro,121); 
end; 
*O" : begin 
writeln; 
write ('Enter new Mark I/1I Order Cost: '); 


V015R: =Get_Real (0.0, 999999. 0); 


Log 


orderCost:=VOI15R;: 
delete (PD82str5,241,15); 


insert (NumToString(V015R),PD82str5,241);; 


end; 
dfsi0  g  Jeyeta) in 
writeln; 
write ('Enter new Max Number of Ouarters Safety Levei Demand: ‘); 
MSLOD: =Get_Real (0.0,999999, 0) ; 
delete (PD82str5,31,15); 
insert (NumToString(MSLOD), PD¥2str5,31): 
end; 
HiGM & joyero| \iq' 
writeln; 
write (‘Enter new Procurement Method: '); 
NO2SE: =Get_Real (0.0,999999, 0); 
delete (PD82str3,121,15); 
insert (NumToString(DO2S5E) ,PD82str3,121): 
end; 
‘R' : begin 
writeln; 
write (‘Enter new Procurement Shortage Cost: ‘): 
V1034:=Get_Real (0.0,999999, 0); 
shortCost:=V1034; 
delete (PD82stro,1360,15); 
insert (NumToSt:1ng(V1034) ,PD82stro,136); 
end; 
°‘S* : begin 
writeln; 
write (‘Enter new System Reorder Level Low Limit Qty: '); 
B020:=Get_Real (0.0,999999, 0); 
delete (PD82st12,91,15); 
insert (NumToString(B020) ,PD82str2,91); 
end; 
‘yo ee begun 
writeln; 
write ('Enter new Reorder: Level Constraint Rate: '); 
V295:=Get_Real (0.0,999999, 0); 
delete (PD82stro,160,15); 
insert (NumToString(V295),PD82str6, 166); 
end; 
‘U' + begin 
writeln; 
write (‘Enter new Storage Cost Rate: '); 
SCR: =Get_Real (0.0,99999.0); 
storRate:=SCR; 
delete (PD82st15,181,15); 
insert (NumToSt: r1ng(SCR),PD82str5,181); 
end; 
‘V's begin 
writeln; 


write (‘Enter new Time Preference Rate: '); 


160 








‘Ww’ 


UNG 


V101A:=Get_Real (0.0,99999. 
discRate:=VIOIA; 
delete (PD82stro,106,15); 
insert (NumToString{V101A) ,PD82str6, 106) 
end; 
begin 
writeln; 
write (‘Enter Today''s Date (YYJJJ): '); 
TD: =Get_Real (0.0, 99999.0); 
delete (PD82st15,211,15); 
insert (NumToSt1ring(TD) ,PD82st15,211); 
end; 
begin 
writeln; 
write (‘Enter PLT sigma to mu ratio: ‘'); 
rat 1OPLTSTDMU: =Get_Real (0.0,10.0); 
end; 
begin 
writeln; 
write (‘Enter number of years of economic retention: 
numYrsERR:=Get_Real (0.0,numYrsOH) ; 
end; 
begin 
writeln; 
write (‘Enter current inflatron rate: ‘'); 
infRate;:=Get_Real (0.0,1.0); 


end; 


chr(13): done: =TRUE 


end; 


until done=TRUE; 


holdFrac:=B057 + VIO1A + SCR; 


aeeran (outfile, 'pad82in.fil'); 


rewrite (out file); 


writeln(outfile,PD82str1, PD82str2, PD&82str3, PD&2str4, PD82strs5, 


PD82str7, PD82str8); 


close (outfile); 


else; 


end; 


procedure InitPD86File; 


var infile, outfile:text; 


PD82str1: 
PDOS2str2, 


PD82st1¥: 


PD8ostrl: 


PHSesti2, 


string[24]; 
PD82sti3, PD82str4, PDB2str5, PD82stro, PD82str7, 


string[255]; 


string[24]; 


PD8ostir3, PD¥estr4, PD¥ostr5, PDBostroe, PD8bstr7, 


Lisi 


Ns 


PD82stro, 


PD8ostr#: string(255]; 


PH8ostr9: string[e0}; 


C003,CO001W: string([2]; 
CO01B,LASTIN,COOIT! ,CO01T2,RPRIN,ONEWAY: chaz; 
FILLER: string[5]; 


D046D:string[9]; {NIIN) 


B0O11A,B073,FMLTCNT, FMLYEXP, FMLYGRS, FMLYMNM, FMLYSYSORI), FMLYSYSRO, 
FMLYOPAST, FMLYPLT, FMLYRPRSRV, FMLYRTAT, FMLYROSIZ, FSOPPR1,FSOPPR2,FSOPPR3, 
FSOPPR4 , FSOPPRS, FSOPPR6, FSOPPR7 , FSOPPR8 , FSOPPR9, FSOPPR10,FSOPPR11, 
PSOPPR12,FSOPPR13,FSOPPR14,FSOPPR15, FSOPPR16, FSOPPR17, FSOPPR18,FSOPPR19, 
FSQPPR20, FSOPPR21 , FSOPPR22, FSOPPR23, FSOPPR24, FSOPPR25,FSOPPR2b6, 
PSOPPR27,FPSOPPR2¥, FSOPPR29, FSOPPR30, FSOPPR31,FSOPPR32, FWO, B023D, HRZNLNGTH, 
MEANNONZR ,B001B,B019A,B019B,B019C, B021,B019,B021A,OPAST, PLTPPR, BO12F, PPV, 
PPVO, BRLDCU, F009, B012E,RSV,SOPPR1,SOPPR2,SOPPR3,SOPPR4,SOPPR5,SOPPR6, 
SOPPR7 ,SOPPR8, SOPPR9, SOPPR10,SOPPR11,SOPPR12,SOPPR13,SOPPR14,SOPPR15, 
SOPPR16,SOPPR17,SOPPR18,SOPPR19,SOPPR20, SOPPR21,SOPPR22,SOPPR23, SOPPR24, 
SOPPR25, SOPPR26,SOPPR27,SOPPR28, SOPPR29, SOPPR30,SOPPR31,SOPPR32, 
SYSBO,SYSRCR,A023B,TRPR,TSDRS,BO055, F007, Z0BS, EXPDEFRS, EXPDEFRSR, 
EXPDEFSDR, FEXPDEFRS, FEXPDEFSDR , PROJADDBO, PROJADDVRBL, PROJSMAVRBL, 
PROJSSADDBO, PROJSSADD), PROJSSSMA , ROSHRTRND, ROSHRTYR, VLBUYS, VRBLHRSR, 


VRBLHRSO, UNITSHRTP, UNI TSSHRTR: real; 


begin 
assign (infrle,‘md82out. fil’); 
reset (infile); 
read(infile,PD82strl, PD82str2, PD82str3, PDS2str4, PD82str5, PD82stro, 
PDS2str7, PD82strs); 
close {infile); 
C003: =* 1H". 
COO1LBS => 
LASTING = "3.- 
DO46D:='000000000'; {NIIN} 
COOlT1 =="; 
SCOOT 2 =." 
CO0O1W: =' a. 
RPRIN=="N"; 
ONEWAY:='N'; 
FILLER: =‘ ee 
strTemp: =copy (PD82str2,46,15); BOLIA:=StringToReal (StrTemp) ; 
strTemp: =copy (PD82str3,31,15}); BO73:=StringToReal (StrTemp) ; 
PMLTCNT: =0.0; FMLYEXP: =0.0; FMLYGRS: =0.0; FMLYMNM: =0.0; FMLYSYSORD: =0.0; 


FPMLYSYSRO: =0.0; PMLYOPAST: =0.0; FMLYPLT: =0.0; FMLYRPRSRV: =0. 0; FMLYRTAT: 


il} 

© 

© 
” 


FMLYROSIZ:=0.0; FSOPPR1:=0.0; FSOPPR2: =0.0; FSOPPR3 : =0.0; FSOPPR4:=0.0; 
FSOPPRS:=0.0; FSOPPRO:=0. 0; FSOPPR7:=0.0; FSOPPR8:=0. 0; FSOPPR9: =0.0; 
FSOPPR10:=0.0; FSOPPR11:=0.0; FSOPPR12:=0.0; FSQPPR13:=0.0; FSQPPR14:=0.0; 
FSOPPR15:=0.0; FSOPPR16:=0.0; FSOPPR17:=0.0; FSOPPR18:=0.0; FSOPPR19:=0.0; 


PSOPPR20:=0.0; FSOPPR21:=0.0; FSOPPR22: =0.0; FSOPPR23:=0. 0; FSOPPR24:=0.0; 


EGZ 


FSOPPR25: =0.0; FSOPPR260:=0.0; FSOPPR27: =<0.0; FSQPPR28: =0.U0; FSOPPR29: =! 


FSOPPR30:=0.0;FSOPPR31:=0.0; FSOPPR32:=0.0; FWO:=0.0; 


strTemp:=copy (PD8$2st12,121,15); BO23D:=StringToReal {StrTemp) ; 


HRZNLNGTH: =0.0; MEANNONZR: =0.0;B001B:=0.0; 


strTemp:=copy{PD82st12,76,15); BO19A:=StringToReal (StrTemp) ; 


BO19B:=0.0;B019C:=0.0; 


stiTemp: =copy (PD82sti7,226.15); BO21:=StringToReal (St1Temp) ; 


strTemp: =copy (PD82st17,196,15); BO19:=StringToReal (StrTemp) ; 


Beery. =0.0;OPAST: =0.0;PLTPPR: =0.0;BO12F:=0.0; 


strTemp 
PPVO: =0 
st1Temp 
F004; =0 
RSV: =0. 
SOPPRI1: 
SOPPR3: 
SOPPR®: 
SOPPR14 
SOPPR19 
SOPPR24 


SOPPR29 


:=ecopy {(PD82str15,91,15); PPV: =StringToReal (St: Temp); 


0 WE 


’ 


:=copy (PD82st17,121,15); BRLDCU:=StringToReal (StrTemp) ; 


0 


UG 


=0. 


=0 
=0 


;BO12E:=0.0; 


0; SOPPR2: =0.0; 
.0;SOPPR4:=0.0;SOPPRS:=0.0;SQOPPR6: =0.0;SOPPR7:=0.0;SOPPR8:=0.0; 


-0;SOPPR10:=0.0;SOPPR11:=0.0;SOFPR12:=0.0;SOPPR13:=0.0; 


fo oOrrR)]5:=0.0;SOPPR]6:=0.0; SOPPR17:=0.0;SOPPR18:=0.0; 


>=0.0;SOPPR20: =0.0;SOPPR21:=0.0; SOPPR22:=0.0; SOPPR23:=0.0; 


:#0.0;SOPPR25:=0.0; SOPPR260:=0.0;SOPPR27:=0.0;SOPPR28:=0.0; 


o=0.0; SOPPR30:=0.0;SOPPR31: =0.0;SOPPR32:=0.0; 


SYSBO: =0.0;SYSRCR:=0.0; 


strTemp: 
strTemp: 
strTemp: 


strTemp: 


FO007:=0 


0 


Scomvirbezstr2, 1,15); A023B:=StringToReal (StrTemp); 
=copy (PD82str5,226,15); TRPR:=StringToReal (StrTemp) ; 
=copy (PD82st1:5,2260,15); TSDRS:=StringToReal (St1Temp) ; 


=copy {PD82str2,181,15); BO55:=StringToReal (StrTemp) ; 


ZOBS =O). O!s 


EXPDEFRS: =0.0; EXPDEFRSR: =0.0; EXPDEFSDR: =0. 0; FEXPDEFRS: =0.0; FEXPDEFSDR: =0.0; 


PROJADDBO: =0.0; PROJADDVRBL: =0.0; PROJSMAVRBL: =0.0; PROJSSADDBO: =0. 0; 


PROJSSADD: =0.0; PROJSSSMA: =0.0;ROSHRTRND: =0. 0; ROSHRTYR: =0.0; VLBUYS: =0.0; 


VRBLHRSR: =0.0; VRBLHRSO: =0.0; UNITSHRTP: =0.0;UNITSSHRTR: 


I) 5 (0)9 


{create PD86 input frie} 


PD8ostr1:=C003+ CO01B+ LASTIN+ DO46D+ COO1IT1+ CO01T2+ COOI1W+ RPRIN+ ONEWAY+ 


FILLER; 


PD8ostxr2:=NumToSt: ing(BO11A) +NumToStr ing (B073) +NumToSt r1ng(FMLTCNT) + 


NumToSt ring(FMLYEXP) ¢+NumToString (FMLYGRS) +NumToString (FMLYMNM) + 
NumToSt ring (FMLYSYSORD) +NumToStr ing (FMLYSYSRO) + 

NumToSt ring (FMLYOPAST) +NumToSt1ring (FMLY PLT) + 

NumToSt ring (FMLYRPRSRV) +NumToSt ring (FMLYRTAT) + 

NumToSt ring (FMLYROSIZ) +NumToString(FSOPPR1) + 


NumToString(FSOPPR2) +NumToString(FSOPPR3) +NumToStr ing (FSOPPR4); 


PDSostr3:=NumToString(FSOPPRS) +NumToString(FSOPPRe) +NumToString(FSOPPR7) + 


NumToSt11ng (FSOPPR8) +NumToString(FSOPPR9) +NumToStr ing (FSOPPR10) + 
NumToSt1ing(FSOPPR11)+NumToString(FSOPPR12) + 
NumToSt11ng (FSOPPR13)+NumToString(FSOPPR14) + 
NumToSt ring(FSOPPR15)+NumToSt ring(FSOPPR16)+ 
NumToString (FSOPPR17) +NumToSt: ing(FSOPPR1 8) + 
NumToString(FSOPPR19) +NumToSt ring (FSOPPR20) + 


NumToSt: ing(FSOPPR21); 


Gs 


PD8ost14:=NumToString{ FSOPPR22) +NumToSt ring {FSOPPR23) + 
NumToSt ring (FSOPPR24)+NumToString{FSOPPR25) + 
NumToSt1 rng (FSOPPR26) +NumToSt ring {FSOPPR27) + 
NumToSt:1 ing {FSOPPR28) +NumToString(FSOPPR29) + 
NumToSti ing (FSOPPR30) +NumToSt: ing (FSQOPPR31)+ 
NumTosSti11r1ng( FSOPPR32) +NumToSti31ng (FWO) + 
NumToSt: ing (B0230)} +NumToSt ring (HRZNLNGTH) + 
NumToStz 1ng {MEANNONZR) +NumToSt ring (B061B)+NumToString({BO019A); 

PD8ostr5:=NumToStr ing (B019B) +NumToString(B019C) +NumToString(BO21)}+ 
NumToSt ring (B019)+NumToString(BO21A)+NumToSt ring {OPAST} + 
NumToSt1 ing ({PLTPPR)+NumToSt1 1ng(BO12F) +NumToStiing (PPV) + 
NumToSt ring (PPVO) +NumToSt1ing (BRLDCU) +NumToStz ing (F009) + 
NumToStiing (B012E)+NumToSt ring (RSV) +NumToString(SOPPR1) + 
NumToSt ring (SOPPR2) +NumToString (SOPPR3); 

PDSostr16:=NumToSt: ing (SOPPR4) +NumToString{SOPPRS) +NumToSt11ng(SOPPR6) + 
NumToSt: ing {SOPPR7) +NumToSt ring (SOPPR#) +NumToSt1ing (SOPPR9) + 
NumToSt:1ng(SOPPR10)+NumToSti ing (SOPPR11)+NumToString(SOPPR12) + 
NumToStz:1ng(SOPPR13)+NumToSt ring (SOPPR14)+NumToString(SOPPR15)+ 
NumToSt ring (SOPPR16) +NumToSt ring (SOPPR17)+NumToSt11rng(SOPPR18) + 
NumToSt ring (SOPPR19) +NumToSt:z 1ng(SOPPR20) ; 

PD86sti7:=NumToSt ring (SOPPR21)+NumToSti1ng (SOPPR22) + 
NumToSt ring (SOPPR23) +NumToSt ring (SOPPR24) +NumToSt ring (SOPPR25) + 
NumToSt ring (SOPPR26) +NumToSt ring (SOPPR27) +NumToSt ring (SOPPR28) + 
NumToSt1ing (SOPPR29) +NumToString (SOPPR30) +NumToSt ring (SOPPR31)})+ 
NumToString(SOPPR32) +NumToSt1ing(SYSBO) +NumToSt1ing(SYSRCR) + 
NumToString(A023B) +NumToSt ring (TRPR) +NumToString{TSDRS) ; 

PD8ost2r8:=NumToString(B055)+NumToString{F007)+ 
NumToString (ZOBS) +NumToSt ring (EXPDEFRS) +NumToSt ring (EXPDEFRSR) + 
NumToSt ring (EXPDEFSDR) +NumToSt ring (FEXPDEFRS) + 
NumToSt ring (FEXPDEFSDR) +NumToSt ring ( PROJADDBO) + 
NumToSt ring {PROJADDVRBL) +NumToSt ring (PROJSMAVRBL) + 
NumToSt ring (PROJSSADDBO) +NumToSt ring (PROJSSADD) + 
NumToSt ring {PROJSSSMA) +NumToSt ring {ROSHRTRND) + 
NumToSt ring (ROSHRTYR) +NumToSt ring {VLBUYS) ; 

PD86stzr9:=NumToSt ring (VRBLHRSR) +NumToSt ring (VRBLHRSQ) + 


NumToSt ring (UNITSHRTP) +NumToSt ring (UNITSSHRTR) ; 


assign {outfile, ‘pd8oin.fil'); 
rewrite (outfile); 
writeln(outfile,PD8ostrl, PD8ostr2, PD8ostr3, PD8ostr4, PD8ostr15, PD8ostro, 
PD8ost17, PD8b6str8, PD8ostr9); 
close {outfile); 


end; 


End. {unit pdunit} 
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unit PQueue; 


inte) face 
const MAXPQUEUESIZE=300; 


type dataRkecoid = record 
Oty:intege.; 
Week: integer; 
end; 
HeapArrayType=array [1..MAXPQUEUESIZE] wf datarecoid; 
PriorityQueueType = record 
heapSize:integei ; 
heapA1rray:HeapAri ayType 


end; 


{must be called before the priority queue is first used} 
{also resets the priority queue so it is empty} 


procedure InitralizePriorityQueue (var pQueue: PriorityQueueType) ; 


{error 1f called when it already has MAXPQUEUESIZE elements} 


procedure Insert PriorityQueue (var pQueue:PriorityQueueType; data:datarecord) : 


{returns the element with the largest value} 
{error if no elements in the priority queue} 
function Cur1rWeek (pQueue: PriorityQueueType) :integer; 


function CurrOty (pQueue: PriorityQueueType)}) :integer; 

{removes and returns the element with the largest value} 

{error if no elements in the priority queue} 

function ExtractOty {var pQueue: PriorityQueueType) : integer; 

function ExtractWeek {var pQueue:PriorityQueueType) : integer; 

function EmptyPriorityQueue (pQueue:PriorityQueueType) : boolean; 

function SizePriorityQueue (pQueue: PriorityQueueType) :integer; 
implementation 

{error if the binary trees that are children of the index do not satisfy the 
heap property) 


procedure Heapify (var pQueue:PriorityQueueType; i:integer); 


var left,right,smallest:integer; 


tempVar:dataRecord; 


begin 


with pOueue do begin 


11g. 5) 





eile fats SSO A )) sp Ih 
smallest:=1; 
1f (left -= heapSize) then begin 
if (heapArray [left] .Week » heapArray{1].Week) then begin 
smallest:sleft 
end 
end; 
1f (right<=heapSize) then begin 
rf ({heapArray[right].Week ~ heapArray[smallest].Week) then begin 
smallest:=:rght 
end 
end; 
rf smallest ~ -: then begin 
tempVar:=heapArray [1]; 
heapArray[1]:=heapAriay[smallest]; 
heapArray [smallest] :=tempVar; 
Heapify (pQueue,smallest) 
end 
end {with] 


end; {procedure} 


{removes and returns the element with the largest value] 
ferro: 1f no elements in the priority queue} 


function HeapExtractWeek (var pQueue: PriorityQueueType) : integer; 


begin 
with pQueue do begrn 
HeapExt1i actWeek: =heapArray [1] .Week; 
heapArray[1]):=heapArray [heapSize]; 
heapSrze:=heapS1ze-1; 
Heapify {pQueue, 1) 
end {with) 


end; (procedure) 


{removes and 1eturns the element with the largest value} 
{error if no elements in the priorrty queue) 


function HeapExtractQty (var pQueue:PriorityQueueType) : integer; 


begin 
with pQueue do begin 
HeapExtiactQty: =heapArray[1]).Qty; 
heapArvay [1] :=heapArray [heapSize]; 
heapSt ze: =heapSirze-1; 
Heapify (pQueue,1) 
end {with} 


end; {procedure} 


{erior if called when it already has MAXPQUEUESIZE elements] 


166 





proceduie HeapiInsert (vat pQueue: PriorityQueueType; data:datarecord) ; 


vat index, parent:integer; 


done:boolean; 


begin 
with pQueue do begin 
done:=false; 
heapS1ze: =sheapSize+1; 
index:=heapsize; 
parent:=1ndeyx div 2; 
rf parent=0 then begin 
done: =TRUE 
end else tf (heapArray [parent ].Week «= data.Week) then begin 
done: =TRUE 
end; 
while {index - 1) and (not done) do begin 
heapArray [r1ndex] : =heapArray [parent }; 
index:=parent; 
parent:=index div 2; 
if parent=0 then begin 
done: =TRUE 
end else if (heapAiray[parent] .Week «<= data.Week) then begin 
done: =TRUE 
end 
end; {while} 
heapArray [tndex)]:=data 
end {with} 


end; {procedure} 
procedure InitializePriorityQueue (var pQueue: PrrorityQueueType) ; 
Val index: integer; 
begin 
pQueue.heapSize:=0 
end; {procedure} 
procedure InsertPriorityQueue (var pQueue:PriorttyQueueType; data:dataRecord); 
begin 
HeapInsert (pOQueue, data) 
end; {procedure} 
function Cur11Week (pQueue: PriorityQueueType) :integer; 


begin 


CurrWeek: =pQueue.heapArray [1] .Week; 


iow 


end; {function} 
function CurrOQty (pOueue: PriorityQueueType) :integer; 


begin 
CurrOty:=pQueue.heapArray [1] .Oty; 


enci; {function} 


function ExtractOty {var pOQueue: PrirorrtyQueueType) : integer; 


begin 
ExtractOty:=HeapExtractOty (pQueue}) 


end; { function) 


function ExtractWeek (var pQueue: PriorityQueueType) : integer: 


begin 
ExtractWeek:=HeapExtractWeek (pQueue) 
end; { function} 
function EmptyPriorityQueue (pQueue: PriorityQueueType) : boolean; 
begin 
EmptyPriorrtyQueue: =pQueue.heapS1ize=0 
end: {function} 
function SizePriorirtyQueue (pQueue: PriorityQueueType) : integer; 
begin 


SizePriorityQueue: =pQueue.heapSize 


end; { function) 


end. {unit PQueue} 
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NOTE: The veritcal line in each graph indicates the reorder point 
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